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AYGEN 


High performance, low cost, and relative 
ease of storage and handling make liquid oxygen 
not only the best oxidizer available but guarantee 
it a prominent place among high performance 
oxidizers in the future. 


Liquid oxygen has certain advantages when 
comparing oxidizers: 


1. A potential supply unlimited 

2. Production possible at source 

3. Inexpensive 

4. Non-toxic, non-corrosive, smoke-free 
combustion products 

5. Easily ignited and smooth burning 

6. Few limitations on component ma- 
terials of construction 


7. Not affected appreciably by ambient 
temperatures 


Arthur D. Little recently designed and con- 


ARTHUR OD. LITTLE, 


Mechanical Division 


a 30 MEMORIAL DRIVE + CAMBRIDGE, MASS. 


structed an air-transportable plant for the sepa- 
ration of oxygen from air. This plant was designed 
to produce 10 tons of liquid oxygen per day. 


Experience in the technology of producing, 
handling and storing liquefied gases aided the 
Mechanical Division of Arthur D. Little in the design 
of this plant where minimum weight, small size 
and purity of product were the prime objectives. 
The Mechanical Division has further used its experi- 
ence with liquid gases to design and produce the 
ADL Helium Refrigerator which is for the storage 
of liquid gases to reduce evaporation loss to zero. 


Although liquid oxygen is produced in war- 
time chiefly for use as a propellant it may also be 
used for welding operations and breathing aids 
when converted to gas by the ADL Liquid-Oxygen 
Pump. 


Send for catalog R-2, describing these and 


other products and services of the Mechani- 
cal Division. 
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ROCKET ENGINEERS 


RESEARCH DEVELOPMENT + DESIGN « TEST 
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A long range program of research 
and development in guided missiles 
“~~ has created unlimited opportunities 
in all phases of rocket engineering. 


Engineers with advanced degrees are needed for 
positions in Combustion Research and Physical 
Chemistry. 


Engineers with or without advanced degrees are 
needed as: 


RESEARCH ENGINEERS... for studies in heat 
transfer and Thermodynamics 


DESIGN ENGINEERS... for design phases of © 
liquid rocket power plants, thrust chambers, 
gas turbine pumps 


FIELD ENGINEERS... for coordination of activi- 
ties at field test sites 


TEST ENGINEERS ... for development and pro- 
duction testing of liquid rocket power 
a plants and their components 


| al COMPLETE ROCKET TESTING FACILITIES 


Openings also for Design Draftsmen and Technicians 
Send complete resume to: Manager, Engineering Personnel 


P. O. BOX 1, BUFFALO 5, NEW YORK 


Z 
/ 
/ 0. /~ 
As 
O 
og ify, 
— 
| < 
5 
=, 
| 
= 
| 
< 
| 
> 
| 
5 


To Engineers 
with Families: 


Housing conditions are ex- 
cellent in the Los Angeles 
area. More than 40,000 ren- 
tal units are available in the 
Los Angeles area. Huge 
tracts for home ownership 
are under construction now. 
Thousands of homes have 
been built since the last war. 
Lockheed counselors help 
you get settled. 


CALIFORNIA HAS 
AISED ENGINEERS 
SALARIES 


a The step up to Aircraft Engineering isn’t 
as steep as you might expect. Aircraft expe- 
rience isn’t necessary. Lockheed takes your 
knowledge of engineering principles, your expe- 
rience in other engineering fields, your aptitude, 
and adapts them to aircraft work. You learn 
to work with closer tolerances, you 
become more weight conscious. What’s 
more, Lockheed trains you at full pay. 
You learn by doing—in Lockheed’s 


Lockheed also offers 
these extra benefits: 


Generous Travel allowances * 
Outstanding Retirement Plan 
* Vacations with pay * Low 
cost group life, health, acci- 
dent insurance « Sick Leave 
with pay * Credit Union, for 
savings and low-cost financ- 
ing ¢ Employees’ Recreation 
Clubs ¢ Regular performance 
reviews, to give you every op- 
portunity for promotion On- 
the-job training or special 
courses of instruction when 
needed. 


on-the-job-training program. When neces- 
sary, you attend Lockheed classes. It depends 
on your background and the job you are 
assigned. But always, you learn at full pay. 


gineers 
d the good news: 


Lockheed Aircraft Corporation has raised 
engineers’ salaries recently, in recognition of the 
importance and excellence of their work on 
both military and commercial aircraft. 


The substantial increases make Lockheed engineers 
among the highest paid in the aircraft industry. 


Engineers who join the Lockheed staff 
receive the benefits of these pay boosts. 


In addition to increased salaries, 
Lockheed also offers engineers: 


1. A “‘bonus” every day in better living —just 
because you live in Southern California, 
in an area where the climate is beyond compare. 


2. A future that offers both security and advance- 
ment, helping create planes for defense, 
planes for the world’s airlines in Lockheed’s 
long-range development program. 


3. Better personal working conditions among 
men who have built a reputation 

for leadership in aviation. 


Send today for illustrated brochure describing 
life and work at Lockheed in Southern Cali- 
fornia. Use this handy coupon. 


ENGINEER TRAINING PROGRAM 
Mr. M. V. Mattson, Employment Mgr. Dept. RT-4 


Lakh 
nd AIRCRAFT CORPORATION 


Dear Sir: Please send me your « 
brochure describing life and work at Lockheed. 


Burbank, California 


My Name 


My Field of Engineering 


My Street Address 


My City and State 
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KOLLSMAN 


for products of precision 
and dependability 


Aircraft Instruments and Controls Miniature 
AC Motors for Indicating and Remote Control 
Applications Optical Parts and Optical Devices 
Radio Communications and Navigation Equipment 


Today, Kollsman works without pause toward the 
fulfillment of America’s defense needs. And to our 
nation’s research scientists, the skill, ingenuity and 
creative drive of Kollsman Research Laboratories 
are available for the solution of instrumentation and 
control problems. 
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Fluctuations in a Spray Formed by Two Impinging Jets 


MARCUS F. HEIDMANN' and JACK C. HUMPHREY! 


An investigation of the spray formed by two impinging 
jets of water indicated that upon impingement a ruffled 
sheet of liquid was formed which disintegrated inter- 
mittently to form groups of drops which appear as waves 
propagating from the point of impingement. The effect 
of various injector parameters on the frequency of wave 
formation and on the spray pattern was investigated be- 
cause of its possible relation to combustion instability. 


Introduction | 


OMBUSTION instability has been encountered 

during various phases of development of liquid-_ 
propellant rocket engines. This combustion insta-— 
bility is characterized by sustained oscillations in 
combustion chamber pressure over a wide frequency 
range. This problem has been the subject of a number 
of investigations (1, 2, 3).2. The common concept of 
the analytical treatment of combustion oscillations 
is that the oscillations occur within a feedback loop. 
In general, it is assumed that oscillations in com- 
bustion chamber pressure result from a coupling be- 
tween chamber pressure and energy addition to the 
chamber. In other words, a change in chamber pres- 
sure causes a change in energy addition which in turn 
causes an additional change in chamber pressure. The 
criteria for instability in such a system are that energy 
addition is in phase with the pressure change, and that 
the energy addition results in a pressure change of 
sufficient magnitude to overcome damping factors in 
the system. Such criteria necessitate a change in 
energy addition proportionately larger than the change 
in pressure, or that the gain factor in the feedback loop 
is greater than unity. 

Low-frequency oscillations (in the range from approxi- 
mately 20 to 200 cycles) appear to be a resonance in 
a feedback loop where propellant feed system and rocket 
combustion chamber are dynamically coupled. High- 


Presented at the Annual Convention of the AMERICAN ROCKET 
Socrety, Atlantic City, N. J., November 30, 1951. 

' Aeronautical Research Scientist. 

2 Numbers in parentheses refer to the References on page 167. 
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frequency oscillations (>200 cps) appear to be a feed- 
back cycle within the combustion chamber itself, and 
hence is closely associated with propellant atomization, 
vaporization, burning rate, and related factors. Be- 
cause of the role of the injector in influencing the burn- 
ing, a study of the hydraulic characteristics of the in- 
jector as a possible source of high-frequency instability 
has been undertaken at the NACA Lewis Laboratory. 
The results of the initial investigation where spray 
fluctuations were observed in two impinging jets of 
water (4) are presented herein. 


Impinging Jet Spray Characteristics 
Spray Fluctuations 


In the apparatus for study of the spray characteris- 
tics of two impinging jets, water was forced by regu- 
lated gas pressure through two precision-bore glass tubes 
with a L/D of approximately 75 and arranged so that 
angle of impingement and stream length could be con- 
veniently varied. Several sets of glass tubes were 
used to vary stream diameter. In this manner the 
spray of the impinging jets as a function of jet velocity, 
jet diameter, impingement angle, and distance of the 
two jets to the point of impingement was observed. 
Microflash photographs of two impinging jets of water 
(Fig. 1) show that upon impingement a ruffled sheet of 
liquid is formed perpendicular to the plane of the two 
jets. The liquid sheet disintegrates intermittently, 
forming groups of drops which appear as waves propa- 
gating from the point of impingement. This wave 
propagation is better shown in Fig. 2 which shows a 
series of high-speed motion pictures of the same jets. 

This phenomenon of wave propagation appears to 
be conducive to instability and may be an additional 
gain factor in a feedback loop. The photographs in- 
dicate that the flow rate in the spray may vary from a 
minimum of zero between waves to a maximum of a 
possible 200 per cent of average. If the fluctuations 
could be induced by small variation in chamber pres- 
sure or propellant flow, they would exhibit themselves 
as a very large gain factor in the feedback loop, which 
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may be sufficient to cause sustained oscillations at 
high frequencies. 
In view of the probable relation of spray fluctua- 
tions to combustion instability, a further investigation 
of the characteristics of the fluctuation was made. 
_ The fluctuations in a spray formed by two impinging 
jets of water were studied to determine the nature of 
the fluctuations and to evaluate the effect. of injector 
‘parameters on the fluctuations. The frequency of 
wave formation in the spray was determined by allow- 
: ing the wave to interrupt a beam of light focused on a 
photoelectric cell, and the resulting signal was analyzed 
by using either a frequency meter or recording oscillo- 
scope. The measurements were made at a point where 
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FIG. | MICROFLASH PHOTOGRAPHS OF FORMATION OF SPRAY BY 
TWO IMPINGING JETS OF WATER 
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FIG. 2. SUCCESSIVE FRAMES FROM HIGH-SPEED MOTION PICTURES 

OF FORMATION OF SPRAY FROM TWO IMPINGING JETS OF WATER 
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FIG. 3° TYPICAL WAVE AMPLITUDE AND FREQUENCY CHARAC- 
TERISTICS FOR FLUCTUATIONS RESULTING FROM IMPINGEMENT OF 
TWO JETS OF WATER 
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complete disintegration of the spray occurred, which was 
generally one inch from the point of impingement. 

It will be noted in Figs. 1 and 2 that the groups of 
drops propagating from the point of impingement are 
of various sizes with irregular spacing between groups. 
A quantitative evaluation was made of wave intensity, 
the quantity of liquid per unit volume, and the fre- 
quency from an oscilloscope trace obtained for a typical 
spray. The evaluation was made assuming propor- 
tionality between oscilloscope deflection and quantity 
of liquid per unit volume intercepting the light beam 
directed at the photoelectric cell. The results of this 
evaluation are shown in Fig. 3, which shows frequency 
and intensity on an accumulative basis. For example, 
the frequency at an intensity of 3 is 450 eps, where 450 
is the frequency based on the number of waves of in- 
tensity 3, plus all those larger than 3. The intensity 
scale used is linear, but the magnitude is arbitrary. 
Fig. 3 indicates that there is an abundance of low- 
intensity waves. Approximately one half the waves 
have an intensity of !/;. the maximum intensity. Al- 
though on first impression the fluctuations in the spray 
appear to be random in nature, on the basis of fre- 
quency and intensity the fluctuations do have defined 
characteristics, as illustrated in Fig. 3. Under con- 
stant operating conditions, these characteristics are 
reproducible and the frequency of wave formation over 
a finite time interval is constant. Although Fig. 3 
presents the frequency distribution in the spray, in 
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FIG. 4 EFFECT OF JET VELOCITY AND WAVE FREQUENCY FOR 
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FIG, 5 EFFECT OF IMPINGEMENT ANGLE ON WAVE FREQUENCY FOR 
0.040-IN. DIAM JETS AT VELOCITY OF 40 FPs 
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subsequent references to spray frequency the frequency 
based on all waves formed will be used, which, in the 
case of Fig. 3, is 2580 eps. 


Effect of Injection Parameters on Wave Frequency 
Jet Velocity and Impingement Angle 


The effect of jet velocity on wave frequency over a 
jet velocity range of 20 to 80 fps for impingement angles 
of 60, 80, and 100 deg was determined for 0.025-in.- 
diam jets. The relation between frequency and velocity 
at the various impingement angles is shown in Fig. 4. 
Fig. 4 shows that frequency increases with velocity, 
with the relation approaching direct proportionality. 
The frequency gradient is approximately 30-cps velocity. 
Although no data were taken below a velocity of 20 fps, 
the nature of the curves indicates extrapolation to the 
origin. 

The trend of increasing frequency with decreasing 
impingement angle, indicated in Fig. 4, is more clearly 
In this case, the effect of impinge- 
ment angle on frequency was determined, using 0.040- 
in.-diam jets. The frequency was determined over the 
maximum impingement-angle range of the apparatus, 
including the frequency measured in the nonimping- 
ing jet, at a point where complete disintegration of the 
jet exists. The curve shown in Fig. 5 is a theoretical 
curve of frequency based on the change in resultant 
spray velocity with impingement angle (assuming 
conservation of momentum). Good agreement with the 
theoretical curve was obtained at impingement angles 
of 50 deg and above. The frequencies obtained at 
low impingement angles are reproducible; however, 
some of the deviation from the theoretical curve may 


shown in Fig. 5. 
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be attributed to the indistinct nature of the waves at 
low impingement angles. 

The variations in the spray pattern resulting from 
changes in jet velocity are shown in Fig. 6. The pat- 
terns for velocities of 20, 40, 60, and 80 fps at an im- 
pingement angle of 60 deg are shown for 0.025-in.-diam 
jets. These photographs illustrate that the wave pat- 
tern is not substantially altered by jet velocity. Smaller 
drops, however, are formed at the higher velocities. 

The wave patterns for impingement angles of 30, 60, 
and 90 deg at a jet velocity of 60 fps for 0.025-in.-diam 
jets are shown in Fig. 7. These photographs indicate 
greater definition of wave pattern as well as a shorter 
distance for complete disintegration of the liquid 
sheet as impingement angle increases. 


Jet Diameter 


& 

Frequency measurements were made over a velocity 
range of from 20 to 80 fps and at an impingement angle 
of 80 deg for jet diameters of 0.025, 0.040, and 0.057 in. 
These results are shown in Fig. 8. An increase in fre- 
quency with decreasing diameter at constant velocity 
is indicated; however, the effect of diameter is compara- 
tively small. A diameter range of from 0.025 to 0.057 
in. at constant velocity represents a flow range of 
approximately 5 and resulted in a frequency variation 
of from 200 to 300 cps. Inasmuch as velocity based 
on flow calibration was used in the comparison, the 
small variations in frequency are insufficient to predict 
the effect of diameter beyond that the effect is small in 
comparison to the effect of impingement angle and veloc- 
ity. 

The effect of diameter on the spray pattern is illus- 
trated in the photographs shown in Fig. 9. The spray 
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FIG. 7 EFFECT OF IMPINGEMENT ANGLE ON SPRAY PATTERS. 
2 JET DIAM, 0.025 1n.; JET VELOcITY, 60 FPS 


patterns for the 0.025, 0.040, and 0.057-in.-diam jets 
at a velocity of 50 fps and an impingement angle of 
70 deg are shown. Beyond that, the size of the result- 
ing spray changes with diameter, and only small 
differences are detectable. 
of successive enlargements of the same picture. 


The length of the jets before impingement was varied 
under several conditions to determine the effect of this 
parameter on wave frequency. The effect of changing 
jet length over a range of 6 to 80 diam had a negligible 
effect on spray frequency. The effect of jet length, 
however, is more pronounced in spray pattern. Fig. 10 
portrays the effect on the 0.025-in.-diam jets ata velocity 
of 40 fps and angle of 60 deg. Stream lengths of 10, 
20, 40, and 60 diam are shown. It will be noted that 
dispersion of the spray increases as stream length in- 
creases. At a stream length of 60 diam, the jet has 
become very unstable and almost complete separation 
occurs. Beyond this length it is possible for a fraction 
of liquid from one jet to pass through the path of the 
other jet without interference or impingement. The 
spray patterns perpendicular to those in Fig. 10 are 
not presented because the spray pattern in this plane 


was substantially unaffected. “~~ 


Discussion 


Jet Length 


The results of this initial investigation of fluctua- 
tions in a spray formed by two impinging jets indicate 
that the major effect on the fluctuations was obtained 
by a change in the velocity of the sheet formed by the 
If, in Fig. 4, the jet velocity and sheet veloc- 


two jets. 
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CONDITIONS: 


Fig. 9 gives the appearance 
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FIG. EFFECT OF JET DIAMETER ON RELATION BETWEEN 
FREQUENCY AND JET VELOCITY AT IMPINGEMENT ANGLE OF 
80 DEG 


ity are assumed to be directly proportional for a 
constant impingement angle, the frequency is 
shown to vary proportionately with sheet velocity. 
Also in Fig. 5, as impingement angle is changed, 
the sheet velocity should vary as the cosine of one 
half the impingement angle (based on conservation 
of momentum). 

This relation applied to the experimental effect 
of impingement angle is shown in Fig. 5. The correla- 
tion appears good for impingement angles between 50 
and 100 deg. and the frequency of a nonimpinging jet 
corresponds to the frequency for a 0-deg impingement 
angle. Lack of correlation at small impingement 
angles, however, indicates that the frequency is not 
a unique function of the sheet velocity. One other 
factor apparently affecting the ruffling of the sheet 
and the subsequent frequency is the irregular nature 
of the jets before impingement. This effect is indicated 
by the photographs in Fig. 10. In this case, increasing 
the jet length amplifies the irregularities in the jets, 
and a greater dispersion of the spray is obtained. Other 
factors may also influence the frequency; however, 
a complete theory on the formation of these fluctuations 
in the spray can be obtained only by a continuation of 
the investigation of the parameters affecting the spray 
formation. 

In determining the effect of injector parameters on 
the spray formation, it was found that changes in the 
injector parameters had no effect on the basic mechanism 
of spray formation. Under all conditions, the spray 
was formed from a ruffled sheet of liquid that disinte- 
grated intermittently, and the changes obtained in 
the measured frequency of wave formation with injector 
_ parameters we re comparative ly uniform. 

The frequency of wave formation is important to 
the problem of combustion instability, for it indicates 
that the propellant flow rate in the spray can vary 
with sufficient rapidity to promote high-frequency 
pressure fluctuations, screaming, in combustion. The 
relation of the fluctuations in the spray to combustion 
instability, however, appears to be dependent on the 
problem of inducing fluctuations of a definite frequency 
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FIG. 9 


n the spray. Limited experiments with sound waves 
iave shown that fluctuations in a spray formed by two 
mpinging jets of water cannot be readily induced. 
Sound waves produced by an oscillator and speaker 
und concentrated on the spray gave inconclusive results. 
\t lower frequencies the spray could be disrupted, but 
his appeared to result from the air turbulence caused 
by the high concentration of sound energy. Waves 
were not induced but, rather, the spray or sheet was 
virtually torn apart independent of the wave formation. 
At induced frequencies approaching the natural fre- 
quency in the spray, the sound waves had no apparent 
effect on spray formation. Experiments with liquid 
pressure pulsators and with mechanical vibrators 
have also indicated that fluctuations in the spray are 
not readily induced. 

Preliminary experiments on the effect of viscosity 
and surface tension on spray formation have shown 
that, although spray formation changes considerably 
with higher viscosity, the effect of surface tension is 
small. The effect of viscosity is presently being ex- 
plored to provide data on a theory of spray formation 
and the resulting fluctuations. Such a theory will 
be useful in determining whether fluctuations can be 
induced and their relation to combustion instability. 
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FIG. 10 
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SFFECT OF JET DIAMETER ON SPRAY PATTERN. 
PINGEMENT ANGLE, 70 DEG 


CONDITIONS: JET VELOCITY, 59 FPS; IM- 


Summary of Results 


From the investigation of the fluctuations in the 
formation of a spray from two impinging water jets, 
the following results were obtained: 

1 Upon impingement of two jets, a ruffled sheet 
of liquid was formed, perpendicular to the plane of 
the two jets. The liquid sheet disintegrates intermit- 
tently, forming groups of drops which appear as waves 
propagating from the point of impingement. 

2 The intermittent disintegration of the liquid 
sheet resulted in irregular spacing between waves and 
in variable wave intensity. There was an abundance 
of small waves with the number of waves above a given 
intensity decreasing as the intensity increased. 

3 The frequency of wave formation was constant 
over a finite time interval under constant operating 
conditions. The frequency varied between 1000 and 
4000 cps for the range of test conditions used in this 
investigation. 

4 The spray parameter having the major effect on 
spray frequency was found to be the velocity of the 


sheet formed on impingement. The relation between 


(Continued on page 167) 
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HYDRAULICS RESEARCH LABORATORY 


whe 4 comparison of various types of thrust chamber injec- 
tion systems is outlined. The applicability and relative 
merits of several common injector heads are emphasized. 
The use of dyed-water flow analysis and high-speed 
photography in injector head spray research is described. 
It is shown that dye flow permits distinction of the “‘pro- 
pellants” in the spray and enables the designer to take 
corrective steps if poor distribution is indicated. Methods 
of analyzing droplet size and atomization rates are re- 
viewed and compared. The cause and effect of hydraulic 
“flip”? in orifice jets are discussed, as are the effects of 
cavitation and feedback in the hydraulic system. 


Introduction 
IQUID-FUEL rocket motor systems are character- 
ized by flows of large quantities of liquids at high 
pressures and velocities. This characteristic inevitably 
results in anomalous hydraulic behavior and severe 
design criteria for the various components of the motor 
complex. The propellants must begin to flow from a 
reservoir (tank) through a maze of piping and valves, 
possibly through a cooling jacket, ultimately through 
an injector head, into the environment of the chamber. 
During its travels, the liquid must fight an uphill 
battle (energywise), against pipe, fitting, valve and 
jacket friction, the injector head orifice holes, and, 


Presented at the Annual Convention of the American Rocket 
Society, Atlantic City, N. J., November 30, 1951. 
1 Rocket Research Engineer. ‘Member ARS. 
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finally, the high chamber pressure. It can derive its 
initial energy from pressurized tanks, or it may start 
out at a low energy level and get a boost along the way 
by a pump. 

It is difficult to treat any one element distinctly 
and separately from the rest of the hydraulic circuit; 
however, some distinction is possible if interaction is 
neglected. 

One of the more important components of the rocket 
motor system, the injector head, will be discussed 
briefly, assuming that it is independent of the remainder 


of the hydraulic circuit. 4] 


Basically, the injector head has the functions of intro- 
ducing the propellants into the thrust chamber in the 
most efficient and auspicious manner (i.e., it should 
expedite the transfer of the propellant from the liquid 
to the gaseous phase as quickly as possible). The 
injection can either be random or indiscriminate, leav- 
ing the functions of atomization and mixing to the 

vagaries of the chamber; or, the injection system can 
prepare the propellants for combustion by performing 
these functions, utilizing some of the kinetic energy of 
the moving liquids. 

A brief classification of injection systems, using the 
criteria of mixing and atomization, is as follows: 

Random (Fig. 1): A typical example of this injec- 
tion system is the ‘random showerhead” injector. It 
usually consists of a matrix of holes drilled in the flat 
face of the injector. Oxidizer and fuel are injected 
into the chamber through adjacent orifices. The rate 
and ultimate atomization of the liquid streams depend 
upon the chamber environment, such as the chamber gas 
density, temperature, and injection velocity. The 
mixing depends upon the collision between the droplets 
and the subsequent mingling. 

Generally, random injectors 
ciency with nonhypergolic propellants, probably due 
to their poor mixing and distribution characteristics. 
The mixing may be improved by increasing the number 
of orifices; however, the manifolding design and con- 
struction problems increase disproportionately. Since 
the combustion with this injector tends to be smooth 
with common propellants, further research is warranted 
to improve the combustion efficiency of this head, while 


Injector Head Design 


have shown low effi- 


retaining its stable combustion characteristics. 
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Splash (Fig. 2): The principle of this type of 
injection is primarily one of the impingement of the 
propellant stream or streams upon the walls of the 
chamber or some protuberance of the injector head, 
such as an annular ring or “splash plate.”’ This plate 
serves to break up the stream for rapid atomization. 
The inelastic collision of the stream also abruptly 
changes its velocity and direction, and results in burn- 
ing close to the injector head. This mechanical stream 
breakup permits low injection pressures. 

The German “Enzian” injector head is a classical 
example. This configuration has been tried with vary- 
ing success in this country. It may yield high speci- 
fic impulses with relatively smooth combustion at differ- 
ential injection pressures up to 50 psi. However, the 
combustion pattern tends to be uneven with severe 
localized heat rejection, in extreme cases almost 5 Btu 
sq in./sec. Also, the position and tolerance of the 
splash ring are critical and deviations from design can 


result in wide variations in performance. This condi- 


common nonhypergolic propellants, the well designed 
impinging showerhead can consistently yield specific 
impulses with values 90 per cent of theoretical; the 
value of the heat rejection may be reduced by designing 
for a spray envelope (usually fuel rich) whose mixture 
ratio differs from the rated value. Some evidence exists 
which indicates the injector is sensitive to sonic and 
transverse oscillations which influence the fluid stream 
impingement resulting in rough or unstable combustion. 

This configuration is a close approach to a “‘standard”’ 
injector in the rocket motor industry. 


b. The Hypoid Showerhead 


(Fig. 4): Besides the normal impingement angle, 
the fluid streams are subjected to another degree of 
direction or “twist” by the alignment of the orifice 
holes so as to have not only impingement mixing, but 
also further forced mixing of the propellants, due to the 
hypoid action. The fishtails of mixed liquid impinge 
to form a sheet. Injector heads of this type can give 
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, FIG. | RANDOM SHOWERHEAD FIG. 2 


INJECTION 


SPLASH INJECTION 


tion has been remedied to some extent by knurling or 
serrating the splash ring face in order to break up and 
distribute the spray more efficiently. 

Liquid-Liquid (sometimes known as “forced action 
injection”): This injection method comprises the 


largest segment of current injector practice. Its 
chief virtue lies in its ability to break up and mix the 
! propellant streams, using some of the kinetic energy 
: of the liquids. In general, a well-designed head gives 
! good mixing and distribution at any engine size from 


100-Ib thrust to existing upper limits (1).? 
Some typical examples of injector heads based on 
this principle are: 


a. The Impinging Showerhead 


(Fig. 3): This is a favorite configuration, generally 
consisting of concentric rows of holes for fuel and oxidi- 
zer. It has proved advantageous to bring the propel- 
lant impingement points as close to the injector face 
as possible. This permits maximum utilization of the 
thrust chamber volume and results in smoother com- 
bustion, apparently due to the absence of unreactive 
vapor pockets near the injector face. With presently 


? Numbers in parentheses refer to References on page 138. 
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FIG. IMPINGING (LIQUID- FIG. 4 HYPOID INJECTION 


LIQUID) INJECTION 


a specific impulse at least as good as that obtained with 
“normal” impinging injectors, but with lower propellant 
injection pressures. These perform 
equally well with hypergolic or nonhypergolic propel- 
lants (which is not true for many other types of injec- 
tors). However, they are subject to face burning and 
occasionally give rise to high-frequency oscillations; 
furthermore, the hypoid configuration is not easily 
machined with common shop methods. 


injectors can 


c. Concentric Ring 


(Fig. 5): This isa variation of the impinging shower- 
head, depending for mixing and atomization upon the 
impingement of liquid sheets emitted by annular slots 
in the injector head. This configuration gives excellent 
extraneous-head mixing, since the effect of the impinging 
sheets is that produced by an infinite number of single, 
adjoining liquid streams (2). This contrasts with the 
mixing produced by double jet impingment, which 
may be erratic due to jet misalignment and oscillation. 
However, the considerable difficulty of manufacture 
of the concentric ring injector prevents its wide applica- 
tion. Also, excessive heat transfer rates have been 
experienced with this type of injector. 
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There are obviously many permutations of the above 
configurations. The designer can let his imagination 
revel at will, and produce types such as the like-on- 
like (fuel impinges on fuel, ete.) injector, or injector 
heads with inserts of swirler or common oil-burner 
nozzles (Fig. 6). These latter spray nozzles which are 
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CONCENTRIC RING 
(ANNULAR SLOTS SECTION 


MIXED CYLINDRICAL SHEET 


SPRAY NOZZLE 


FIG. 5 CONCENTRIC RING FIG. 6 SPRAY NOZZLE 


commercially available, have been favorites of design- 
ers since the early days of rocket motor manufacture. 
It seemed simple enough to insert these nozzles into an 
injector-head face according to some empirical geometri- 
cal pattern. The wide angle of spray and good atomiza- 
tion were counted on to give good performance. How- 
ever, this type of injector head is not widely applied. 
This may be due to several reasons, the principal one 
being that the specific impulse is somewhat less than 
for the other types. 

An element of the commercial spray nozzle, the swirler 
or vane insert, which imparts angular momentum to 
the spray, is occasionally inserted into the plenum of 
normal head orifices in order to retain the atomizing 
features of the spray nozzle without sacrificing the 
flexibility and ease of manifolding of the conventional 
injector. Some success has been attained with differ- 
ential injection pressures of 200-300 psi, obviously too 
high for practical applications. 

Forced action, or “stream impinging” injectors will 
probably remain the primary injectors for some time to 
come, for motors up to five tons thrust. The configura- 
tion is flexible enough to permit ready adaptation to 
various sizes and types of thrust chambers. However, 
this type of injector must be accurately drilled, both 
for impingement angle and directness of impact, and for 
smoothness of the holes. Many otherwise well-imping- 
ing showerheads are discarded because of scored and 
burred holes, which affect the stream geometry. 

Premix (Fig. 7): Premixing the propellants before 
introduction into the chamber would seem to be a 
most logical and obvious answer to the ills aceompany- 
ing other injection methods. 

In a typical premix injector nozzle, the propellants 
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are individually injected into a plenum chamber and are 
partially mixed. They then flow through a swirl 
chamber, in which mixing is completed and an angular 
momentum is imparted to the liquid in order to produce 
a conical spray. If a conical spray is not desired, pre- 
mixing is simplified, and is often accomplished by having 
two concentric annular nozzles, with the liquids mixing 
just before entry into the chamber. The V-2 injectors 
were, in a sense, premix types, with the alcohol and 
oxygen mingling in a small ‘‘ante-chamber.”’ 

Since premixing is a liquid-liquid interaction of com- 
bustibles, the flame front can progress upstream to a 
point where the mixture is inadequate to maintain 
burning. This usually happens, and if the burning 
reaches a confined domain, an explosion occurs. There 
are two ways to prevent this: either by increasing the 
injection pressure until the fluid velocity is greater than 
the flame velocity (whatever that may be), or by increas- 
ing the ignition lag, so that the liquid can leave the 
nozzle before burning occurs. Neither alternative is 
desirable or easily accomplished. 

Despite these drawbacks, the theoretical advantages 
of a premix injector are attractive enough to warrant 
further research, expecially for very large thrust cham- 
bers, using nonhypergolic fuels where “standard” 
injectors would be too heavy, cumbersome, and expen- 
sive. 

To summarize, it may be said that no single injector 
type or system is ideal or applicable to all circumstances. 
However, inevitably, greater standardization will occur, 
when more data on the significance of the injector are 
available, and when many of the empiricisms have been 
reduced to scientific foundations. 


Injector Research 


It is difficult to assess the absolute importance of the 
injector head when considering the over-all rocket motor 
system. 

An “‘ideal’’ injector should result in the following 
desirable thrust chamber characteristics: (a) A high 
specific impulse; (b) stable combustion over a wide 
range of chamber pressure; (c) smooth start-up; (d) a 
low rate of heat injection, within the working limits 
of existing coolant flows and materials; (e) an absence 
of streaky or localized combustion; (f) minimum throat 
erosion; (g) efficient operation through wide variations 
of ambient temperatures. 

This ideal head should 
| also possess the following 
features: (a) Ease of manu- 
facture and reproducibility ; 
(b) a minimum of face ero- 
sion; (c)a low injection pres- 
amen sure; (d) adaptability to 

various propellant combina- 
tions: (e) a minimum of 
critical material; (f) insensi- 
tivity to vibrations of vari- 
ous types; (g) steady fluid 
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FIG. 7 PREMIX jets. 
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Little quantitative data are available which directly 
relate the hydraulic injector constants with the above 
parameters. For this reason, experimental research is 
continuing and new methods of evaluation are constantly 
being tried. 

In practice, each new injector head is first tested 
hydraulically, and, if satisfactory, is then checked in a 
carefully controlled firing test. 

Hydraulic behavior, of course, is not necessarily a 
criterion of chamber behavior; however, hydraulic 
studies are at least a beginning and have shown that 
some correlation between test-stand performance is 
The following are some of the injector head 
or orifice spray characteristics which seem to be impor- 
tant, and methods of their analysis. 

Distribution of Propellants: A knowledge of the rela- 
tive locations of oxidizer and fuel streams in the injector 
spray is often very useful. An oxidizer envelope may 
presage excessive chamber erosion or heat rejection. 
Discontinuities or streakiness may result in localized 
erosion or overheating of chamber or nozzle. <A satis- 
factory qualitative analysis of this parameter has been 
evolved by the use of dye injection into the oxidizer or 
fuel water flow lines. Many dyes have been tried, but 
ultimately aniline “Wool Yellow’ and aniline ‘Niagara 
Sky Blue’’’ proved the best, for photographic and visual 
observation. A further aid is a ‘“strobolux’”’ or other 
flashing stroboscopic light which slows down the spray 
velocity enough to permit detailed visual observation. 
A mechanical sampling method which has been a 
standby since the early days of spray studies, is the 
test-tube collector (Fig. 8). A large number of tubes 
are exposed to the spray and then checked for liquid 
level individually, to correlate with the geometrical 


possible. 


FIG. 8 TEST-TUBE SPRAY SAMPLER 


’ Made by National Aniline Div., Allied Chemical Corp., 
Buffalo, N. Y. 


* Equivalent dyes, soluble in water, can be ae 
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location of the tube. This method is cumbersome and 
rather difficult to interpret for large flow rates. 

Mixing of Propellants: Actual tests have shown that 
rapid and efficient mixing favorably influences combus- 
tion efficiency and stability. Good mixing depends 
upon a number of factors, such as the miscibility of the 
liquids, the velocity and area of impact (in impinging 
streams), the size and surface tension of the droplets 
(if the streams are atomized before collision) and con- 
vective recycling of the vaporized fluids in the hot thrust 
chamber. The last factor can only be surmised, con- 
sidering existing data and equipment. The other fac- 
tors and the actual mixing can be evaluated mechani- 
cally with sampling tubes, using water, carbon tetra- 
chloride, or other mixtures, acid-base mixtures for pH 
measurements of the tube contents, or color analysis with 
a colorimeter. In the later case, a small sampling tube 
is traversed across the dyed water spray and a sample 
is taken by momentarily opening and closing a little 
‘ap valve at the mouth of the tube. The colored 
sample is then compared with a standard reference 
Thus a polar plot of mixture 
ratio across the spray pattern may be made. Attempts 
have been made to reduce surface tension and promote 
mixing by the addition of surface active agents. So 


sample in a colorimeter. 


far, results have been inconclusive. 

Velocity: 
have a bearing on the mixing rate; furthermore, it 
influences penetration into the chamber before burning, 
rate of ligament formation in the stream, atomization 
rate, and the direction of the resultant stream (after 
impingement). This latter factor is actually a result 
of the momenta of the two liquids. The velocity can 
be easily computed from standard hydraulic equations 
relating discharge coefficient, weight flow, orifice area, 
and injection pressure. Droplet velocity may also be 
measured by high-speed photography with calibrated 
grids or timing lights. 

Atomization Rate and Droplet Size: If the propellant 
mixture is finely atomized, the total available surface 
area exposed to the combustion process increases thou- 
sands of times. Similarly, a rapid rate of atomization is 
desirable since it disperses the liquids into this desirable 
state with the least delay, and burning can proceed 
efficiently, making the most of the available chamber 


The stream velocity, as stated before, may 


volume. 

Since the surface area of a droplet decreases more 
slowly than its volume, with decreasing diameter, the 
available surface area-to-volume relationship of finely 
atomized droplet sprays is quite favorable. This factor 
expedites the rate of transformation from the liquid-to- 
gas-to-combustion state. On the other hand, the rate 
of vaporization is influenced at once by the high chamber 
pressure, and the convective and radiative heat trans- 
fer to the droplets. Experiments have indicated that a 
high chamber pressure tends to retard the rate of evapor- 
ation more than the high temperature accelerates it. 
Further research is required to substantiate this and to 
develop techniques of direct measurement of the liquid 
evaporation (3). 
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OPTICAL PHOTOMETER FOR DROPLET ANALYSIS OF FINE SPRAYS 
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FIG. 9 SPRAY DROPLET PHOTOMETER DIAGRAM 
OPTICAL” SHADOW" SPRAY ANALYSIS 
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ABOUT 10” FROM SPRAY 


FIG. 10 OPTICAL SHADOW SYSTEM 


FIG. 1l SPRAY SHADOW PHOTOGRAPH 


ss Much work has been done on the observation and 
measurement of spray droplet sizes and the factors 
affecting the rate of formation of drops from liquid 
streams, and the ultimate drop diameter. A favorite 
mechanical method of sampling the droplets has been 
by their impingement upon magnesium-oxide covered 
glass slides; the drops are assumed to leave an impact 
pattern analogous to their size (and distribution). 
The impact cavities, or the actual drops, with some types 
of slides, are then laboriously counted and measured 
with a calibrated microscope (4). 

Another, indirect, method of measuring droplet 
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- tion remains to be established. a 


sizes is an optical technique; a beam of light is passed 
through the spray and is received by a photo detector 
which measures the intensity of this transmitted beam, 
as compared to some reference beam. This is in effect 
a photometer (Fig. 9) which relates light transmission 
loss to droplet diameter and concentration (5). 

The above methods are unsuitable for heavy, coarser 
rocket injector sprays; here photography is a useful 
method, especially high-speed photographs (5 microsec 
or less) which can ‘freeze’ the spray pattern. If this 
is not too dense, the individual droplets can readily be 
counted and measured at A variation of 
normal photography is the ‘‘shadow” method (Fig. 10). 
Light from a pinhole source is paralleled by a concave 
mirror (or convex lens) and sent through the spray. 
A photographic plate may be placed at some distance 
beyond the spray and exposed by a high-speed flash 
light. The resultant image is full size and an exact 
silhouette of the droplet (Fig. 11). 

No definite working relationships have been estab- 
lished for the prediction of spray droplet size and 
behavior in a thrust chamber (or even in the atmos- 
phere). Considering the many existing propellants 
and variations in combustion conditions, such relation- 
ships may be impossible to establish rigorously. 
Factors, such as high injection and chamber pressure, 
low viscosity and surface tension, have been widely 
proved to favorably affect a rapid atomization and small 
droplet formation. A typical relationship established 
by Triebnigg (6) combines some of these factors (Fig. 


12). 
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ae FIG. 12 EQUATION RELATING SPRAY FACTORS 

The relation was determined for a diesel fuel injec- 
tion system, but is typical of the conditions encountered 
in rocket practice. 

It is difficult to measure the above quantities under 
actual working conditions; also the relative importance 
of atomization rate and droplet size in rocket combus- 

Other Hydraulic Factors 

Although the “normal” hydraulic behavior of the 
injector head is readily analyzed and plays a large part 
in rocket motor hydraulic research, other important 
fluid flow phenomena occur in the motor system and 
injector head which are not so easily evaluated. Some 
of these are: 

Hydraulic “Flip” (Figs. 13and14): This condition 
of unstable orifice flow occurs in orifice holes with length/ 
diameter ratios greater than those of knife-edged ori- 
fices, i.e., typical holes in an injector head. 
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INCREASING INJECTION PRESSURE 


FIG. 13 


HYDRAULIC “FLIP” SEQUENCE SHOWING 
CHANGING STREAM CHARACTERISTICS 


PLEXIGLAS ORIFICE, HYDRAULIC “‘FLIP”’ 


INCREASING INJECTION PRESSURE A 


INCREASING INJECTION PRESSURE 


BELOW FLIP TO ABOVE FLIP 
FIG. 14 HYDRAULIC “FLIP”? SEQUENCE 


At definite injection pressures, related to particular 
length/diameter (L/d) ratios, the water jet flowing 
through the orifice suddenly changes its characteristics 
as the injection pressure is increased slightly. Flow 
becomes quite turbulent and then smooth, with a rapid 
uncontrollable increase in pressure drop but no increase 
in the flow rate (Fig. 15). 

This phenomenon has been adequately treated in a 
previous report (7); it suffices to say that the cause 
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ABOVE FLIP TO BELOW FLIP 


seems to be the breaking-away of the stream from the 
orifice walls except at the inlet edge, effectively estab- 
lishing a knife-edged orifice. A reversible flip transi- 
tion time of '/329 sec has been measured with water. 
If this transition time occurred with actual motors 
and if a reversible flip action took place continuously, 
a combustion frequency of approximately 320 cycles 
second might be expected. Actually, frequencies near 
this value have been detected with injectors which had 
a hydraulic flip near the rated flow. 

Impinging Jet Stability: It has been observed that 
impinging fluid jets tend to form a “Christmas tree” 
or “pagoda” pattern of fluid ligaments (Fig. 16). This 
formation seems to be caused by a combination of: 
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FIG. 17 GAS POCKETS IN INJECTOR MANIFOLDS 


change in kinetic energy, surface tension, and air (or 
chamber gas) friction (8). This latter factor does not 
seem too important. Work done at Bell shows no 
change in frequency or number of these ligaments with 
low (vacuum) chamber pressure and 30-psig chamber 
pressures. 

Some partial evidence exists relating 90 cycle pagoda 
frequencies (measured with a photoelectric counter) 
with combustion stability frequencies of 90-100 cycles. 

Cavitation: Cavitation existing in pump-fed systems, 
and its effect is a common problem and need not be 
elaborated here. However, other instances of gas 
entrapment have been discovered in injector head 
manifolds. The gas pockets were detected by inserting 
a transparent back-plate in the manifold and observing 
the flow under stroboscopic illumination (Fig. 17). 

It was assumed that these pockets could act as elastic 
cavities which would induce pulsations in the orifice 


instability. 


jet and act as a source of burning instability. Attempts 
have been made to tap injector manifolds, and cooling 
jackets, in order to “bleed” off these trapped gas cavi- 
ties and eliminate them as a possible source of burning 
instability. While some success has been obtained, it 
is apparent that further research is required in order to 
establish this form of cavitation as a definite source 
Also, further research may show other, 
hitherto unsuspected examples of gas pockets in the 
propellant feed system. 


The author wishes to acknowledge the collaboration 
of Mr. F. Stratton in the photographic techniques, and 
the assistance of Mr. C. Kurland in the Hydraulics 
Laboratory, and invaluable contributions by various 
members of the Rocket Power Plant Division. 
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Transfer Functions of Rocket Nozzles 


- Daniel and Florence Guggenheim Jet Propulsion Center, al 
California Institute of Technology, Pasadena, Calif. 


The 


ENTRANCE THROAT EXIT 


The transfer function is defined as the fractional oscil- 
lating mass flow rate divided by the fractional sinusoidal 
pressure oscillation in the rocket combustion chamber. 
Chis is calculated as a function of the frequency of oscilla- 
tion. For very small frequencies, the transfer function is 
approximately 1 with a small “lead component.” For 
very large frequencies, the transfer function is considera- 
bly larger than 1, and is approximately 1 + (y7M)) 
where y is the ratio of specific heats of the gas, and M, is — 
the Mach number at entrance to the nozzle. 


Vv 


4 

ECENTLY, the problem of combustion instability 

of rocket motor has been studied by several — 
authors (1, 2, and 3).? In these investigations, it is — 
assumed that the percentage increase of the mass rate | i 
of flow through the nozzle is equal to the percentage — 7 / 
increase of pressure in the rocket cylinder. It is, how- | | | 


ever, not certain whether this assumption is correct. | 7 >X 
Since this flow condition enters in a direct manner into — = % 


the instability calculation, the relation between flow | 
variations and the pressure variations should be deter- — 
mined more carefully. It is the purpose of this paper 
todo this. The result of the present study is 


FIG. 1 ROCKET NOZZLE WITH LINEAR VELOCITY DISTRIBUTION 


equal to the product of density and velocity, the frac- 


as the transfer function of the rocket nozzle, i.e., the 
ratio of the fractional increases in mass flow sen the 
chamber pressure as function of the frequency of the 
oscillation. It indicates that the transfer function is 
rather a complex function of the nozzle geometry and 
the frequency of oscillation, and that the previous very 
simple assumption is only justified for very small 
frequencies of oscillation. 


Flow Conditions 


The flow in the nozzle will be considered as one- 
dimensional, i.e., at each nozzle section the conditions 
are taken to be uniform and the independent variables 
of the problem are then the time ¢ and the distance x 
along the nozzle axis in the direction of flow (Fig. 1). 
Let p be the pressure, p the density, and u the velocity. 
The primed quantities are the oscillating quantities; 
thus p’ is the density oscillation. Similarly, the un- 
primed quantities are the steady-state or undisturbed 
quantities. Therefore p’/p is the oscillating pressure 
at x expressed as fraction of the steady-state pressure at 
x. These fractional quantities are assumed to be 
small so that only first order terms need to be con- 
sidered. Since the mass rate of flow per unit area is 
Received January 2, 1952. 

! Robert H. Goddard Professor of Jet Propulsion. 
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tional increment of the mass rate of flow is ( ao ) 
p u 


Hence the purpose of this paper can now be stated as 


u'\ | p’ 
simply to compute the ratio - at the en- 
p u//\p 


trance to the nozzle as a function of the oscillation fre- 
quency w. This ratio will be called the transfer function 
of the rocket nozzle. 

The conditions at the entrance to the nozzle are fixed 
by the plausible assumption that the temperature of 
the combustion gas is not changed by variations in 
pressure. Let the gas be considered as a perfect gas, 
and let the subscript 1 denote the entrance to the 
nozzle; then 


Since the flow in the divergent part of the rocket 


‘nozzle is supersonic, the obvious additional condition 


for the complete determination of the flow is simply 
that the propagation of oscillations must be toward 
the downstream direction. For instance, if U(2) is the 
velocity of the propagation of density oscillations at the 
section x, then for an observer who travels with 
the speed U, the density is stationary, i.€ 


= 
rhis is the equation for determining the velocity of 
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The condition of downstream propaga- 
Thus, in the nozzle the signs 


propagation. 
tion means positive U. 


of o(? ) /ot and () az must be opposite. Similarly, 


the signs of and )joz must also be 


posite in the nozzle. we 
Formulation of the Problem in 
-_In the nozzle, if A is the cross section of the nozzle 
at x, then the continuity equation is 


+ 2 + + = 


The steady-state flow or the undisturbed flow satisfies, 
however, the simpler continuity equation 

Apu = const 
By eliminating A from Equation [3] by this equation, 
and by retaining only the first-order terms of the 
oscillating quantities, one obtains the linearized con- 
tinuity equation as 


Similarly, the dy namic equation is 


p u 


du po 
For any fluid mass, the entropy is maintained. There- 


fore, if y is the ratio of specific heats of the gas, 


Equations [4], [5], and [6] are three equations for the 
three dependent variables (p’/p), (p’/p), and (u’/u). 

It is convenient in the following calculations to intro- 
duce a specific nozzle shape such that the steady veloc- 


ity in the nozzle increase simply linearly with z 
(Fig. 1). The simplest way to specify this is 


where the subscript 1 denotes again quantities at the en- 
trance to the nozzle. 2, is thus the z-coordinate at the 
entrance to the nozzle. Therefore, the origin of the 
x-axis is not generally at the entrance to the nozzle; it 
is there only if the steady velocity at entrance is equal 
to zero. With Equation [7], Equation [6] becomes 


Therefore, if the entropy oscillations at the entrance 
to the nozzle is ee’, i.e., 


then, according to Equation [8] in general 


By eliminating (p’/p) between Equations [5] and [10], 


= € exp E (: _ log =) 


the resultant equation together with Equation [4] 
constitute a system of two equations for the two un- 
knowns (p’/p) and (u’/u). 

Now let 22 be the z-coordinate at the exit of the 
nozzle, or if 1 is the length of the nozzle, 


-_-Define a nondimensional independent variable as 


Then 
& = = +1) = + [13] 


The oscillating quantities with frequency w can now 


be written as 
[14] 


By aes these expressions into the differential 
equations for p’/p and u’/u, equations for f and g are 
derived. The results can be, however, written more 
compactly if one uses the reduced frequency 8, defined 
as 


uy Us 


where us is the exit velocity of the nozzle. £8 is then 
the ratio of the angular frequency w and the velocity 
gradient u,/z; along the nozzle axis. Then from Equa- 
tions [4], [5], and [10] the equations for the f and 9 are 


obtained as 
df 4% 
[oi + + ipf = 


and 
7 
dé 


+ 


where M is the local Mach number of the undisturbed 
flow at z or £. In fact, for the assumed linear velocity 
distribution in the nozzle, 


£,?) 


M, is the Mach number at the entrance to the nozzle. 
By eliminating g() from Equations [17] and [18], a 

single second-order equation for f(€) is obtained. 

The result can, however, be reduced to more con- 


venient form by using a new independent variable 
defined as 
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It is easy to show that z is actually the square of the 
ratio of u to the so-called critical sound speed. Thus 
z = 1 at the throat of the nozzle. In terms of z, the 


differential equation for f is _* 
: 218 df ip(2 + iB) 
1 — ip’ — 
° — i(B/2) 2y 2 ip 
is (2) + 1) + i; . [21] 


The relation between f(z) and g(z) is 


(2 + iB)g(z) = — 1) + — + 2) 


z\-6/)F — 1) + 


2. is, of course, the value of z corresponding to & as 
viven by Equation [13]; i.e., 


1 
71 = - 
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The initial condition needed to solve Equation [21] 
is specified by Equations [1] and [9]. That is, at 
= 2,0rz = 


2 


€ 
f(a) = 
The complete solution requires the cognizance of the 
fact that the propagation of oscillations must be toward 
the downstream direction as discussed in the previous 
section. 

When f(z) is determined, one can compute g(z) by 
Equation [22]. 
the velocity oscillations are determined as functions 
proportional to the amplitude of entropy oscillation e. 
Since the point of interest is the ratio of the oscillations, 
the arbitrary « does not really enter into the final 
result. In fact, according to Equation [1], the ratio 
of the fractional increase in mass rate of flow and that 
of pressure, or the transfer function G() is - 


[(2). + =1 +o [25] 


Therefore G(8) is unity as assumed by previous in- 
vestigators of combustion instability only if there 
is no velocity oscillation at the entrance, or (u’/u); = 0. 
[It is perhaps worthwhile to point out the fact that the 
present definition of transfer function G(8) corresponds 
to the practice of servomechanism analysis, if the 
pressure oscillation (p’/p), is taken as the “input” and 
the mass flow oscillation as the “‘output.’”’ Hence the 
transfer function G(8) will be useful also in the syn- 
thesis of the systems for servocontrol or servostabiliza- 
tion of the rocket motor, as suggested by W. Bollay (4). 


Solution for Small Frequency 


Then the density oscillations and 


The basic differential equation with linear velocity 
distribution in the nozzle as given by Equation [21] can 


be reduced to the hypergeometric differential equation 
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by the substitution f(z) = z-w(z). w(z) is then the 
hypergeometric function. However, the calculation is 
rather difficult. For a first discussion, only cases of 
either very small frequency or very large frequency 
will be considered here. 

If the frequency 8 is very small, the functions f(z) 
and g(z) can be expanded in terms of this parameter 
f(z; B) = f(z) + BFO(z) +... 
giz; B) = g(z) + Bg(z) +... 


These expressions are appropriate if the interval of z 
concerned does not include z = 0. Fortunately, the 
physical problem generally specifies nonvanishing wu at 
the entrance to the nozzle. Therefore, & ~ 0 and thus 
zi ~ 0. The restriction on the validity of the ex- 
pansions is thus generally complied. By substituting 
Equations [26] and [27] into the basic Equations [21] 
and [22], and then equating terms of same powers in 8, 
one has 


and 
2(1 dz? 22 dz te dz 7 + f 
Furthermore 
df™ 
= — WF 2) — + DU 2) + - [30] 
(1) 


The initial condition of Equation [24] now becomes 


= 0. [33] 
A straightforward calculation will give the following 
relation for the propagation velocity U for the density 
fluctuation defined by Equation [2], 
U_ _ _ 


L ia[fO(z) + Bf Mz) 


u df ~ dfo df ——— [34] 
z 2 E 


where u is, of course, the undisturbed flow velocity at 
the section concerned. Therefore, in order that the 
propagation speed be finite even for small 8, f(z) 
must be a constant. Hence the appropriate solution 
of Equation [28] together with the condition of Equa- 
tion [32] is 


€ 


Then Equation [30] gives immediately 
Up to this order of approximation then, G = 1, as 


shown by Equation [25]. Therefore the assumption 
of the previous investigators of the combustion in- 
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i stability of rocket motor is indeed correct if the fre- 


quency is vanishingly small. 


f(z), Equation [29], can be written as 


1 


Therefore the solution is 


log z é 1 
1 1 C 

: The quantity inside the square bracket is finite and 

positive for z > 0, while the last term is infinite at the 

throat of the nozzle where z = 1. Hence, according to 
[34], the condition for finite positive propaga- 
speed U = 0. Therefore® 


“dz 


_logz 


This equation together with Equations [31] and [33] 
then gives 


( 
—11 


44 [38] 


= 


Therefore, up to the first order in 8, the transfer function 


Therefore, when the frequency 8 is not exactly zero, 
the transfer function has a small “lead component,” 
proportional to the frequency. aS 


Solution for Large Frequencies 


If the value of 8 is very large, one must find the 
asymptotic solution of Equation [21]. The domin: ating 
terms in Equation [21] are rs, 


(6/2) 171 
4y 
For the particular integral f*, take 
— (ip/2 
= 
where Z(z) is a function of z not involving 8. There- 
fore, by retaining only the highest order term, 


dz 4 2 


’ Professor L. Crocco has very kindly pointed out that the 
condition for positive finite propagation speed U is equivalent 
to the suppression of all singularities at the sonic point. This 
observation may simplify the analysis for general values of 8. 


With f(z) given by Equation [35], the equation for 


is 
1 
log — 
=1 + fla) 4y 1-2 


( — (ip/2) 


By substituting these derivatives into Equation [40], 


it is found that 


er 
thus 


— ¢ /2) 


find the complementary function, let 


= 


f(z) = etBX(z) 
then 
4 


dz? ~ dz 


By substituting these into the homogeneous equatio1 
corresponding to Equation [40], one has 


Thus 
and 
dz y+11-: 


d(21) =0 ....{42 


The complete solution for large 8 is thus 


gue ; log (2 JJ m2) | 
t+ 


f(zjeiat = — * 


i 


The first term is the particular integral; it satisfies the 
condition of downstream propagation of disturbances. 
The same condition requires, however, that A(z) be 
negative and increasingly negative as z increases from 
z;. A,(z) does not satisfy this condition and hence has 
to be rejected. Therefore, B = 0. D is fixed by the 
initial condition of Equation [24], and is — e/y(y — 1). 


Therefore, finally 
1 
- —e | [44] 
2 — 
where Ao(z) is given Equation [42]. 


By taking the dominating terms of Equation [22] 
for large 8, one has 


f(z) = 


df 
ae) = (+ 


with f(z) as given Equation [44], then g(z) is deter- 
mined as 


g(z) + 2 
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a (Continued on page 162) 


Therefore, the transfer function G(8) for large 6 is 

GB) = 1+ 


But 2, is specified by Equation [23], therefore 


GB) =1+ - [46] 


yMy’ 
Therefore, for large frequencies the transfer function 
is again real and positive, i.e., the mass flow oscillations 
and the pressure oscillations at the entrance to the 
nozzle are again in phase. However, the mass flow 
oscillations are considerably larger than the pressure 
oscillations; the ratio is given by Equation [46]. 


Concluding Remarks 


In the previous section the calculation of the transfer 
function of the rocket nozzle is considered as a problem 
in nonsteady one-dimensional flow. For very small 
frequencies of oscillations, the transfer function is 
indeed approximately equal to unity, as assumed by 
previous investigators of combustion instability. _How- 
ever, the present calculation also shows that the transfer 
function is a complex number and has a lead component 
increasing with increase in frequency. For very large 
frequencies, the analyses give results indicating larger 
mass flow oscillations than supposed before. But the 
transfer function for infinite frequency 6 is again a 
positive real number. Therefore, the Nyquist diagram 
of the transfer function G(B) is a curve, not a line. As 
an example, for y = 1.22, and M, = 0.2, according to 
Kquations [39] and [46] 

G(B) 1 + 1.43878, 


for B< 1 
for B>1 


underlying assumptions. First of all, the nozzle shape 
is such as to give a linear velocity distribution in the 
nozzle. Then the gas is assumed to be perfect and 
nonreacting. 
and residual chemical reactions in the exhaust gas will 
modify the result. 
short nozzles, the wave length may be small in com- 


In reality, viscosity, heat conduction, 
At large frequencies or for very 


parison with the nozzle length, then the one-dimensional 
flow assumption may introduce appreciable error. For 
these reasons, it is desirable to measure the transfer 
function of the rocket nozzle experimentally. Here, 
however, one must remember that the flow phenomenon 
in connection with the combustion stability is pro- 
foundly influenced by the entropy oscillations gener- 
ated in the combustion Direct stability 
measurements on rocket motor with artificially oscil- 
lated propellant flow will be necessary. The present 
analysis indicates, however, that the nozzle could be 
cut short with practically no divergent part as the 
transfer function is not influenced by the particular 
exit velocity so long as the exit velocity is supersonic. 
The propellant flow rate will then be smaller than a full- 


processes. 
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H. S. TSIEN? 


When a thin cylindrical shell of uniform thickness is 
very rapidly heated by hot high-pressure gas flowing inside 
the shell, the temperature of material decreases steeply 
from a high temperature at the inside surface to ambient 
temperatures at the outside surface. Young’s modu- 
lus of material thus varies. The purpose of the present 
paper is to reduce the problem of stress analysis of such a 
cylinder to an equivalent problem in conventional cylindri- 
cal shell without temperature gradient in the wall. The 
equivalence concept is expressed as a series of relations be- 
tween the quantities for the hot cylinder and the quanti- 
ties for the cold cylinder. These relations give the similar- 
ity law whereby strains for the hot cylinder can be simply 
deduced from measured strains on the cold cylinder and 
thus greatly simplify the problem of experimental stress 


analysis. 


E cylinder of a solid propellant rocket is sub- 
et to very rapid heating during its short dura 
_ tion of operation. The temperature distribution across 
the thin cylindrical wall, although approximately the 
same in every section, is not linear. This condition is 
most severe at the end of the combustion of the pro- 
pellant grain. From the point of view of a materials 
engineer, this case is distinguished from others by the 
time rate of heating, which is so large as to not allow 
sufficient time for appreciable change in the structure 
of the material. The strength of the wall material 
under this operating condition is quite different from 
Bion under slow heating. This fact is clearly and con- 

— ¢clusively shown by R. L. Noland in a recent paper 
. (1). From the point of view of a stress analyst, the 
rational design of a solid propellant rocket cylinder is 
: - thus complicated by the very large thermal stress and 
variable Young’s modulus of the material across the 
wall as a result of the large temperature gradient. 
Furthermore, experimental stress determination under 
actual firing tests is rather difficult due to the short test 
time available and the high temperature. 

It is believed that for these reasons the only case 
which has been analyzed by reliable rational method is 
_ the case of rocket cylinder under uniform internal pres- 


sure. The bending stresses due to canted nozzles, the 
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Walled Cylinders 


and C.M. CHENG? 


Daniel and Florence Guggenheim Jet Propulsion Center, 
California Institute of Technology. Pasadena, Calif. 


stresses due to end enclosures, mounting lugs, etc., are 
only estimated by very rough methods. The purpose 
of this paper is to improve this situation by suggesting 
an approach which will reduce the general stress prol)- 
lem of hot cylinder to a problem of an equivalent cold 
cylinder. The equivalent problem can then be solved 
either analytically by the conventional method or d:- 
rectly by experimental stress determination. In either 
choice, the problem is believed to be greatly simplified. 
This law of equivalence between hot cylinder and col: 
cylinder can be called the similarity law. 


Stresses and Strains of a Thin- Walled Cylinder 


The fact that the thickness of the cylinder is small in 
comparison to its radius and length allows a great sim- 
plification in the strain analysis. To wit, the deforma- 
tion of every point of the cylinder can be described sut- 
ficiently accurately by the displacements of the points 
on a single surface within the wall of the cylinder. This 
surface is the so-called median surface. The position 
of the median surface is so determined that a bending of 
the median surface will not induce net extensional forces 
in the plane of the median surface, across the thickness 
of the wall. When Young’s modulus is a constant, 
as is the case for a cold cylinder, the median surface lies 
midway between the outer and the inner boundary sur- 
faces of the cylinder. When Young’s modulus is 
not a constant but decreases with increase in tempera- 
ture, the median surface is displaced toward the cold 
side, as will be seen presently. 

Let x, 6, z be the coordinate system with origin on 
the cylindrical median surface such that x points in 
the axial direction of the cylinder, @ is in the circumfer- 
ential direction, measured on the median surface, and 
z is normal to the median surface, pointing toward the 
axis of the cylinder. Let U, V, and W be displacements 
of a point (2, 6) on the median surface in the directions 
v, 6, and z, respectively. They are thus functions of .c 
and @; but not of z. Then the above-mentioned fun- 
damental simplification of thin shells can be stated as 


follows: If the direct strains in the x and 6 directions 
are e, and ég, and the shear strain y,»9, then 
oU 
W_ A ow 


ARS Journal. 


The 
the 
[1] 


amy 


N, = 


and 


M, 


whet 


May 


7 d 
> 
Pat 
‘vA 

~e 

E 

ol 
la: 

are 
stre 
the 
ben 
\ 
7 
: ~ 
j 
h 


where RF is the 


radius of the median cylindri 


‘al surface, 


or the “radius of the cylinder.”” This result is some- 
times described as the Kirchhoff bending assumption: 
Plane normal to the median surface before bending re- 
mains so after bending. 
The significant stresses in a thin shell are the direct 
stresses a, and o¢ in x and 86, and the shear stress 7,». 
All other stresses are small in comparison to these three. 
Now let 7 be the temperature of wall above a reference 
temperature, say the room temperature. 7’ is assumed 
to be only a function of thickness ordinate z, but not of 
xand @. Thus the heating of the cylinder is assumed 
ae to be uniform over the entire surface of the cylinder. 
7 This is generally very closely approximated in reality. 
% If a is the coefficient of thermal expansion, the thermal 
* expansion strain is then a7’. By Hooke’s law one has 
ed 
= 5, — voy) + al 
E 
1 
ler = vor) + al 
4 (1 + ») 
where # is Young’s modulus, and vy is Poisson’s ratio. 
T E is, of course, a function of temperature or a function 
* of z. v, however, will be assumed to be constant for 
n- | lack of definite information. By solving for the stresses, 
a- | one obtains from Equation [2], 
u- 
ts Ox Ate) + ve) — (1 + v)aT(z)} 
E 
* i + vee) (1 + v)aT(z)} [3] 
of _ 
eS 21+) 
— For thin-walled cylinders, the equilibrium equations 
M, | are expressed in terms of “sectional averages” of the 
eS] stresses given in Equation [3]. That is, one speaks of 
the normal forces NV, and No, the shearing force N the 
S| bending moments M,, Mo, and the twisting moment 
My. They are related to and 7, by the following 
equations 
On * S Ng = S Nw S 
in = — Sf o.2dz M, = —Sf ogdz = — 1 
[5] 
he he integrals in the above equations all extend across 
is the thickness of the wall. By substituting Equations 
e {1] and [3] into Equations [4] and [5], one has for ex- 
ample | 
tre = . 
and 
oU voV W 
4. 
ow vy OV ow 
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E(z)T(z)dz 
— 


The integrals again extend across the thickness of 
wall. It is evident from the above expressions 
N, and M, that considerable simplification can 
achieved by choosing the median surface in such a way 
that 


be 


[10] 


This is actually the condition to locate the median sur- 
face. Since Young’s modulus FE decreases with in- 
crease in temperature, it is seen from Equation [10] 
that the median surface is nearer to the cold boundary 
surface than to the hot boundary surface. For a rocket 
cylinder, hot inside but cold outside, the median sur- 
face is near to the outside surface. With this choice of 
the median surface, the forces and the moments are re- 
lated to the displacements by the following simpler 


equations: 
OW 
W. ow 
= De (= +R 
D: (Sa Re 26? + Mr 
aw. ow 
J = 4 
1 ew 1 ov 


Mw =(1- ( +) 

It is worthwhile to point out two facts: Firstly, the 
choice of reference temperature is quite arbitrary. | 
Changing the reference temperature will change the 
value of Nr, Equation [8]. But a corresponding ad- 
justment in the normal strain will leave the normal 
forces N,, Ne given by Equation [11] unchanged. | 
Therefore the physical problem is quite independent of | 
the choice of reference temperature. Mz is independent — 
of the reference temperature due to Equation [10]. 
Secondly, the present system of equations reduces to — 
those of the conventional thin shell theory if tempera-_ 
ture gradient is absent, or if Young's modulus is | 
independent of the temperature of material. In that — 
case the median surface is midway between the bound- 
ary surfaces, and D2, is simply the usual flexural rigidity — 
of the shell. 


R? 0100 Ox 


* 
Do = — Jf Elz)dz 
E(z)zdz 
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Nondimensional Quantities and Equations of With specified loading p, there are also eleven un- 
Equilibrium knowns, U, W, Ng, No, Neg, Mz, Me, Meg, aNd go. The 


It is useful to introduce nondimensional quantities system of equations is thus complete. 


by taking R as the reference length. Thus 


Infinite Cylinder Under Uniform Internal 


u = U/R, vy = V/R, [14] 
The simplest special case in the present general prol)- 
lem is the case of very long cylinder under uniform in- 
[15] ternal pressure. If the rocket cylinder is long in com- 
and parison with its diameter, the actual stress system dur- 
m, = M:R/D:, — mg = MgR/Dz, = Mygh/Dz, ing operation is approximated by this idealized simple 


Mr = MrR/Dz...[16] case. In this problem of infinitely long uniformly 
loaded cylinder, the forces n;, ne, and moments m,, 7 
are constants independent of ~ and 6. The shearing 
force neg and the twisting moment mj vanish. w is 


Therefore Equations [11] and [12] become now 


Ou . 
proportional to é or dE is a constant, say hy. v van- 


w+ nT 


ishes. w is a constant, and is negative in the present 


1—v/fo 
coordinate system, say Then Equations [17], 
and [18], and [20] give 
n,° =k, + vke — nr 
o*w ov + o*w +m [18] m,° = = mr 
tag ton T)..(18 
ioe « 1 - ») (So + w =) where the superscript 0 denote the quantities in the 
a 0£08 present simple stress system. When the temperature 
The equations of equilibrium in terms of forces and distribution and the material properties are specified, 
- moments are here exactly the same as the conventional Equation [23] gives the strains k, and kz in terms of the 
theory (2). The only innovation is to write them in internal pressure p, and the axial load n;. If the axial 
- nondimensional form. For this purpose, one has to load is produced by the same internal pressure, then it 
- _ define the nondimensional quantities q, qe, and p of the can easily be shown that 
dimensional sectional shearing forces Q,, Qe and normal 
pressure loading P, against z-direction as follows: 
: It is of interest to note that the bending moment m; 


Q:/Do = Qe/Do p = PR/Dp...... [19] 


and me are equal to mr and are independent of the con- 
Then the equilibrium equations of forces are ditions of loading. 


: By solving Equation [23] for k; and ke, one has 


On, =. ov 

ov ow hk, = (n° — + nr 

& Mg ov ko = - nr 

— — g l-y 1+ 


If the design condition is the maximum strain of the 
material, then Equation [25] gives the criterion directly 
from the pressure and temperature loading. 


_ ow 
Ne oe + 


=p 


The equations for the equilibrium of moments are 


Linearized Theory for General Secondary Loading 


As stated in the previous section, the actual stress 
system in a rocket chamber is approximately that of 


+ + meg - + Bq, = 0 infinitely long cylinder under uniform internal pres- 

am ‘: a (21) sure. This stress system can be called the primary 

May de + dE = stresses. The deviations from the primary stress system 

= 0) are results of bending due to canted nozzle, to end en- 

elosures, to mounting lugs, ete. These additional 
where 

; ee stresses or the secondary stresses are, however, only a 

fraction of the primary stresses. Therefore it is justified 

a ai is thus a quantity of the order of (R/b)?, b the thick- to consider the second order terms of additional stresses 

ness of the wall. There are eleven individual equations and deformations as small in comparison to the first 


order terms, and thus negligible. In other words, 
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in the of (7), (18], [20], and [21]. 
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u =ke+tw w + w’ 
a, = a! + ne = + ne’ Neg = Neg’ 
me = mr +m! me = mr + me’ = [26] 
and 
Ge = do = p=p'+p’ 


where k, and kk. are given by Equation [25]. The 
primed quantities are then the secondary deformations 
and the secondary stresses, they are considered to be 
small in comparison to the primary deformations and 
stresses. From Equations [17] and [18], one has the 
following relations between the deformations and the 


stresses, 


problem of general loading on a cold cylindrical shell, 
and can thus be treated by the known methods de- 
veloped for this conventional problem. In fact, the only 
difference bet ween the hot cylinder and the cold cylinder 
is the appearance of the term mr in Equation [28]. 
However, even this difference is trivial: The reason is 
the very large magnitude of 8B as shown by Equations 
[6], [7], and [22]. In faet, if, as is generally the case, 
Nr is of the same order of magnitude as N,", then the 
above-cited equations show that the ratio of mr/8 and 
n,° is at most of the order of b/R, where 6 is the thick- 
ness of the shell. Since the shell is considered to be 
thin, or b/R < 1, the terms involving mr in Equation 


[28] can be dropped without impairing the accuracy of 
 l—v dv’ , du the present theory. When this is done, then in the 
7 2 \o¢ 06 system of equations given by Equations [27] and [28], 
the effects of thermal stresses and variable Young’s 
: \ [27] modulus are not explicit. As far as the nondimensional 
om! = de +o +9] equations are concerned, the problem of hot cylinder is 
= mat ie identical to the problem of cold cylinder, and the basic 
Oe? de? equations are now essentially the same as that adopted 
4 =) by L. H. Donnell for his study of the stability of thin 
: 008 = 0 cylindrical shells (4). This is the basis of the similarity 
By substituting Equation [26] into the equilibrium law discussed in the following section. 5. 7 
equations, Equations [20], and dropping out the second ' 


order terms of the primed quantities, a system of lin- 
earized equations is obtained. 
further reduced by substituting the and q»’ obtained 
from the last two equations of bending moment equilib- 
rium into the third equation. The final result is the 
following system of three equations: The first is the 
equilibrium of forces in the axial direction; the second 
is the equilibrium of forces in the circumferential di- 
rection; and the third is the equilibrium of forces in 
the radial direction: 


1 — dw’ _ 0 ) 
oe? 2 2 ” dE 
1 — 0%’ 1+» Ow! 0 
06? 2 2 
and 
2 —— v t+ = 
06! fe) | 4 
mr\ Ow’ , Ow’ 
— Bp’ + (ne - de + Bp' oe? J 


These equations have been simplified on the basis that 
8 is a large quantity of the order of the square of the 
radius-thickness ratio. 

In the latter equation, p’ is the secondary load im- 
posed on the cylinder expressed as a distributed pressure 
over the surface of the cylinder, directed radially out- 
ward. If the load is a concentrated load, then it has to 
be expanded into a product of Fourier series and Fourier 
intergral as done by 8. W. Yuan (3) in his treatment of 
concentrated load on a cold cylindrical shell. Other 
types of loads can be similarly expanded. Then Equa- 
tion [28] is a system of three equations for three un- 
knowns wu’, v’, and w’. The forces and moments are 
related to these displacements by Equation [27]. 

The problem of general secondary loading as expressed 
by Equations [27] and [28] is very similar to the 
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This system can be © 


Similarity Law for General Loading 


If the problem of secondary loading is to be solved 


analytically, the results of the previous section show 
that it can be reduced to an equivalent problem of cold 
cylinder and solved accordingly. However, a more use- 
ful application of the equivalence concept lies in the 
possibility of experimentally determining the stress and 
the strain on the equivalent cold cylinder and then 
using the similarity law to determine the stress and the 
strain in the hot cylinder. There are mainly two ad- 
vantages of this semi-experimental approach: (a) 
The experiments on a cold cylinder can be done more 
easily and much more accurately than possible on a hot 
cylinder. The test period can be as long as desired and 
not limited to the short burning time of the rocket. 
(b) The stresses induced by mounting lugs, etc., are very 
difficult to approximate by simple load systems amena- 
ble to theoretical calculations. For instance, the loads 
from a mounting lug are not really a concentrated force 
To take them as a con- 
would 


and a concentrated moment. 
centrated force and a concentrated moment 
grossly overestimate the actual stress. Such difficulties 
disappear if the loading is done experimentally. 

With such experimental stress determination in mind, 
it will be convenient to have the hot cylinder and the 
equivalent cold cylinder of same general sizes. Thus 
the radius R and the length L will be the same for both 
cylinders. In order that the nondimensional differen- 
tial equations, Equation [28], be the same for the hot 
cylinder as for the cold cylinder, the parameters in these 
differential equations must be the same. That is, if a 


quantity of the cold cylinder is denoted by a bar over 
the quantity 
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The condition of Equation [29] can be satisfied by mak- 
ing the thickness 5 of the cold cylinder smaller than the 
thickness 6 of the hot cylinder. This is of course to be 
: i expected, since Young’s modulus of the hot material 
is smaller and hence material is “softer” than the cold 
material. When the thickness 5 is determined from b 
2 by using Equation [29], Dp can be computed. Then 
Equation [30] gives the internal pressure P® and the 
axial load V,° from the specified P® and N,° for the 
hot cylinder. These steps then fix the geometry of the 
-__ ¢old cylinder and the primary system of loads. 
er For the additional secondary loads, the fact that 
f Equations [27] and [28] are linear equations can be 
utilized to introduce an added freedom in specifying the 
_ A loads. Linear relations are not altered by multiplying 
= = * the variable by a constant. Therefore for additional 


loads and additional displacements, the nondimensional 
quantities for the cold cylinder and the nondimen- 
sional quantities for the hot cylinder need not be 


identical, but differ by a factor «. Thus 
(a’, b’) = du’, v', w’) 
and 


Nes Neg me's mes meg’) > en’, 
Mz", Moz’, Mg’) ey 


But the nondimensional pressure loading p is related to 

‘the actual pressure loading by Equation [19]. There- 

Ee _ fore the secondary pressure loadings P’ for the cold 

__ eylinder and the secondary pressure loading P’ for the 
hot cylinder are related through 


The ratio of the pressure load is then Doe/Do. Since 
the radius R and the length L of both cylinders are the 
- same, other types of loads such as concentrated force, 
_ or moment for the cold cylinder and for the hot cylinder 
- must also bear the same ratio. Needless to say, the 
loads for the cold cylinder must be applied at corre- 
sponding points for loads in the hot cylinder. 

The additional forces N,’, N,9’, the additional shear 
Q,’, and the additional moments N,’, M,,’ at the ends 
of the cylinder are controlled by Equations [15], [16], 
and [19]. It is easily seen that because of Equation 
[29], the ratio of these quantities for the cold cylinder 
and the hot cylinder is again Doe/Do. 

Therefore, knowing the load system on the hot cylin- 
der, one can find the corresponding load system for the 
cold cylinder. The factor ¢ for the secondary loads can 
be chosen at convenience of the experimenter. For 
instance, « might be so chosen as to make the ratio 
Dye/Do equal to unity. Then the secondary load sys- 
tem for the cold cylinder is exactly the same as the hot 
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cylinder. 
is selected and the corresponding strain on the cold 
cylinder determined by strain gages, the inverse equiva- 
lence problem is then to find the strain in the hot cylin- 
der from the test data on cold cylinder. 
Take for instance, the axial strain e,(z). For the 
cold cylinder, according to Equations [1] and [25], 
=k nr (= - z [34] 
where Z is the value of z measured from the median sur- 
face of the cold cylinder, midway between the boundary 
surfaces. Now let é, be the average of the measured 
axial strains on the outer surface and on the inner sur- 
face of the cold cylinder, and let Aé, be the difference 
of the measured axial strains on the outer surface and 
on the inner surface of the cold cylinder. Then from 
Equation [34] 


, 


re) 
By eliminating and 


from Equations [35] and 
[36], one has 


ef2) = (1 [37] 


The value of z is measured radially inward from the 
median surface of the hot cylinder and thus is larger in 
magnitude for the inside surface than for the outside 
surface. 


1 
ti 


Similarly, 
ez) = (1 - *) ke + nr + ae). 
and 


1 
D 


where @¢ is the average of the measured circumferential 
strain on the outer surface and the inner surface of the 
cold cylinder, Aés is the difference of these strains; 
¥.9 and Ay,» are the corresponding quantities for the 
shearing strain. In Equations [37] and [38], k; and 
ke, and nr are the primary strains computed from Equa- 
tions [25] and [8]. Therefore these equations allow the 
calculation of the strains in the hot cylinder from test 
results from cold cylinder, and thus complete the de- 
sired similarity law. 

For a stress analyst, the next step is perhaps the cal- 
culation of the principle strains at each value of z in the 
shell, and examine whether the larger of these principal 
strains exceeds the design limit of the material at the 
temperature prevailing at that point. 
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When the proper load for the cold cylinder 
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Example of Dimensioning the Equivalent Cold 
Cylinder 


As an example of the procedure outline in the previous 
section, the data given by Noland (1) will be used to 
find the equivalent cold cylinder. The temperature dis- 
tribution in the wall is taken from Fig. 2 of that paper 
and is reproduced as Fig. 1 here. The material is as- 
sumed to be 19-9 DL, and the variation of Young’s 
modulus with temperature is plotted in Fig. 2, again 
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FIG. 2. YOUNG’S MODULUS AS A FUNCTION OF TEMPERATURE 


using Noland’s data. First, the position of the median 
surface will be determined, using Equation [10]. This 
is found to be 0.588 b from the inside surface. Next, 
by taking the cold cylinder to be at 100° F, the ratio of 
the thicknesses of the hot cylinder and the equivalent 
cold cylinder is computed by using Equations [6], [7], 
and [29]. It is found that 
b = 0.936b 


Thus the equivalent cold cylinder of the same material 
is 93.6 per cent as thick as the hot cylinder. 

The ratio of the loads on the cold cylinder and the hot 
cylinder is controlled by Do/Do. This is computed as 


Do/Do = 


By taking a = 10-*/°F, and 6 = 0.095 in., the values 
of thermal stresses are 


Nr = 22400 lb/in., 
Mr = 192 )b. 
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1.29 


nr = 0.958 X 107? 
mr = —0.124 


no = = 0.72 x 10-3 
Now 
2 
8 = 12 (7) 
b 


If the radius of the cylinder is 2.25 in., and if the internal 
pressure P® is 1500 psi, then 
= 1.44 x 10-3 


If the axial tension N,° is due to the same internal pres- 
sure, then 


Therefore the ratio of mr/8 and n,° is 


mr ‘Bn,° = —0.0224 


This is indeed smaller than the thickness-radius ratio 
b/R = 0.0422. Hence the surmise that m7/8 is negli- 
gible against n,° is now justified by numerical calcula- 
tion. 


Junction Stress Between Cylinder and Head 


of the similarity law for 
and the cold cylinder is 


The previous formulation 
stresses in the hot cylinder 
based upon the assumption that the secondary load 
system is specified. This is not true for the junction 
stresses induced by fitting, say, a hemispherical head 
to the cylinder. Such stresses are determined by the 
equality of deformations of the head and the cylinder 
at the junction. If the semi-spherical shell has the same 
thickness as the cylindrical shell, then the temperature 
distribution in the spherical shell will be the same as in 
the cylindrical shell. An analysis shows that the 
“similar” cold test specimen can be made also of uni- 
form thickness 5, determined by Equation [29]. How- 
ever, the similarity of both the primary loading and 
secondary loading now requires an additional restric- 
tion: must now be unity. When these conditions are 
fulfilled, the similarity of junction stresses in the hot 
cylinder and the cold cylinder will be assured, and the 
relations previously developed for computing the 
stresses in the hot cylinder from test data on the cold 
cylinder remain valid. 


Ring Stiffener Around Cylinder 


To strengthen the thin cylindrical shell against con- 
centrated loads from the mounting lugs, a ring stiffener 
is often attached to the outside surface. It remains cold. 
Young’s modulus of the ring material thus is that 
of cold material. To determine the dimensions of a 
“similar” ring for the cold cylinder, the conditions to 
be satisfied are those of Equations [31]. In other words, 
for a ratio ¢ of deformation of the ring on the cold cylin- 
der to that on the hot cylinder, the force required must 
bear the ratio e«D)/Do. This means that the ratio of 
the stiffness of the ring for the cold cylinder to that on 
the hot cylinder is Do/Do. If the ring is rectangular in 
section, the correct ratio of the width of the two rings 


(Continued on page 167) 
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This paper presents the solution of the heat flow prob- 
lem in composite walls under heat transfer conditions 
which are typical of uncooled rocket engine walls. Ana- 
lytic expressions in the form of Fourier sums are obtained 
for the temperature distribution in a composite wall con- 
sisting of an inner (refractory) medium and an outer (me- 
tallic) medium under Newtonian heat transfer into the 
first medium with negligible heat transfer from the 
second medium to the exterior. The expressions obtained 
are based on a plane parallel composite slab as a represen- 

tative model for relatively thin cylindrical walls, with 
_ thickness-to-radius ratio not exceeding 0.2. The general 
results for the composite slab are simplified for the limiting 
eases of a thin refractory shield with a thick shielded me- 
_ dium and a thick refractory shield with a thin shielded 


Nomenclature 


Most of the symbols listed here are more fully defined i in the text 
as they are introduced. 
specific heat 


d = wall thickness 
h = heat transfer coefficient 
k = thermal conductivity 
t = 
= = space coordinate 
a? =k/pc = t ermal diffusivity 
> \ = characteristic constant of partial differential 
equation 
A = Fourier amplitude 
i = temperature above ambient 
T, = flame temperature above ambient 
A = characteristic constant 


The subscripts 1 and 2 are used with symbols pertaining to the 
shield and shielded media, respectively. Additional subscripts 
m,n, designate running indices in Fourier sums. 


1 Introduction 


: 4 N MANY rocket engineering applications the prob- 
i lem of transient heat transfer through composite 
walls has been gaining increasing importance. The 
typical composite wall consists of an insulating cylin- 
drical shell of refractory material within a cylindrical 
metal envelope, which, for structural reasons, must not 
attain excessive temperatures during the firing period. 
The insulator receives heat by convective heat transfer 
from the rocket flame. Some of the heat transferred 
is conducted through the interface between the insu- 
lator and envelope, and is accumulated in the envelope. 
Heat transfer from the envelope to the exterior is gen- 
erally negligible until the envelope attains excessive 
temperatures. 
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Heat Flow in Composite Slabs 


E. MAYER ' 


The M. W. Kellogg Company, Jersey City, N. J. 


The results obtained for heat flow problems in com- 
posite cylindrical shells (1, 2)? are, as a rule, too compli- 
cated to be of practical value in application to the com- 
posite rocket wall, especially because of uncertainties 
in data on material properties and heat transfer coeffi- 
cients. For typical wall thickness-to-diameter ratios it 
is sufficient to consider a composite plane parallel slab 
as a representative model of the cylindrical wall. This 
paper contains a mathematical description of the tem- 
perature transients in a composite slab corresponding 
to typical heat transfer conditions in rocket motors. 
While the general results obtained here are still rela- 
tively complicated, the solution is reduced to two 
simple limiting cases of considerable importance; 
namely, the “thin” insulator and the “thick” insu- 
lator. 


2 Heat Flow in a Simple Slab 


In order to refer the results obtained for the composite 
shell to a basic problem which has been studied in con- 
siderable detail (1, 3) we recall here briefly the tem- 
perature distributions obtained in a simple slab receiv- 
ing heat through one boundary by convective heat 
transfer and insulated at the other boundary. The 
slab is contained between the planes x = 0 and « = 
—d, as shown on the following sketch. The tempera- 
ture distribution within the ma- 
terial is governed by the partial 
differential equation 


ar _ ar 


in which the thermal conductiv- 
ity k, the specific heat c, and 
the density p are assumed 
constant. With Newtonian heat 
transfer at the flame side x = 
—d, the boundary condition at 
this plane is 


—d, =h[T, — T(—d, t)].......... [2] 


where h is the heat transfer coefficient and 7’, the flame 
temperature, both assumed constant. At the other 
boundary the adiabatic condition is assumed to hold: 


; t) 


For convenience, all temperatures 7’ are measured above 


2 Numbers in parentheses refer to the References on page 158. 
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2 
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ambient, and the initial condition satisfied by the solu- 
tions of [1] is taken as 


The method of separation of variables (4) applied to 
|1 | leads to solutions of the form 


+.C... 


in which C is a constant and )? is a separation constant. 
With the form [5] substituted in [1], the time terms can- 
cel and the spatial dependence of the temperature is 
-ontained in 


where a? = k/pe designates the thermal diffusivity of 
the material. Thus, the spatial function is sinusoidal 


f(z) = A coe + Bein ™..............16] 
a 
ws 
lhe boundary conditions [2] and [3] lead to 


[4 sin + B cos = 
a Qa 


f a 

Qa 


From the time-independent term in [2’] it follows that 
( = T,. Furthermore, since from (3’) B must be 
zero, Equation [2’] reduces, after cancellation of the 
of the factor Ae~**, to an eigenvalue equation for \: 


Ad dd hd = 


a Qa 
Introducing the dimensionless parameters 


a k 


we put the eigenvalue equation into the form 


the solutions of which are tabulated in Ref. 1 for the 
range 0 < NV < »,. Each value of N determines an 


infinite set of discrete eigenvalues A, = A,d/a, n = 
a Sere The associated functions f in [6] 

An 

f(r) = A, cos = A COS 
a d 

form a complete orthogonal set satisfying : 

0 

dr = 0 n#~m 


which permits the determination of the amplitudes A, 
by Fourier analysis of the initial condition [4] ol ; 


T(z, 0) = A,, cos X,, 
n d 


0 xr 0 xr 
A,/T, = An ‘os? A, > dr = 
0 cos A cos? A 


—4 sin A,/(2A, + sin 2A,) 


With these results the transient temperature dis- 
tribution for the simple slab receiving heat by New- 
tonian heat transfer is represented as a Fourier sum 
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sin A, 


T(z, t) = T, A. + A. “cos A, [5’ 
It may be noted here that in practice the above Fourier 
sum converges rapidly, one or two terms in the sum 
being sufficient for computation of 77/T, except for 
very short times when 7/T, << 1. Detailed curves 
of 7’ vs. t for fixed positions x in dimensionless variables 
appear in the Gurney-Lurie (3) charts, which are based 
on the Fourier solution corresponding to conditions 
when only the first term of the sum in Equation [5’] 
is of importance. 


3 Heat Flow in Composite Slabs 


The mathematical description of the transient tem- 
peratures arising in a composite wall consisting of two 
plane parallel slabs of different media under New- 


tonian heat input is obtained with reference to the co- 
ordinate system shown on the following sketch. 
One of the slabs is con- 
between the planes 

x = —d,,andx =0,the 

other between x = 0 and 

x = d,. The temperature 

T(x, t) is governed by the 

heat flow equation 


x= -d, 
where the thermal diffusivity is constant in each me- 
dium having the values - 
(as? = ke/pocs for < xr < de 
With Newtonian heat transfer at the flame side k = 
-—d,, the boundary condition at this plane is 


h{T, — T(—d, = 


where the heat transfer coefficient h and the flame tem- 
perature 7’, are assumed constant. At the outer bound- 
ary of the composite wall, « = d, the condition corre- 
sponding to vanishing heat transfer is employed: 


While actually there may occur some radiation losses 
to the exterior as the temperature at the outer bound- 
ary rises, the latter heat transfer rates, as estimated 
from radiation formulas, are negligibly small in prac- 
tice. Choosing the temperature scale so that the ini- 
tial uniform temperature of the composite wall is 
ero, we have as the initial condition 


[11] 


In brief, then the mathematical problem at hand is 
the solution of the partial differential Equation [8] sub- 
ject to the boundary conditions [9] and [10] and the 
initial condition [11]. 

The solution of Equation [8] is represented by differ- 
ent functional forms in the two regions on the above 
sketch. Thus: 
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T(z, t) = t) 


4 


At the interface x = 0, continuity conditions on the 
temperature 7 and the heat flux kO7'/Ox lead to the 
following connections between F; and F»: 


[14] 
dF,(0, t) _ , t) 


7 Further conditions on the F’s are obtained from the 
boundary and initial conditions on 7, Equations [9], 
{10}, and [11]: 


—k — Fi] at «x = —d...... [16] 
0 at [17] 

F\(z, 0) = Fx, 0) [11’] 


Seeking solutions by separation of variables we try 
solutions of the form 


F(z, t) = e—tf(x) + Ci F(x, t) = + C2 
which lead to 
d? 
with the results 
fi(z) = A; cos — AS 4. B, sin ~ 


4 
AL 
4 fox) = Az cos + By, sin 


From the boundary condition [16] we conclude that 
C, = T,, so that the time dependent terms can be 
canceled yielding 


= ky df = —hf; at 


dz 


[167] 
On comparing terms in Equation [14] we find C,; = 
C, = T, and 


The remaining boundary conditions become, upon 
cancellation of the time terms, 


3 


The four boundary conditions lead to a system of 
homogeneous linear equations for the amplitudes Aj, 


at [15’] 


at | (17’] 


B,, Ag, B: 
Ali B,| tan +5] +0 
r 


Bin _ Bin _ tam _ 
Bn _ Ban 


For nontrivial solutions to exist (i.e., solutions other 
than A; = A; = B, = B, = 0), the determinant of the 
coefficients must vanish yielding the eigenvalue equa- 
tion for A, 


ky 

tan Ay + > 

1 hd, A, tan A; 
where: ‘ 

c= dd; he = Ad» 
a ae 


For given heat transfer coefficient h and the com- 
posite wall data (material properties /, pc and thick- 
nesses d, of each medium) the eigenvalue equation 
possesses an infinite set of ’s 


as its solution. With each X, there is an associated 


function f, and the solution for the temperature dis- 
tribution is a Fourier sum of the form 


t) = +T, -d<r<90 
T(z, t) = 
n 


where the representations of f, in the two media are 


fin = Ain COS Ain 4 -d<r<0 
foal ) d, d, 
nit) = 
fon = Asn COB Asn + Bon sin Az» 


In the expressions for f, only one of the four ampli- 
tudes Aj,, Aen Bin, Bey is independent, the ratio of any 
other to this one being determined by the boundary con- 
ditions, Equations [14’] to [17’]. Thus we may put 


An = Asn =A, 


= — = tan Aon = pen ; 


and obtain 
d, d, 
= [20] 
fon = A, (cos Aen d> + Hen Sin Az, i) 


in which the unknown amplitudes A, are determined by 
Fourier analysis on the initial condition 


T(z, t)= at ¢=0.. 

It is to be noted at this point that the f,’s are not 
orthogonal over the range —d; < x < de, and hence to 
perform the indicated Fourier analysis we must first 


14") 

167) 
ae 


~ 
5 ial equ 
| 
> 
| Ar] 
: 
k df, 
dx 
| 
Ad, 
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obtain the orthogonal functions associated With the rn?) | dz = 0 
composite wall system. These functions turn out to — 
be given, as demonstrated below, by Now since A, # A», for m # n, we conclude that the 
ams errr integral must vanish and the orthogonality condition 
Jin at ont . 
e(z) = = [21] is established. 
V p2tafon OSs; We now determine the amplitude A, by Fourier an- 
where: alysis on the initial condition [11”] with the aid of the 
here: 
orthogonality properties of the functions Rewrit- 
(xr) = J ing Equation [11”] as ay 
‘To prove the orthogonality of the ¢,.’s 
we multiply both sides by the factor and in- 
=O m#n.........[2]] tegrate over the composite wall: 
we start with the differential equations satisfied by f de = f [22] 
n — 
In the summation on the left the product ¥(x)f,(a) is 
dx , written as ¢,(a) in accordance with the definition of the 
at latter. By virtue of the orthogonality condition [21], 
2 n n o/ 
dz the summation reduces to a single term f,,2(x)dv. 
We multiply the former of these equations by ¢,, the Reintroducing the explicit forms of ¢,, 
latter by ¢,, and subtract one from the other, obtaining éa(z) = 
n rearrangement and integration over the wall thick- . ; 
earrangement and integration over the wall thic V pic: Am [ cos Aim + SiN Aim —d;<z<0 
= V pt Ay, [ cos + sin Aom 0 = ds 
(+ we perform the integrations with the appropriate form 
d, ” dz? "dx? of ¢,, in each of the two regions of the composite wall. 
With the notation 
+ as dx? — dr? ) dx 
D = ‘os A z + A x 1 
Consider now the first integral on the right. After in- eA ae eee 
serting here ¢,, = V pic and ¢, = JV y - ds 2 
ig here @ fim ind pili Sim we ob f [ cos x a 
tain 0 d: 
ky [fn Tre Nee = [os Aim + Him sin Aim | az 
. . ‘ds 
which, by integrations becomes Nem = [ cos Aim + pom Sin Aem | dx 
df, n d m 
[fim fin | we have for Equation [22] 
df; Dim + p2c2Deom| im + p2C2Nom] 
Because of boundary condition [16’] ky = at 
dx 
x = —d,, and hence the bracketed terms cancel at the 4. — 7, + om (23) 
lower limit. Thus, only at x = 0 is there a contribution ais * prc:Dim + prt2Dem 
to the first integral considered. sierra the integrations for D and N we obtain 
Proceeding similarly with the second integral on the 
right we obtain cancellations at the upper limit x = d, Din = in [C1 + pim?)Aim + (1 — prim?) sin Aim COS Aim — 
by virtue of [17’] and therefore we find Quim Sin? Ajm] 
de dhe dfim Dom = {a + 2)Aom + »m?) sin Asam COS Aom + 
—_— ks [fon de” San Nim = a [sin Aim + im (Cos Aim — 1)] 
Employing the continuity conditions at the inter- 
ace 
at by which the amplitudes A,, are completely deter- 
ky ky mined. 
a a The procedure, then, of finding temperature tran- 
we obtain a complete cancellation on the right and sients 7'(z, t) in the composite wall consists essentially 
hence arrive at -— * ON of the following steps. First, with the given material 
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4 The Eigenvalue Equation for Composite Slabs 
NN While the eigenvalues A for the simple slab depend 
FIGURE Te 8 F 
: Tol Nl N ce only on the parameter N = hd/k, as given by Equation 
<7 Nd [7’], the eigenvalues for the composite slab depend on 
N AS three parameters, which can be exhibited by rearrange- 
, xe ment of the eigenvalue Equation [18] in the form 
x 
where 
ANS kopot2 d, ds A, 
CESS \\ 
2 Ne a A simple graphical method of determining the eigen- 
NU values of A; from [18’] involves the following steps. 
The parameters o, and are computed for the 
Ny given heat transfer coefficient and material data. 
SAS \ Next, each side of Equation [18’] is plotted as a func- 
mai 1 aN tion of A;. The branches of the tangent on the left of 
~~ [18’] intersect with the branches of the cotangent on 
Ls CHS AAINAK the right at an infinite set of points having the eigen- 
_ 
values as the abscissae of intersections: A; = Aj, 
Li be S Ayw...Ai,... From these the corresponding and 
As. => rA;,, r,, = = 
3 10 9 —e 30 109 By the above graphical procedure the first two eigen- 
values Ay, were determined for a broad range of the 
constants, dimensions, and heat transfer coefficient parameters Nj, o, and » as shown on the curves, Figs 
we determine the eigenvalues \, or the related eigen- ‘The all 
values Ai, = And from the eigenvalue Equation [18]. was based on the consideration that, in general, the 
- = . shielding material has a lower thermal conductivity 
hese immediately yield the time dependence and the than the envelope material, while the heat capacity 
amplitude ratios r unit volume, pc, for most materials is of the same 
~ pe pe, 
tan - order of magnitude. Thus pic) ~ and < 1 
limit the range of interest of o in rocket applications to 
won = tan Ao, = tan Ay ~ 
ay} 


Finally, the amplitudes A, are computed from [23] 
and the spatial functions f, are constructed as given by 
[20]. In this manner all the terms entering into the 
temperature distribution T(x, t) given by [19] or its 
alternate 


t A n ‘os A n + Min sir A 
0 1 B 1 In In 
—d, < I < 0 


le A a(cos Aon z Hen sin Aon 7) 
O<r<d, 


T(z, t) = 


are uniquely determined in terms of specified material 
and heat transfer coefficient data. 

In the application of Equation [19’], the first two 
terms in the Fourier sums are usually sufficient for an 
accurate computation of the temperature in the com- 
posite slab. The first two eigenvalues necessary for 
computations with these Fourier terms, are given in the 
next section for a wide range of heat transfer and ma- 
terial parameters, 
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5 The ‘Thick’? Thermal Shield 
A considerable simplification in the calculation of 
temperature transients in composite walls can be ef- 
fected when, either because of low shield material con- 
ductivity, k; << hs, or because of large shield thickness 
d,; >> ds, the envelope material has no appreciable tem- 
perature gradients. Indeed it may be readily visual- 
ized that whenever the heat transferred through the 
interface at x = 0 is sufficiently small, the envelope is 
instantaneously at the interface temperature, 7(0, ¢). 
Under these conditions the envelope material behaves 
as though its thermal conductivity were infinite, and 
the temperature distribution within the shield can be 
deduced from the general results of Section 3 by a limit- 
ing process in which kz ©. 
Consider, for example: 


| Aen lim Ande 0 
ke 
lim pen = lim tan Ao, = 0 
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leading to the result that the spatial variation of the 


is wil The eigenvalue Equation [18’] can be simplified by 
noting that for large kz 


1 
cot nA — 

Wa’ 


ki pres de 
ke dy 


a2 


d; ds pci 


so that 
& 
lim o cot 7A; = lim 29,2 [24] 
ko—>o Ai 7 Ay p2C2d2 


and, t 


therefore, Equation [18’] reduces to 


al The first two eigenvalues computed from ™ latter 


1 picid, 


Ay piled 


are plotted on Fig. 2 for the range 
5 k, 
30 
10 
Ae | 
! 
3 
— | 
FIGURE ID 
THE FIRST TWO ROOTS OF THE EIGENVALUE 
EQUATION FOR THICK SHIELDS: tena, + ton". test 
hd, 
30 ] 
3 
5 1.0 Cp de 
Ac 


For the determination of the expressions D,, N,, which 
enter into the amplitude formula [23], the limiting proc- 
ess yields the values 


lim Da» = d, 


kr 


lim Non = d, 

ko 

while D,, and N,, remain as derived previously in 
Section 3, with the additional restriction implicit in 


[24] 


A pices = 


Hence, the amplitude formula [23] is reducible to _ 


temperature within the envelope material vanishes: lim -A, = — 27, ; 
lim fon lim cos Aon = Hen sin Aon = An sin Ain + Hin COS Ain 
ko—> de ver + pin®)Ain + (1 — pin?) Sim Ain COS Ain + 2pin COS? At, 


. [237] 


With these results Equation [19’] represents the tem- 
perature distribution in the composite slab in which the 
shielded medium behaves essentially as on thermal ca- 
pacity of magnitude pected. per unit area, being at any 
instant at the interface temperature 


T(0, t) = Ty + Ane—>ntt 


It may be noted here that, while the limiting case 
considered in this section can be formulated directly as 
a simple slab problem for the shield material receiving 


heat from the flame as in the boundary condition 
o7T(—d,, t 
and transferring heat to the thermal capacity according 
to the boundary condition 
27,0) 

Or 


[3”] 


o7(0, t) 
ot 


the method of separation of variables does not immedi- 
ately lead to the solution because of the time derivative 
present in [3’]. Thus, while the separation of varia- 
bles yields the eigenvalue Equation [18’], the spatial 
eigenfunctions obtained are not orthogonal over the 
range —d, < x < 0 occupied by theslab. [Owing to 
the presence of the time derivative, this boundary con- 
dition is not of the Sturm-Liouville type which assures 
orthogonality (4)]. In order to perform the Fourier ex- 
pansion of the initial condition, one must first construct 
orthogonal functions from the ones yielded by [2’] 
and [3’]. The limiting procedure employed in this 
section obviates the need for the construction of orthog- 
onal eigenfunctions from those obtainable with the 
boundary conditions [2”] and [3”]. a 


6 The “Thin” Thermal Shield 


Owing to the role of the expression d?/a? in many 
heat flow problems, this ratio represents roughly the 
heat diffusion time, 7, across a thickness d of a material 
of diffusivity a?. From this point of view the parame- 
ter 7 is, by definition. related to the diffusion times 
7, and 7: in the two components of the composite slab 


by 
picid;? ke Ai 


For increasing values of n the eigenvalues, A; diminish, 
as can be seen from the curves on Fig. 1, and for the 
first two eigenvalues Equation [18’] may be written in 
the approximation 


tan (4 + #) (1 + = ¢ cot nAy... 


This, in turn, may be put into a form for Az = mAj, 


. [18] 


a analogous to the eigenvalue equation for a simple slab 
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Furthermore, for 7 >> 1, the first two amplitudes A, 
and A, in Equation [23] are given approximately by 
4. Ne 7 Azn COS 93" 
so that in this case the spatial function reduces to 


Axl cos Ao, on SIN Aon i.) = 


4 sin Ao, Cos Aan( 1 


2QAon + sin 2Aon 


lim fz, = lim 
n>>I1 n>> 


which is analogous to the result for the spatial depend- 
ence of the temperature in the simple slab Equation 
[5’]. Here, of course, the ‘‘flameside’’ corresponds to 
the coordinate x = 0, and the insulated boundary is at 
do, 

Thus, for sufficiently large n, the problem of the tem- 
perature distribution in a shielded envelope can be re- 
duced to solving for the temperature distribution in a 
simple slab, with the material data of the envelope and 
a reduced heat transfer coefficient given by [25] pro- 
vided that the first two eigenvalues A,, = Ao,/n com- 
puted from the solution of [18’’’] satisfy the approxi- 
mation in [18”]. 

The expression for an effective heat transfer coeffi- 
cient [25] for thin shielding layers may be inferred 
from a discussion (Ref. 1, p. 15) of thin surface skins 
of oxide, grease, scale, and other materials of poor con- 
ductivity, in which the thermal capacity may be neg- 
lected, and, hence, the temperature gradient may be 
assumed uniform at any instant. With the notation of 
Section 3 
\Fi(z, t) 
(F(x, t) 


-d<2r<0 


we obtain for the instantaneously uniform-temper- 
ature gradient in the shield material 


oF; 


or [T, — Fi(—d,, t)] = const... [16’] 


which, upon integration, yields 
F\(z, t) iv. = F\(—d,, t)] + F,(0, 


Now, the interface boundary condition becomes 


oF 0, t) oF. —d, t) 
oF.(0, 2) 

k .. [157] 


Substituting the value of F,(—d,, t) obtained from 
{16”]: 


t) = 
1 


and employing the second interface condition [14] to 
replace F,(0, t) by F2(0, t), we obtain from [15’] 
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As tan A: = ™ i+ N, [18 ] 

with an effective (reduced) heat transfer coefficient = 

het = [25] 


1+ N, 

which shows that the temperature distribution F»2 in 
the shielded material is subject to a boundary condition 
of Newtonian heat transfer with an effective heat trans- 
fer coefficient as given by [25]. With the boundary 
conditions [15”] and [17], the treatment of envelope as 
a simple slab leads to the eigenvalue equation [18’’’], 
but the validity of neglecting the thermal capacity of 
the shield material need be justified by satisfying the 
approximation in [18”] for the first few eigenvalue 
which enter into the practical computations of the 
temperature distribution. 


7 Application to Composite Cylindrical Walls 


In application to rocket engineering, the tempera- 
ture transients arising in relatively thin cylindrical 
walls are, in practice, computed with sufficient accu- 
racy by treatment of the equivalent plane, parallel slab 
problem. It can be shown by detailed study of the tem- 
perature transients in the simple cylindrical shell under 
Newtonian heat input that, if the thickness-to-radius 
ratio does not exceed 0.2, the equivalent slab treatment 
yields results accurate within 0.03 7, on the conservative 
side from the design point of view, since the cylinder 
temperature does not exceed anywhere, at any instant, 
the corresponding slab temperature. The results of 
the analysis presented here for composite slabs are also 
conservative in application to cylindrical walls, but ap- 
pear sufficiently accurate in practice if the total thick- 
ness-to-radius ratio does not exceed 0.2 
As an example of the application of the thin shield 
formulas we consider the effect of zirconia on a low- 
-arbon steel nozzle in prolonging the duration time be- 
fore allowed temperatures in the steel are exceeded. We 
may limit the allowed temperature, 7',, on the flame 
side to 7,/2, where the flame temperature, 7',, is of 
the order of 4500 F. As a representative value of the 
heat transfer coefficient in the nozzle we set 


Btu/hr 
ft? deg F 


and determine the duration times in a '/;-in. thick, low- 
-arbon steel wall for the following cases: (a) no shield- 
ing; (b) shielding by an 0.02-in. thick layer of zirconia. 

The material properties employed are listed in the 
table below, those for zirconia being taken from Ref. 
(5). 


Property SSS Material— 
Low carbon Insulating 
Thermal conductivity 
ft? (deg F/ft) 
Density p lb/ft® 485 175 
Specific heat 
Btu/Ib deg F — 
Diffusivity 
, ft? 
ec hr 


7 
> 
, 
> 
, 
> 
7 
» 
) 
; 
; 
; 


_ 1500 
k 20 


25 
1.56 


which yields from Equation [7 


A, = 100 As = 3.56 
With these the first two exponents in Equation [5’] be- 
come 
= 6301" = = 79508 


where ¢ is in units of hours. Inserting the above re- 
sults in Equation [5’] corresponding to the flameside 

0, we obtain the equation for the duration time 
tp in units of seconds: 


1— .625¢ —-175tp — .095e—2-2ltp = — = 
with the result tp = 1.25 sec. 

In case (b) the thin shield formulas [18’’’] and [25] 
involve the reduction factor 


: 1500 
1+N, =1+4+- 3 (3) = 9.3 


yielding the effective Biot number and heat transfer 
coefficient 
N h 


= 1.68 hess = 1 => 161 


Btu/hr 


Var = 
ft? deg F 


Use of the simple slab formula [5’] with this V,,, shows 
that only the first term in the Fourier sum is of im- 
portance, and that the first eigenvalue is approximately 
given by 


A, tan A, = A,? = Neg 


whereby [5’] reduces to 


ie 

indicating that the shielded metal now behaves essen- 

tially as a thermal capacity with negligible gradients 

within it. The prolonged duration time tp’ in units 


of seconds is _ 
‘yy 
= 3600 log 24 sec 


The thick shield formulas developed in Section 5 are 
convenient in the design of insulating refractory liners 
for rocket chambers with heating times up to the order 
of a minute, as in some liquid propellant units. Such 
units are usually regeneratively cooled, but the possi- 
bility of using composite chamber walls which are com- 
petitive from the weight point of view may be con- 
sidered. In a composite wall the weight per unit area 

w = pid, + 
may be minimized subject to the following considera- 
tions. The metallic envelope must be sufficiently 


thick to carry the pressure loads during firing, i.e., 


’| the first two eigenvalues 


maintain the inequality in [27]. 


- which minimizes the weight per unit area w. 


where p is the chamber pressure, R the chamber radius, 
and Sz is the yield stress in the envelope material. The 
yield stress S. is, however, a function of the temper- 
ature 7, generally diminishing as 72 increases and, 
therefore, as the duration of the firing tp increases. 
Now, for any desired duration time tp, it is possible to 
select an infinite set of thicknesses (d,, d2) large enough 
In the set (d;, 
there is, in general, a particular pair of thicknesses 
As a de- 
sign criterion we may, therefore, consider the optimum 
choice of (d;, d2) such that, for specified materials and 
heating parameters, the pair (d;, d2) shall make w in 
Equation [26] a minimum and lead to equality in [27] at 


time = tp. 


The comparison of weights for two zirconia-metal 


- composite walls with reference conditions in a chamber 


given by 
h= 300 53 deg F T, = 4500 deg F pk = 1400 in2 ie 
is shown on Fig. 3. The desired duration time has been 
FIGURE IL COMPOSITE WEIGHT PER UNIT AREA w VS, THICKNESS dy 
A 
= = 
= 
ms i 
x 


06 8 
MAETAL THICKNESS, 


taken as tp = 60sec. For a regeneratively cooled unit 
operating under the above chamber conditions, the 
weight per unit area of the stainless-steel jacket wall 
and the included coolant is w = .040 Ib/in.*. From 
Fig. 3 it is seen that the zirconia-molybdenum com- 
bination appears to be weight-wise competitive with 
the regeneratively cooled unit up to 60 seconds of firing 


duration. 
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4 
= 4 o 
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Compas 


An analysis has been made of the expansion of gases 
through a de Laval nozzle under isothermal conditions. 
It is shown that for ideal gases, the increase in kinetic 
energy of the gases is equal to the heat transferred into the 
gas to maintain constant temperature. The initial energy, 
or temperature, enters only in the design of the nozzle and 
determination of the pressure ratio required for the trans- 
ference of the fixed quantity of heat. Unlike the adiabatic 
case, the jet velocity is not limited by the initial energy 
content of the gas but only by the pressure ratio or energy 
content of the reservoir transferring the heat. Formulas 
are derived for isothermal exhaust velocity, mass flow, 
and thrust coefficient which are analogous to those for 
adiabatic flow. 


Nomenclature 
a = isothermal sonic velocity = WR,T 
a’ = adiabatic sonic velocity = 
c* = mass flow constant (chamber pressure per unit throat area 


per unit mass flow rate) 
Cp = specific heat at constant pressure 


cr = specific heat at constant volume 
Cr = nozzle thrust coefficient 

dq = increment of heat energy | 

e = base of natural logs _ 
E = internal energy 

f = nozzle section area 

F = thrust exerted due to flow through nozzle 
F’ = free energy 

H = enthalpy 

m = dm/dt = time rate of mass flow — 
p = free-stream pressure 

R; = specific gas constant = R/M 

S = entropy 

7 = free-stream temperature 
v= mean gas velocity 
= specific volume 
Y = Cp/cy = ratio of specific heats 

= prescribed function of 

= nozzle divergence correction factor 
p = density of gas 

Subscripts 

= adiabatic 

« = chamber or initial 

- = exit or ambient 

i = jsothermal 

o = atmospheric or standard 

= throat 


Received July 9, 1951. 

1 Staff Engineer. Member ARS. 

2 Chief, Chemistry Section. Member ARS. 
3 Numbers in parentheses refer to References on page 162. 


May-JuNE 1952 


of Adiabatic and Isothermal Expansion 
— Processes in Rocket Nozzles 


H. S. SEIFERT' and D. ALTMAN? 


Jet Propulsion Laboratory, California Institute of Technology, 


by convective heat transfer, or most probably by ab- 


4 


A 
‘ 


nozzle 


Pasadena, Calif. 


HE adiabatic expansion of gases through a de 


Laval nozzle is well understood and fully docu- 
mented (1 and 2).* It is of some theoretical interest to 
compare the behavior of such a nozzle with that of the 
less familiar one in which the expansion process takes 
place isothermally. A partial treatment of such a 
nozzle has already been given in this JouRNAL (3). An 
isothermal expansion would be difficult to realize in 
practice, since it implies that energy is supplied continu- 
ously during expansion to the axial translational degree 
of freedom of the gas molecules. It might be achieved 
to a limited extent by the release of chemical energy, 


sorption of radiant energy. An interesting consequence 
of isothermal expansion is that the velocity of the gas 


increases without limit as the pressure ratio increases. 


™ 
Conservation of Energy 


Assume that an isothermal process is possible. 
will be shown that the thermodynamic function le 
free energy F’, which has been found useful in chemical 

‘alculations (4), will play the same part in isothermal 
‘alculations as the we H does in adiabatic 


calculations. The free energy F’ is defined as 


 F’ = total he: ” — unavailable or reversible heat 
= K+ TS 


F’ represents the energy which is available (i.e., free) 


to do external work. ( ‘onsider the equations of energy 


-and momentum. 
ei 


Energy: 


dq = dE + pdV =d(E + pV) 
Momentum: 
From Equation [2] 
nil dp = vdv = . d(v?) 
p 2 


Therefore, Equation [1] becomes 


dq = + pV) +d G 


= 
| 
{ 
is 
d 
| 
| 
+pV +5 — dq =0............[3] 
159 
7 


By definition 


For the isothermal case 7’ = const, and 
T dS = d(TS8)... 


Using Equation [4] in Equation [3] 


equation (cf. Equation 24 of Ref. 5) for v? is 


If numerical values are inserted in Equations [10] 
and [11], with v2 neglected, the following ratios for iso- 
thermal to adiabatic exit velocity v;/v_ result (take y = 
1.25 in the adiabatic case, and R,T, , the same in each 


vy? — = 


l : 
ing a(E + pV +5 *) — d(TS) = 0 
ease): 
(E + pV — 78) + dot = const 20:1 
2 10021 
FY + gv? = const...............-. [5] This tabulation indicates a 16 to 36 per cent so-called 


This equation may be compared with the usual adia- 
batic energy relation pn 


Exhaust Velocity 


The relations [5] and [6] are general and apply to all 
real fluids. To calculate the exhaust velocity for an 
isothermal nozzle, it will be necessary to reduce general- 
ity by considering a fluid; for example, an ideal gas 
whose equation of state can be analytically expressed. 
It is possible to write an equation of state (6) for rocket 
combustion products, which is a good approximation to 
real gas behavior because the temperatures are high, 
pressures are moderate, and corrections due to inter- 
molecular forces are low. This discussion will be lim- 
ited, however, to the ideal nonreacting gas case. 

Assume an equation of state 


Combining Equations [2] and [7] to eliminate p, 


- = —R,T In p + const 


vdv = 


[8] 


Thus it may be seen from Equation [5] that the free 
energy for isothermal expansion of a perfect gas is 


= R,T |n p + const...... [9] 


and that this isothermal process is one in which F’ 

converted into ordered kinetic energy of the gas. If 
subscript c indicates chamber conditions, the integra- 
tion constant in Equation [8] may be evaluated by 


noting that when p = p,, then v = v, and T = 
Thus 
It may be seen from Equation [10] that as p > 0, 


v—> o, unlike the adiabatic case in which v approaches 
a fixed limit v,,., a8 p—> 0. The conventional adiabatic 


‘ The authors are indebted to H. S. Tsien for this derivation of yy 


Equation [5]. 


improvement in exhaust velocity as the pressure ratio 
goes through the indicated range, and the limiting per- 
centage of improvement is infinite. The limiting infi- 
nite value results from the fact that the nozzle absorbs 
an infinite amount of energy from the heat reservoir 
in this case. 

Equation [10] may be regarded as a special case of 
Equation [11] in which y approaches unity. Equation 
[11] becomes indeterminate of the form 0/0 as y — 1.0. 
By applying L’Hospital’s rule and taking the limit, 
Equation [11] transforms into Equation [10]. 

Since cp = yR,/(y — 1), the physical meaning of + 
= | is infinite specific heat. The more nearly y > 
1.0, the smaller the temperature change associated with 
a unit amount of external mechanical work performed 
adiabatically by the gas. The more complex the gas 
molecule, the greater its ability to store energy at a 
given temperature; that is, the greater the number of 
ways it can absorb energy in other than translational 
motion and the nearer y approaches unity. A gas 
with y = 1.0 has in effect an infinite reservoir of heat 
within each molecule and can expand both adiabati- 
cally and isothermally at the same time. A real gas 
with y > 1.0 may expand isothermally, but only if heat 
is supplied to it from external sources. 


ood | 
Mass Flow and Throat Area 


Equation [10] may be expressed in terms of area f 
rather than pressure p by using the equation of state 
(Equation [7]) and of continuity 


to eliminate p,/p; thus 


and Equation [10] becomes 


v 


Another combination of the state and continuity equa- 
tions (Equations [7] and [12]) yields 
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From the form of Equation [15] it is evident that a 
throat will exist (minimum of f) since f > © as v ap- 
proaches either 0 or o. To find the throat values of v 
and f, denoted by a subscript ¢, differentiate Equation 
[15] with respect to v and equate to zero. This opera- 
tion gives 


=elv /2R. = 0... [16] 


or 


Equation [17] shows that the gas velocity v, is also 
the isothermal sonic velocity at the throat. This equa- 
tion is to be compared with the adiabatic case (cf. 
Equation 37 of Ref. 5) for which 


The throat area f, in the isothermal case may now 
be calculated from Equation [15], neglecting v,? rela- 
to v/, 


Equation [19] may also be regarded as a-means for 
‘alculating the mass flow rate m. If it is desired to de- 
fine a characteristic velocity c* = p,f,/m analogous to 
the adiabatic case (cf. Equation 52 of Ref. 5), then 

[20] 
where a, is the sonic velocity for a gas with y = 1. 
This c* may be of value in calculating mass flow and 
as in the adiabatic case may be considered as a measure 
of propellant performance. The corresponding adia- 
hatie definition of c* is 


(y+D/2(1-—y) 


ax —=§ : 
+1 
where a’, is the adiabatic sonic velocity in the combus- 
tion chamber. As in the case of Equation [10], Equa- 
tion [20] may be derived from Equation [21] by taking 
the limit as y > 1.0. 


(y¥+1)/2(1—y) 
ae 


Thrust Coefficient and Expansion Ratio 


In order to calculate thrust F and thrust coefficient 
Cr, customary definitions are used which are unaf- 
fected by the nature of the thermodynamic process in- 
volved, namely, 


and 


F = Sipe (in (pe/Pe) 2+ (pe Do)fe [24] 
2 (p/p iis + 125) 
Pc f: 


Equations [24] and [25] are to be compared with the 
adiabatic relations (cf. Equations 49 and 50 of Ref. 5) 
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F=f [1 - (®) ] + (pe — po)fe 
1 De 


...- (26) 


where 
2 (y+ /2(7-D 
= yi 
i) 

In order to calculate the pressure ratio p,/p, used in 
Equation [25], a relation between area ratio and pres- 
sure ratio is needed. If one takes Equation [10] and 
changes variable from v to f by the successive use of 
2quations [12] and [7], then neglecting v,? 


m\2 R,T.m ) 2 De 
= ( —. = (22,7 .) n= 


Eliminating m with Equation [19] gives 

fe V 2e [In (pe/p))!/2 

which may be compared with the corresponding adia- 
batice equation (cf. Equation 45 of Ref. 5) 


. [29] 


To find the critical pressure ratio p,/p;, one may dif- 
ferentiate f in Equation [28] with respect to p and 
equate the result to zero, or simply set f = f, and solve 
In either case 

Pe/pt = Ve = 1.65... 
which is independent of y (as are all of the other iso- 
thermal relations). Equation [30] is to be compared 
with the usual adiabatic relation (cf. Equation 35 of 
Ref. 5) 


for 


= 1.80 for y = 1.25 
As before, Equation [30] can be obtained from Equa- 
tion [31] by taking logarithms and applying L’Hospi- 

wy 
Conclusions 


tal’s rule as y — 1.0. 


It may be well to reiterate certain conclusions con- 
cerning the isothermal process for ideal gases which 
might otherwise be obscured among the algebra: 

1 Unlike the adiabatic case, none of the performance 
parameters of an isothermal nozzle depends upon spe- 
cific heat (or y). 

2 There is no limit to the theoretical exhaust veloc- 
ity that can be obtained with a finite initial pressure, if 
sufficient heat is available in the reservoir. 

3 The kinetic energy of the exhaust gases is ob- 
tained from the heat Aq which is transferred to the ex- 
panding gases. This heat is shown to be equal to the 


change AF’ in free energy. 
4 The limitations of such a process in practice 
would probably arise from inadequate heat transfer 


mechanisms. 
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4 The Transfer Functions of Rocket Nozzles = 


(Continued from page 143) 


length nozzle, a desirable saving in experimentation. 
The present investigation also gives the characteristic 
nondimensional frequency as the angular frequency 
w divided by the velocity gradient along the nozzle. 
The transfer functions for rocket nozzles of different 
sizes can be correlated through this nondimensional 
frequency or reduced frequency. 

In conclusion then, it is certain that the present 
paper gives only a general outline of the problem of 
transfer function of rocket nozzles. The complete 
detail has yet to be filled in by further studies. 
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Jet Propulsion News 


= C. F. WARNER, Purdue University, Associate Editor 
with the assistance of W. G. BOHL 


Rockets 


HE Russians are working on the further develop- 

ment of the following German rocket-propelled 
missiles: the Wasserfall anti-aircraft rocket; the A9, a 
winged version of the V-2; the A5, a half-scale V-2; 
and a rocket of the V-2 type. The work is being done 
at the reactivated Peenemunde missile center. New 
facilities, including six new launching ramps for a Mata- 
dor-type missile, have been installed by the Russians. 
The airfield at Peenemunde West has been greatly ex- 
panded and over 80 MiG-15 interceptors are stationed 
there for the defense of the area. The base is the most 
tightly guarded “closed area’”’ in the entire Soviet zone 
of Germany. A recent giant explosion within the 
guarded area prompted volunteer firemen of a nearby 
village to rush to the zone. They were immediately 
pushed away by Russian police. 


+ + + 
A PICTURE of the Armstrong-Siddeley “‘Snarler”’ 
liquid propellant rocket motor is presented in Fig. 1. 
This unit, reported in the January-February issue (p. 39) 
of the JouRNAL, consists of an externally driven propel- 


% 


Fic. 1. 


ARMSTRONG-SIDDELEY ROCKET POWER PLANT FOR AIR- 
CRAFT 
lant pumping unit, propellant valve assembly, and thrust 
chamber. The system operates as follows: When the 
pilot switches on his master switch, the fuel and liquid 
oxygen pumps start torotate. The liquid oxygen pump 
primes, then the fuel pump primes, and the igniter 
valves open automatically. The two liquid propellants 
flow to the igniter chamber, and the spark plug fires. 
This causes the igniter to light and to provide a pencil of 
flame in the combustion chamber. The main bipropel- 
lant valves now open automatically and the propellants 
flow into the combustion chamber where they are ig- 
nited by the pencil of flame. Finally, the rocket motor 
fires, the current to the spark plug is cut off, and the flow 
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of propellants to the igniter chamber ceases. The mo- 
tor may be operated at part-thrust, also. 
+ + + 
DURING a recent test conducted by U.S. Army Ord- 
nance officers, the guided missile rocket, Nike, inter- 
cepted and “destroyed’”’ a remote-controlled crewless 
B-17. The Nike has been ordered into mass produc- 
tion. 
+ + + 
RECENT news stories have aroused considerable 
interest in the guided missiles being developed in this 
country. Only the names and general field of applica- 
tion of these missiles are available at this time. 
DEVELOPMENT 
NAME AGENCY 
Navy Bu Ord. 


MANUFACTURER SERVICE 


Anti-aircraft de- 
fense 
Sparrow Navy Bu Aer — Sperry Air-to-air missile 
Nike Army Ord. Bell Telephone, Anti-aircraft de- 
Douglas fense 
G. L. Martin Co. Medium-range 
bombardment 


Terrier Convair 


Matador Air Force 


Air Force Northrop Air- Long-range bom- 


craft bardment 


Snark 


Hermes Army Ord. General Electric Short-range bom- 


bardment 

Loki Army Ord. Bendix Anti-aircraft de- 
fense 

Raseal Air Force Bell Aircraft arma- 
ment 

Lark Navy Bu Ord. — Fairchild Guided Anti-aircraft de- 
Missile fense 

+ + + 


ANOTHER assault on the world’s altitude record for 
single-stage rockets will be made next spring by a Navy 
Martin Viking. The forthcoming Viking will be some- 
what different in appearance from its predecessors, 
however. The redesigned Viking No. 8 is greater in di- 
ameter, slightly shorter, and has fins three-quarters the 
size of the first seven Vikings. This new size will al- 
low considerably more fuel to be carried, thus providing 
a longer thrust period. The power plant of the present 
Viking, made by Reaction Motors, Inc., develops ap- 
proximately 20,000-lb thrust. 

+ + + 

GEORG M. Lorenian of West Germany has made two 
patent applications to the Swiss Federal Patent Office 
at Berne covering the use of liquid propellant rockets to 
brake the descent of men and equipment dropped from 
aircraft and to provide a certain amount of steering 
control. An unmanned model has been built which 
weighs 120 lb. Nitric acid and a mixture of crude xyli- 


dine and triethylamine (“Tonka 250’) is foreed by 
compressed air from the supply tanks into the two com- 
bustion chambers. 


Failure of an air valve reduced the 
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first model to scrap metal during the first drop. A sec- 
ond model is now undergoing bench tests. 
+ + + 

ARMY Rocket Research fights TB by voluntarily 
deferring delivery of hydrazine, a chemical compound 
showing promise as an improved rocket propellant. 
Hydrazine is also an essential part of a new family of 
antituberculosis drugs which show encouraging results 
in treatment of tubercular patients. Only large-scale 
manufacture of hydrazine has been sponsored by Army 
Ordnance Corps; release of the compound into com- 
mercial channels will assure continuation of clinical 
tests. Although many TB cases which appeared to be 
hopeless have shown improvement after treatment, the 
drugs are not yet considered a cure. 

+ + + 

THE Navy has taken off the secret list a World War 
II gadget called the ‘Gesundheit mine,”’ developed 
under Project Sneeze. Its details may be of interest to 
rocket engineers with under-water applications in 
mind. 

Both the Allies and the enemy had proximity mines 
during the war which would explode without actually 
being struck by a ship. The trouble was that the prox- 
imity mine was no respecter of nationality. It was as 
effective against friendly ships as against those of the 
enemy. 

Project Sneeze was an attempt to ‘‘make the mines 
choosy.”’ The heart of the modified mine was a small 
projectile launcher which, aimed by magnetic attraction 
at a passing ship, would throw aboard a canister loaded 
with red pepper. 

“A sensitive microphone pickup on the mine did the 
rest,’’ the Navy reported. “If a crew member on the 
ship said ‘Gesundheit’ the mine would explode; if he 
said ‘God Bless you’ the mine would remain unarmed.” 


New Use for Rockets 


A NEW type of burner, similar in operation to the 


rocket motor, has been developed by Linde Air Products 
Division of Union Carbide and Carbon Corp. The 
Linde Jet Piercing Machine, as it is known, is used to 
cut holes in rock. The flame produced by burning 
kerosene in oxygen leaves the burner with a velocity of 
6000 fps and literally disintegrates the rock upon which 
it is directed. One of the machines is now being used 
in Minnesota in the mining of the iron-producing taco- 
nite ore. Because of its extreme hardness, taconite is 
expensive to mine by conventional drilling methods; 
however, it is hoped that the Linde-machine will greatly 
reduce the mining cost. Mechanical drilling machines 
now drill from 1 to 2 ft per hr in taconite, whereas a 
jet-piercer can do 18 ft per hr, and has done as much as 
35 ft per hr. 


Pulsejet Engines 


THE French have used six Escopette pulsejets to 
power an Emouchet glider for over 200 test flights. 


The Escopette, because of its low weight, low speed 
range, and good fuel consumption, has previously been 


developed in low power sizes. However, SNECMA is 
now developing a number of higher speed and higher 
powered models. Each of the Escopette units develops 
approximately 22-lb maximum thrust and weighs 11 Ib. 

The Escopette, a valveless pulsejet, illustrated dia- 
grammatically in Fig. 2, operates upon a pulsating cycle. 


RECUPERATOR 


INLET DIFFUGER TAIL PIPE 


— 


COMBUSTION 


FIG. 2. ESCOPETTE PULSEJET DESIGN 

When in operation, air entering the inlet diffuser is 
mixed with fuel as it enters the combustion chamber. 
Autoignition occurs and the pressure rises to approxi- 
mately two atmospheres. A major portion of the hot 
products of combustion are then discharged through the 
tail pipe giving the unit a forward thrust. A small frac- 
tion of the products is also discharged through the dit- 
fuser inlet, thus causing an undesirable thrust tending 
to retard the forward motion of the unit. Because of 
this, the Escopette has been fitted with a recuperator 
which directs these products rearward, resulting in an 
increased over-all thrust. The elimination of the me- 
chanical valves gives the unit long life and simplifies the 
injection of fuel. The actual design of the inlet diffuser 
is the result of many tests by SNECMA and represents 
the most original feature of the Escopette. The pas- 
sage is shaped so as to retard the flow of gas outward. 

The unit is started by introducing short periodic in- 
jections of fuel and compressed air into the mouth of 
the inlet diffuser. Ignition is obtained by means of an 
electric spark plug. As soon as ignition takes place the 
spark is switched off. It has been found that the pulse- 
jet may be restarted with ease in flight at a speed of ap- 
proximately 47 mph. 

Fuel pressures between 0.3 and 1.4 atm are used, and 
it has been found that the unit will not operate at pres- 
sures outside of this range. The thrust can be varied 
from 6 to 22 lb with fuel consumption from 2 to 5 gram / 
sec by changing the type of fuel injector. When three 
units were mounted close together, the thrust of each 
was decreased by 10 per cent at the same fuel consump- 
tion. 


Radio to the Moon 


SCIENTISTS at Sterling, Va., recently trained a huge 
radio receiver on the moon and received the message in 
Morse code: “What hath God wrought.’’ The mes- 
sage originated at the Collins Radio Company’s trans- 
mitter at Cedar Rapids, lowa, was reflected from the 
moon’s surface to the Sterling receiver. This message, 
the first intelligible one to be transmitted by lunar re- 
flection, required only 2'/. sec to complete the 480,000- 


ARS JourNAL 


ai 
| 
> 
| 
‘ 
t 
a 
G 
W 
Ww 
| el 
sv 
ra 
i 
45 
f 


mile round trip. The possibilities opened up by the 
success of this venture are: Accurate radio control of 
satellites, and unjammable radio guidance—when the 
moon is right—of long-range guided missiles. 


= Airecraf 


FIG. 3 is the first photograph released by the Navy of 
the Chance Vought F7U-3 Cutlass twin-jet fighter. No 
performance data are available at this time. 

+ + 

THE Army Ordnance Department has revealed that 
its one-man helicopter weighs about 300 Ib, can be 
folded into a 5 X 5 X 14-ft container, and has a top 
speed of 80 mph. The helicopter carries a 90-min fuel 
supply of gasoline, kerosene, or diesel oil, and is powered 
by two 16-lb pulsejet engines mounted at the tips of the 
rotor blades. 


FIG. 3. CHANCE-VOUGHT CUTLASS” JET FIGHTER 
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+ + + 
A GRUMMAN F9F-6 Cougar has undergone recent 
flight tests by the Navy. This plane, rated in the 
“over 600 mph” class, is the swept-wing successor to the 
Navy’s Panther. 
+ + + 
A FULL-SIZED delta wing section of molded plastic 
will be displayed in London shortly. The new molding 
technique was developed at the British government’s 
Royal Aircraft Establishment. Britain plans to build 
a complete plastic supersonic aircraft within two years. 
THE revolutionary design of the Hawker-Siddeley 
Gloster GA5 aircraft is shown in Fig. 4. The delta 
wing jet is powered by twin Sapphire engines, each of 
which has a thrust of over 7000 Ib, or more than the four 
engines of a heavy Superfortress bomber. 
+ + + 
McDONNELL Aircraft Corporation has received an 
Air Force contract to build its Voodoo XF-88A twin-jet 
swept wing fighter. This plane, having the longest 
range of any U.S. jet fighter, will be known as the F-101. 
+ + + 
NORTH American Aviation Inc. has developed a new 
45-deg swept wing fighter, designated the F-100, to be 
powered by the Pratt & Whitney J-57 jet engine. 
+ + + 
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THE Air Force plans to procure a large number of new 
Douglas twin-jet-engine bombers, designated the B-66. 
This plane, a modified version of the Navy’s swept 
wing A3D carrier-based bomber, will weigh over 77,000 
Ib and will be powered by the 9500-lb thrust General 
Electric J-73 turbojet with an afterburner. 

+ + + 
_ EXPERIMENTAL flight tests are being conducted 
on the SO06026 by the French Company SNC A du Sud- 
Ouest. The S06026, powered by a Hispano ‘‘Nene,”’ is 
equipped with an auxiliary rocket engine in the rear of 
the fuselage. 

+ + + 

THE Russians are reported to have a new twin-jet 
night fighter having a wing span of about 26 ft and a 
length of about 28 ft. It is believed that axial flow M- 
012 turbojets having a rating of 6600-lb thrust, power 
this fighter. 

+ + 

UNITED States machine tools combined with French 
engineering under the Mutual Security Program are 
making possible the quantity production of a French 
jet fighter, the Mystére, said to outperform the Ameri- 
can Sabre jet. The Mystére MD452, manufactured 
by the Bordeaux plant of the Marcel Dassault Aircraft 
Company, was test-flown in a “duel’’ with the Sabre 
jets, and the American pilot of the Mystére was greatly 
impressed by its greater maneuverability. The French 
plan to be making thirty of these planes per month by 
the end of the year and to use them in the French Air 
Force. 

+ + + 

THE Swedish jet fighter, the Saab J-29, is reported to 
have attained a flight speed of about 657 mph. The J- 
29 is supposed to have a rate of climb superior to the 
Sabre and to have a longer range than the MiG-15. 

+ + + 

WHEN completed, a fourth prototype of the French 

Leduc ramjet aircraft will be used as basis for a super- 


FIG. 4. HAWKER-SIDDELEY LONG-RANGE JET FIGHTER 
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sonic aircraft development. The third machine is pow- 


ered by two Turbomeca turbojets mounted on the 


tips. 


+ + + 
A NEW turbojet unit of 13,230-lb thrust is to be de- 
signed and constructed by Turbomeca for new French 
fighters. 
+ + 
BOEING Airplane Company has been awarded an Air 
Force contract “for an engineering study of the applica- 
tion of nuclear power plants to aircraft,’’ the Depart- 
ment of Defense announced today. A special group of 
scientists and engineers is engaged in this new project. 
The contract, under which Boeing is working in close 
co-operation with Pratt & Whitney Aircraft of East 
Hartford, Conn., was awarded last year. No details 
may be revealed as to the scope of the program, or as to 
specifications of the power plants or airplanes under 
study. 
THE first Series-2, long-range de Havilland Comet 
fitted with four Rolls-Royce Avon engines (6500-lb 
thrust) has made a successful 2-hour maiden flight. 
This Comet is capable of carrying 44 passengers, mail, 
and freight for over 7500 miles nonstop. 


Meetings of Interest 


THE Applied Mechanics Division of The American 


Society of Mechanical Engineers will hold a Shock and 
Vibration Instrument Symposium at Pennsylvania 
State College, June 19-21, 1952. The program will 
include papers in the following categories: 1. Survey 
and historical; 2. One-degree-of-freedom instruments; 
3. Sources of error and secondary effects; 4. Peak 
reading gages; 5. Calibration; 6. Applications. It 
is expected that the transactions of the symposium, will 
be published in a single volume. 
7 + + + 
~ A TRAINING course in Statistical Quality Control 
will be held at the University of Colorado, Boulder, 
Colorado, June 17-27. 

The Fourth Combustion Symposium will be held at 
the Massachusetts Institute of Technology, September 


THERMAL indicators made of paper, capable of 
determining temperatures from 115 F to 500 F instantly, 
have been developed by the Army Quartermaster 
Corps. They consist of white pigment coatings on 
black paper, each coating designed to melt at a definite 
temperature. When heated to the designated tempera- 
ture, the coating melts and disappears into the black 
porous paper revealing a black background. 

+ + + 

BELL Aircraft Corp. has constructed a vibration test 
tower 40 ft high and 15 ft square. Aircraft and missile 
components weighing up to 3000 Ib can be vibrated at 
from | to 15 eps with a maximum amplitude of 5 ft. 

+ + + 

THE General Electric Co., at its Thomson Labora- 
tory, Lynn, Mass., has developed a new method of fabri- 
cating ges-turbine stator and rotor blades. The method 
is said to require less production equipment, save up to 
39 per cent of critical material, and to cut production 
costs by 55 per cent. Two turbojets equipped with 
these blades have undergone 500-hr endurance tests and 
the blades showed no undesirable conditions or signs of 
failure. Details of the fabrication process are available 
through the USAF to other manufacturers in the nation. 

+ + + 

PUNCH-CARD calculators are now being used at the 
Purdue University Rocket Laboratory to make all per- 
formance test calculations. Raw calibration and test 
data are marked by hand upon the punch cards. The 
punch-card machines make all of the desired calcula- 
tions and print the results upon a set of performance re- 
sult sheets. 

+ + + 

ENGINEERS of the General Electric Company are 
using a new “curved light’’ quartz probing instrument 
to study flames in jet engines. The heart of the new in- 
strument, known as a flame detector, is a metal-encased 
pliable quartz filament which transmits light through 
its length. The visible light from the flame is trans- 
mitted through the quartz to a phototube whose output 
is amplified and recorded. The General Electric engi- 
neers use this detector to study the flames in jet engine 
combustors and tail pipes, and to observe flame leakage 
into the turbine blades. 


New Instruments New Facilities 


SCIENTISTS at ONR’s Naval Research Laboratory 
have developed a “‘Vortex’”’ thermometer which allows a 
correct measurement of ambient air temperature from 
planes flying 500 mph. In the NRL thermometer, a 
tube extends axially into the airstream. Within the 
tube a spiral vane produces a vortex motion. The cen- 
ter of the vortex is a region of lower pressure, which has 
a cooling effect exactly cancelling the heating effect of 
the compression caused by the plane’s motion. It is 
said that the new thermometer will record the correct 
air temperatu re regardless of plane speed. 
+ + + 


ON March 19, Lieut. Gen. Edwin W. Rawlings spoke 
at the dedication of the new General Electric Company’s 
gas turbine plant at Lockland, Ohio. In his dedication 
speech, Gen. Rawlings stressed the surprising advance 
in the design and production of aircraft made by Russia, 
and the great need for expanded effort on the part of 
U.S. industry. The Lockland plant, now employing 
7500 in the production of the G. E. J-47, will eventually 
employ 12,000 workers. 

+ + + 
AIR flow studies at Mach number 10 (ten times the 


speed of sound) are being made in the hypersonic wind 
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tunnel at the Naval Ordnance Lab., White Oak, Md. 
The tunnel has a 5-inch working section and is continu- 
ous in operation. Air is supplied at 3000 psi, and 
heated to 900 F in order to prevent liquidation under 
operating conditions. Both pressure and_ optical 
measurements of air flow are made, in study of air flow 
characteristics at high speeds and basic aerodynamics 
for later use in missile and projectile design. 
Personalities 
7% 
SDWARD H. HEINEMANN, Chief Engineer of Douglas 
Aircraft Company’s El Segundo, California Division, 
has received the Sylvanus Albert Reed Award of the 
Institute of the Aeronautical Sciences. Mr. Heinemann 
received this award for his aeronautical investigations 
xccomplished by the Skyrocket in 1951. 
+ +> + 
Ricuarp W. Powe. has been promoted to the posi- 
tion of Chief Technical Supervisor of the Physical 
Measurements Division of Aerojet. Mr. Powell has 
been with Aerojet since 1950 and was formerly asso- 
‘iated with G. M. Giannini & Company. 
+ + + 
Masor-GENERAL THOMAS F. FARRELL, assistant gen- 
eral manager of the U. 8S. Atomic Energy Commission 
since August 1951, has resigned to become Director of 
ARO, Ine., which manages and operates the Arnold 
Engineering Development Center at Tullahoma, Tenn. 
+> + + 
MaJsor-GENERAL KENNETH D. NicHois, Deputy Di- 
rector of Guided Missiles within the office ot the Secre- 
tary of Defense, has been named Army service member 
of the Department of Defense Research and Develop- 
ment Board, in addition to his other duties. 
+ + + 
JouN P. LONGWELL, member ARS, a project leader in 
the Process Division of the Standard Oil Development 
Company, has been appointed Chairman of the NACA 
Sub-Committee on Combustion. ALBERT J. BLACK- 
Woop, assistant director in the Research Division of the 
Development Company, has been appointed to the 
NACA Sub-Committee on Aircraft Fuels. 


Fluctuations in a Spray Formed by Two Impinging Jets 
(Continued from poge 131) 


sheet velocity and frequency approached a direct 
proportionality. 


= 
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A Similarity Law for Stressing Rapidly Heated Thin- Walled 
Cylinders 


(Continued from page 149) 


for complete similarity is then Do/Do. The ring on the 
cold cylinder is thus wider than on the hot cylinder. 
When this condition on the ring dimension is satisfied, 
the simple similarity relations for the stresses in the 
shell are again correct. 
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American Rocket Society News 


H. K. WILGUS, Associate Editor 


During the Second International Astro- 
nautical Congress held in London, Septem- 
ber 3-8, 1951 (see the November 1951 issue 
of the JouRNAL) the groundwork was laid 
for the formation of an International 
Astronautical Federation. A constitution 
has now been drafted and the National 
Board, after reviewing it carefully, has 
agreed unanimously that the American 
Rocket Society should become the Ameri- 
can voting member of this organization. 
The other voting members are the British 
Interplanetary Society and the German 
Gesellschaft fiir Weltraumforschung. 
Some of the General Aims, Objects, and 
Activities of the IAF are: (1) To promote 
and stimulate the achievement of space 
flight, and to attempt to insure that space 
flight shall be achieved as a peaceful, not 
a military project; (2) to do all in its 
power to secure the widespread dissemina- 
tion of technical and other information on 
space flight through the medium of ex- 
change of publications, collaboration on re- 
search, etc., as between its members; (3) 
to stimulate work on astronautical sub- 


ARS Becomes Voting Member of [AF 


jects by international and national re- 
search establishments, universities, com- 
mercial firms, individual specialists, etc.; 
(4) to hold, at one-year intervals if possi- 
ble, or less frequently if so desired by its 
Council, international meetings to be at- 
tended by delegates from among its mem- 
bers; (5) to found and administer an inter- 
national astronautical research institute, 
charged with the task of achieving space 
flight for nonmilitary purposes. 

The present IAF officers elected at the 
London Congress are—President, Eugen 
Sanger; Vice-Presidents, G. Loeser (Euro- 
pean affairs), and Andrew G. Haley 
(American affairs); Secretary, J. Stemmer. 

The permanent secretariat of the IAF is 
to be located in Switzerland. It is believed 
that this decision to affiliate the American 
Rocket Society with the IAF is a progres- 
sive and constructive action, the outcome 
of which must await the test of time. The 
terms of the IAF constitution under which 
the American Rocket Society has agreed 
to become a member permit withdrawal at 
any time. 


Third International 
Astronautical Congress 


The German Gesellschaft fiir Weltraum- 
forschung (Stuttgart-Zuffenhausen) has 
announced that the Third International 
Astronautical Congress will be held from 
September 1-6, 1952, at Stuttgart, Ger- 
many. The G. f. W. will be host to the 
member societies this year. Last year 
the British Interplanetary Society was 
host at the London Congress. 

ARS delegates to the Congress will be 
selected and announced at a later date. 
Members of the ARS who may be able 
to attend the Congress are requested to 
forward their names to Mr. Andrew G. 
Haley, Vice-President (American affairs) 
IAF, 1101 Connecticut Avenue, N. W., 
Washington, D. C. The program for the 
Congress will be published in the next 
issue of the JouRNAL. 

Papers Invited 

The ARS is now considering papers on 
astronautics, space travel, and related 
subjects for presentation before the assem- 
bled national societies at the Third Inter- 
national Astronautical Congress. 

The subjects of these papers may be 
technical or philosophical in nature, but 
should represent original thought-or con- 
cept. In the absence of the author, papers 
will be read by an ARS delegate. 
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To be considered for presentation, 
papers must be received in the New York 
office by August 1, 1952. Address: Secre- 
tary, American Rocket Society, 29 West 
39th Street, New York 18, N. Y. 


2000 Hear von Braun 
Lecture at Washington- 
Baltimore Section Meeting 


The Washington-Baltimore Section 
meeting on March 19, 1952, held in the 
Administration Building Auditorium, 
Naval Ordnance Laboratory, White Oak, 
Md., was a lecture open to the general 
public. The speaker was Wernher von 
Braun, rocket expert, and technical direc- 
tor for guided missile development, Red- 
stone Arsenal, Huntsville, Ala. If the 
response to his subject, ‘‘Let’s Tackle the 
Space Ship,” is any indication, Washing- 
tonians are impatient to get to Venus or 
Mars or waypoints in between. 

More than 2000 persons heard the lec- 
ture, according to H. J. Archer, president 
of the Washington-Baltimore Section. 
The auditorium was filled, loudspeakers 
had been installed in the corridors, and a 
large number of cars were turned away by 
the time Dr. von Braun had commenced 
his talk at 8 o’clock. 

Dr. von Braun described the possibility, 
within the next 10 or 15 vears, of a man- 
made satellite—a 250-ft wide, wheel- 


shaped space station—that could be either 
a force for peace or a weapon of war, de- 
pending on who makes or controls it. This 
artificial moon, which will be carried into 
space, piece by piece, by rocket ships, will 
have specially designed powerful telescopes 


ROCKET SHIP AND SPACE STATION MODELS ARE VIEWED BY THE 


Photograph Courtesy of Naval Ordnance Laboratory 
LECTURER, WERNHER 


VON BRAUN, AND GUESTS AT THE WASHINGTON-BALTIMORE SECTION LECTURE 
(Lejt to right): LCmdr. F. C. Durant, III, ARS vice-president; Dr. Wernher von Braun, 


rocket expert; Willy Ley, author; 


Harry J. Archer, Jr., president of ARS Washington- 


Baltimore Section; and Andrew G. Haley, ARS general counsel 
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attached to optical screens, radarscopes, 
and cameras. It will travel along a celes- 
tial route 1075 miles above the earth, com- 
pleting a trip around the globe every two 
hours, at a speed of 15,840 mph. Nature 
will provide the motive power, a balance 
between its speed and the earth’s gravita- 
tional pull. 

Dr. von Braun also envisioned the huge 
rocket ship capable of carrying a crew and 
some 30 or 40 tons of cargo into the ‘‘two- 
hour” orbit. He described this as a three- 
stage rocket, 24 stories tall, weighing 7000 
tons, and its three power plants driven by 
« combination of hydrazine and nitric acid. 
The first stage will be powered by 51 
rocket motors having a combined thrust of 

1,000 tons, and consuming a total of 5250 
ions of propellants in 84 seconds. The 
-econd stage (middle section) will have 34 
rocket motors with a 1750-ton total thrust, 
«nd burning 770 tons of propellants. The 
third stage (nose section), carrving the 
crew, equipment, and pay load, will have 
5 rocket motors with a 220-ton combined 
thrust. This “body” or cabin stage car- 
ries 90 tons of propellants, including ample 
eserves for the return trip to earth. 

Following the lecture, a panel of experts 
liscussed various aspects of interplanetary 
travel. The moderator was Milton W. 
Rosen, section head, rocket section, 
Rocket Sonde Research Branch, Naval 
Research Laboratory. The  speaker’s 
panel consisted of Dr. Martin A. Mason, 
formerly with the Beach Erosion Board 
and now Dean of the School of Engineer- 
ing, George Washington University, Wash- 
ington, D. C.; Dr. Paul B. Pearson, chief, 
biology branch, Atomic Energy Commis- 
sion; Rear-Admiral Kenneth H. Noble 
(retired), now of Aerojet Engineering 
Corporation; and Esterly C. Page, 
formerly vice-president in charge of engi- 
neering, Mutual Broadcasting System, 
and now head of the E. C. Page Consulting 
Radio Engineers, Washington, D. C. 


X-1 Experimental Rocket 
Plane Described at 
Indiana Section Meeting 


On February 12, 1952, in joint session 
with the Purdue University chapter of the 
Institute of the Aeronautical Sciences, the 
ARS Indiana Section held a meeting at- 
tended by approximately 300 persons from 
the central Indiana area, including ARS 
members from Indianapolis and Kokomo. 

R. Bruce Foster, head of the rocket re- 
search division, Bell Aircraft Corporation, 
discussed the problems in designing, con- 
structing, and testing the six X-1 aircraft 
built to date. He began his talk with a 
film depicting the X-1 in operation to 
familiarize the audience with the airplane. 
The X-1, Mr. Foster pointed out, was con- 
ventional in its design, being a straight- 
winged aircraft built as nearly as possible 
symmetrical about its longitudinal axis. 
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BY WASHINGTON-BALTIMORE SECTION ON 


MARCH 19, 1952 


(Left to right): Dr. Paul B. Pearson, Atomic Energy Commission; 
dean, school of engineering, George Washington University, Washington, D. C.; 


Dr. Martin A. Mason, 
Esterly C. 


Page, ARS member; and Rear-Adm. (ret.) Kenneth H. Noble, of Aerojet Engineering Corp. 


He emphasized the extreme strength of 
the aircraft and the fact that the design 
loading of the structure was 18 times the 
force of gravity. For a comparison, he 
pointed out that the thin wings of the X-1 
could support the entire weight of a B-36 
without structural failure. He concluded 
his talk with another movie of the X-1 in 
action, climaxed by its taking-off and 
climbing to altitude under its own power in 
contrast to the usual procedure of air- 
launching it from a B-29 airplane. 

Mr. Robertson Youngquist, chief de- 
velopment engineer for Reaction Motors, 
Inc., then spoke on the general progress 
and development taking place in the 
rocket power plant field in the United 
States. Of the several related rocket en- 
gine design and installation problems, the 
most difficult, he stated, were the corrosive 
action of propellants, the pumping pres- 
sures involved, the space limitations, low- 
temperature valving, and the high com- 
bustion temperatures and pressures. A 
film showing the growth of the Reaction 
Motors, Inc., from meager beginnings to 
its present facilities at Dover, N. J., was 
narrated by Mr. Youngquist. The film 
included pictures of test-firings and com- 
parisons of rocket engine performance up 
to and including the V-2 power plant of 
57,000-Ib thrust. 


Aerojet Rocket Engineer Speaks at 
Purdue 

Mr. B. L. Dorman, presiJtent of the ARS 
Southern California Section, and chief of 
the rocket test division, Aerojet Engineer- 
ing Corporation, was the guest speaker at 
the Indiana Section meeting on March 7, 
1952. Mr. Dorman’s talk, illustrated with 
slides, centered on present-day problems 
in the rocket propulsion industry and 
some of the goals which are a constant 
guide to improvement. 

The two major problems at present con- 


fronting the industry, Mr. Dorman said, 
are the need for better propellant combina- 
tions or entirely new propellants, and bet- 
ter materials from which to construct 
nozzles, high-speed pumps, and other com- 
ponents. An important phase of the 
search for new materials embraces exper'- 
mental work in the field of ceramics for 
use as nozzles and combustion chamber 
liners. Some ceramic components, he 
pointed out, are now in limited production. 

The four goals which will require con- 
siderable developmental work rather than 
any radical technology changes, Mr. Dor- 
man said, are increased reliability, im- 
proved performance, light weight, and the 
ability to mass-produce needed units. 

Mr. Dorman concluded his talk with 
some comments on the employment pros- 
pects in the rocket industry. Although 
the industry is still small, it is advancing 
rapidly, he said, and can offer rewarding 
and secure employment to all those inter- 
ested in this field of engineering. 


Indiana Section Tours Purdue Rocket 
Laboratory 


On the afternoon of March 15, 1952, 
about 35 members and guests of the Indi- 
ana Section were conducted on a specially 
arranged tour of the Purdue Rocket 
Laboratory. Clair M. Beighley, research 
associate, Purdue University, in charge of 
the tour, explained various pieces of equip- 
ment in the laboratory and described 
briefly the work being done at Purdue Uni- 
versity for the Navy’s SQUID Project. 

Moving pictures of test firings were 
shown to the group following the tour. 
One film, in color, depicted the different 
phases of construction of the laboratory 
buildings and also the initial test firings. 
The climax of the film consisted of shots, 
taken from two different angles, of a spec- 
tacular explosion of one of the rocket en- 
gines in its test cell. 
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Dr. Gamow Addresses 
Washington-Baltimore 


Section Meeting 
wg On April 25, 1952, the ARS Washington- 
Baltimore Section held a meeting in the 
ir Department of Commerce auditorium, 
Washington, D. C. Approximately 500 
members and guests were present to hear 
_ Dr. George Gamow, professor of physics at 
_ George Washington University, Washing- 
ton, D. C., give a lecture on ‘Rockets in 
Relativity Theory.”’ 

The meeting was opened by Harry J. 
Archer, Jr., president of the Section. 
Andrew G. Haley, ARS General Counsel 
and chairman of the Section’s program 
committee, introduced Dr. Gamow who is 
_ known throughout the scientific world for 

his theory of radioactive decay and for his 
: is recent work on the application of 

“oe nuclear physics to the problems of stellar 

evolution. 
_ The main part of Dr. Gamow’s talk was 


a discussion of a fanciful relativistic 
- rocket—one traveling near the speed of 
light. He showed that with this type of 

7 rocket it would be possible to visit the 

_ Andromeda nebula (680,000 light years 

distant) within three or four hours, due to 
‘ the slowing of time with velocity. 
Following the lecture there was a panel 


—_— under the direction of C. J. 
os Libby, project engineer on the Matador 
missile for the Glenn L. Martin Company. 
of the panel included George P. 
_ Adair, formerly engineer of the 
_ Federal Communications Commission and 

- now a consulting engineer in the fields of 
radio, communications, and electronics; 
Dr. Ralph J. Havens, head of the Atmos- 
pheric Physics Section in the Rocket Sonde 
Research Branch, Naval Research Labo- 
ratory; Dr. G. project scientist, 
Air Force Cambridge Research Center, 
and vice-president, International Astro- 
nautical Federation; and Lt. Gen. Ken- 
neth B. Wolfe, USAF (ret’d), president of 
the Oerlikon Tool and Arms Corporation. 
The Washington-Baltimore Section is 
planning a dinner meeting on May 27 at 
the National Press Club, Washington, D.C. 


chief 


Loeser, 


- Northeastern New York 
Section 


The newly organized Northeastern New 
York section of the ARS reports a recent 
election at which the following officers were 
chosen: President, F. M. Cooper; Vice- 
President, Kurt Berman; Secretary, S. W. 
Crawford; Treasurer, A. J. Orsino;. Di- 
rectors, L. R. Bollinger, H. Lopez, and H. 
M. Weber. 

On April 8, 1952, Kurt R. Stehling, 
of the Bell Aircraft Corporation, gave a 
lecture on “Satellite and Space Flight,” 
accompanied by sound and color films of 
the Bell Aircraft’s X-1 rocket. 

The formation of the technical lecture 
series on jet propulsion elements has proved 


a successful undertaking and has at- 
tracted many new members to the ARS 
group. Several Genera! Electric Company 
engineers and scientists have presented 
talks on various jet propulsion topics, 
The series will be concluded with lectures 
by Edwin H. Hull, ARS Fellow member, 
Thomas M. Pettey, Jr., ARS member, and 
Dr. R. W. Porter, ARS Board of Directors. 


Alabama Section Meets 


The second meeting of the recently 
organized ARSsection at Redstone Arsenal, 
Huntsville, Ala., took place on Feb. 20, 
1952. Dr. Harold W. Ritchey acted as 
president pro tem, and the following offi- 
cers were elected to serve until the regular 
annual election in November: President, 
Dr. Martin B. Schilling; Vice-President, 
Hans Hueter; Treasurer, Arthur Andrews; 
Secretary, James L. Stamy. 

After the business meeting an address 
was given by Krafft A. Ehricke, chief of 
the gas dynamics section of the Propulsion 
and Fuels Branch, Guided Missile De- 
velopment Group, Redstone Arsenal. The 
subject of his talk was “The Development 
of Guided Missiles and the Foundation of 
Interplanetary Flight.”” Mr. Ehricke ex- 
plained the application of rockets to space 
travel, pointing out the problems to be 
encountered and their possible solutions, 
and how present developments might be 
applied to space travel. 

On March 26, 1952, the Alabama section 
(formerly Redstone), with its newly chosen 
officers presiding, held a meeting to elect a 
Board of Directors. The members of the 
Board are: Dr. Kurt Debus, Krafft A. 
Ehricke, James J. Fagan, Gerhard Heller, 
George Henderson, Dr. W. M. Mebane, 
Dr. George Messerly, Dr. Harold W. 
Ritchey, and Dr. Wernher von Braun. 

At this meeting a set of By-Laws was 
adopted, and with its acceptance by the 
ARS National Board, the group has now 
become an officially recognized Section 
of the ARS, known as the “Alabama 
Section of the American Rocket Society, 
Inc.” 

The program for the evening was an 


illustrated lecture, “Civilian Help Needed 
for Space Rockets,” given by Pfc. Geo. 8. 
James, assistant technical editor of the 
Guided Missile Development Group, and a 
director of the Rocket Research Institute, 
Redstone Arsenal. 

Mr. James traced the history of the 
rocket engine from the early efforts of Dr 
Robert H. Goddard, through World Wa: 
II, to the present, where a sufficient rocket 
engine technology exists, he said, to con- 
struct a space station. 

Mr. James pointed out that the need fo: 
supplementing present military and tech- 
nological rocket and guided missile de- 
velopment with civilian education could be 
accomplished through the medium of 
rocket societies. As an example of how 
such societies could participate, colo: 
films showing the activities of the Rocket 
Research Institute were presented. 


ARS New York Section 
Visits IBM 


On Friday evening, March 21, 1952, 
some 100 New York Section members, 
through the courtesy of the International 
Business Machines Corporation, met at 
the IBM offices, 590 Madison Avenue, 
New York, N. Y., for an inspection of the 
IBM Selective Sequence Electronic Cal- 
culator, one of the earliest large calculators, 
and smaller electronic computers and other 
machines. After a general explanation of 
the various calculators and a résumé of 
their potentialities, small groups were 
guided around to observe how these 
mechanical calculators advance scientific 
knowledge by relieving mathematicians of 
the necessity for doing masses of computa- 
tions. 

According to brochures distributed to 
the group, the Selective Sequence Control 
Calculator reads numbers involved in the 
problem and reads the instructions for its 
solutions. By means of a “central nervous 
system,”’ the program devised by the scien- 
tist for the problem in hand directs auto- 
matically the sequence of operation, 
selects the proper number from the various 


RECENTLY ELECTED OFFICERS OF NEW ALABAMA SECTION 
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Arthur Andrews, treasurer; Dr. Martin B. Schilling. president; Hans Hueter, 
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memory units or from reference tables, 
directs them to the calculating unit, guides 
the calculating processes, and routes the 
results back to the proper places in the 
memory unit. This busy “brain’’ will tell 
the machine to record the results of its 
work when they are obtained. 

At the conclusion of the evening the 
ARS members were shown a computer 
working on actual rocket problems. As 
the last members were departing, the ma- 
chine was busy computing the answer to: 
“(f we shoot off Rocket X at White Sands, 
vill it land outside the reservation?” 


Copies of BIS Journal 
Still Available to ARS 


Members 


The American Rocket Society currently 
ieceives from the British Interplanetary 
society, on an exchange basis, a limited 
number of copies of its Journal. Since 
insufficient copies are received to permit 
distribution to all ARS members, the 
oard of Directors has decided to dis- 
‘tribute these Journals upon request, within 
‘he limit of copies available. A handling 
‘ee of $1.00 is charged for a subscription to 
all six 1952 issues. Address your request, 
together with check or money order, to 
Secretary, American Rocket Society, 29 
West 39th Street, New York 18, N. Y. 


ARS National Board Sets 
Up Space Flight 
Committee 


The extent to which the American 
Rocket Society should devote its energies 
and techniques to the problems of space 
travel has long been a question facing the 
national organization. Several proposals 
in connection with space flight recently 
have been placed before the National Board 
of Directors, and to meet this growing in- 
terest, C. W. Chillson, ARS president, has 
appointed an Ad Hoe Space Flight Com- 
mittee. This committee, of which Andrew 
G. Haley, ARS General Counsel, is chair- 
man, consists of R. W. Porter, ARS 
Board of Directors, F. C. Durant, III, 
ARS vice-president, Kurt R. Stehling, Bell 
Aircraft Milton Rosen, 
Naval Research Laboratory, John R. 
Youngquist, Glenn L. Martin Company, 
and W. E. Zisch, Aerojet Engineering 
Corporation. 

The Committee is “to investigate and 
recommend to the Board the ways in 
which the American Rocket Society can 
best serve its members and fulfill its aims 
outlined in the By-Laws, by activity in 
the field of space flight, satellite stations, 
astronautics, etc., and to define the extent 
to which the American Rocket Society 
should become active in these fields.” 


Corporation, 


New ARS Publicity 
Committee Organized 


7 At the April 7th, 1952, meeting of the 
ARS National Board of Directors, Mr. 
David A. Anderton, engineering editor 
of Aviation Week, was appointed Chairman 
of a newly formed Publicity Committee. 
The duties of this Committee (as broadly 
framed in a study by G. Edward Pendray, 
ARS National Board) are to implement 
the work of the Society in increasing the 
prestige of rocket and jet propulsion 
engineering as a profession; to extend the 
general knowledge of jet propulsion; and 
to stimulate active interest by present 
members and encourage new members to 
contribute to the work of the American 
Rocket Society. 

Mr. Anderton is considered ideally 
qualified for this work on the basis both of 
interest and training. An aeronautical 
engineering graduate of Rensselaer Poly- 
technic Institute in 1941, he spent 5 years 
in aircraft design, for the most part at 
Grumman on jet aircraft. subse- 
quently spent 4 years with General Elec- 
tric’s Project HERMES, on which he was a 
project engineer. Since May of 1950 he 
has been with Aviation Week. Readers 
of this magazine will recall his many con- 
tributions to this publication in the fields 
in which the Society is active. 
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CoNnTRIBUTORS: 


Technical Literature Digest 


H. S. SEIFERT, California Institute of Technology, Associate Editor 


F. C. Gunther 
H. S. Seifert 
F. H. Wright 


D. Altman 
D. I. Baker 
R. B. Canright 


Book Review 


HELIcopTtER ANALysis, by Alexander A. 
Nikolsky, John Wiley & Sons, Inc., New 
York, 1951, 340 pp. $7.50. 


Reviewed by Max L. Wiiuiams, JR. 
California Institute of Technology 


This book is presumably an outgrowth 
‘oe and expanded version of Professor Nikol- 
sky’s original “Notes on Helicopter Design 
Theory” (Princeton University Press, 
_- 1944), and as such contains much of the 
original work. The main difference is that 
ey nearly half of this later version is devoted 
to dynamic stability and control analysis, 
and it is in this respect that the book will 
appeal to the practicing engineer. Longi- 
ps tudinal stability in forward flight is con- 
ez sidered in detail, and the method for ana- 
a _ lyzing the lateral case is indicated. Auto- 
matic stabilization is also considered. For 
the most part the book is directed to stu- 
dents, and the author has made it very use- 
ful for this purpose by employing clear no- 
tation, adequate detail, illustrative exam- 
_ ples, and class problems. 
The bibliography, while not extensive, is 
_ representative, the latest reference extend- 
ing into 1950. 
Letailed performance methods are not 
_ discussed in detail nor is the effect of vari- 
ous rotor configurations on the analysis 
dwelt upon at any great length. The 
< : aut hor has also chosen to exclude any dis- 
cussion of the alternate various types of 
_ power supply, such as the pressure jet or 
 pulsejet. 
While a chapter on modern power-plant 
- possibilities and their effects upon helicop- 
7 ter analysis might have been included, the 
book is thought to fulfill a definite need in 
this rapidly expanding field by introducing 
students to helicopters and in providing 
engineers with a ready reference for sta- 


bility and control work 


Books 

Thermodynamics of Irreversible Proc- 
esses, by S. R. DeGroot, Interscience Pub- 
lishers, Inc., N. Y.; North Holland Publish- 
ing Co., Amsterdam, 1951, xvi + 242 pp. 
$5. 


Bibliography of Books and Published 
Reports on Gas Turbines, Jet Propulsion 
and Rocket Power Plants, by E. F. Fiock 
and C. Halpern, U. 8. Nat. Bur. Stand- 
ards Circ. 509, 1951. 25 cents. 

Automatic Feedback Control, by W. R. 
Ahrendt and J. F. Taplin, McGraw-Hill 
Book Co., Inc., N. Y., 1951, 14 + 412 pp. 

The Magnetic Amplifier, by J. H. Rey- 
ner, Stuart & Richarts, London, 1950, 119 
pp. + 18 plates. 15 shillings. 

The Preparation of Programs for an 
Electronic Digital Computer, by M. V. 
Wilkes, D. J. Wheeler, and S. Gill, Addi- 
son-Wesley Press, Inc., 1951, 167 pp. $15. 


Jet Propulsion Engines 


Measurement of the Efficiency of Turbo- 
Reactor Combustion Chambers by Deter- 
mination of the Carbon Dioxide in the 
Exhaust Gases, by J. Rappeneau, Recherche 
Aéronaut. (Paris), vol. 21, 1951, pp. 19-23. 

Flame Speeds of Methane-Air, Propane- 
Air, and Ethylene-Air Mixtures at Low 
Initial Temperatures, by G. L. Dugger 
and S. Heimel, NACA Tech. Note 2624, 
Feb. 1952. 

Tests on 2 Working Model Ram Jet in 
a Supersonic Wind Tunnel, by J. R. Sing- 
ham, R & M Report 25€8, 1951, 16 pp. 

S.N.E.C.M.A. ‘“‘Escopette” Pulse-Jet 
Unit, unsigned, Aircrafi Eng., vol. 24, 
Jan. 1952, p. 17. 

Escopette Pulsejet, unsigned, Flight, 
vol. 61, Jan. 25, 1952, pp. 101-102. 


Engines 


Rockets for Aircraft Propulsion, by S., 
Allen, The Aeroplane, vol. 81, Dec. 7, 1951, 
pp. 726-727. 

Rockets of the Navy, by W. H. Keigh- 
ley, Ordnance, vol. 36, Jan.—Feb. 1952, pp. 
577-578. 

The Viking Rocket, by W. G. Purdy, 
Aeronaut. Eng. Rev., vol. 11, Jan. 1952, pp. 
16-20. 

Red Rockets, by D. A. Anderton, Avia- 
tion Week, vol. 56, Jan. 14, 1952, pp. 37- 
38. 

Rockets for Bazookas, by J. Geschelin, 
Automotive Ind., vol. 106, Jan. 1, 1952, pp. 
34-36. 


Epitor’s Note: This collection of references is not intended to be comprehensive, but is 
rather a selection of the most significant and stimulating papers which have come to the 
attention of the contributors. The readers will understand that a considerable body of 


literature is unavailable because of security restrictions. 


We invite contributions to this 


department of references which have not come to our attention, as well as comment on how 
the department may better serve its function of providing leads to the jet propulsion appli- 


cations of many diverse fields of knowledge. 
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Heat Transfer and Fluid 
Flow 


Steam Heat Transfer to Small Drops o 
Water, by S. Weinberg, Brit. J. Appl 
Phys., vol. 2, Dec. 1951, pp. 363-366. 

Engineering Aspects of Liquid Metal: 
for Heat Transfer, by Y. Trocki, Nucleon 
ics, vol. 10, Jan. 1952, pp. 28-32. 

The Free Vibration of a Gas Contained 
Within a Spherical Vessel, by H. G. Ferris, 
J. Acoust. Soc. Am., vol. 24, Jan. 1952 
pp. 57-60. 

The Passage of Turbulence Through 
Wire Gauzes, by A. A. Townsend, Q. J 
Mech. and Appl. Math., vol. 4, Sept. 1951, 
pp. 308-320. 

Experimental Study of Isothermal Wake 
Flow Characteristics of Various Flame 
Holder Shapes, by G. G. Younger, NACA 
Research Memo. E61K07, Jan. 28, 1952, 
45 pp. + 6 plates. 

Design of Two-Dimensional Channels 
With Prescribed Velocity Distributions 
Along the Channel Walls. II—Solution by 
Green’s Function, by J. D. Stanitz, NACA 
Tech. Note 2595, Jan. 1952. 

The Propagation in a Compressible 
Fluid of Finite Oblique Disturbances with 
Energy Exchange and Change of State, by 
F. W. Ross, J. Appl. Phys., vol. 22, Dec. 
1951, pp. 1414-1421. 

On Shock-Wave Phenomena; Refrac- 
tion of Shock Waves at a Gaseous Inter- 
face, by H. Polachek and R. J. Seeger, 
Phys. Rev., vol. 84, Dec. 1, 1951, pp. 922 
929. 

A Note on Prandtl’s Formula for the 
Wave-Length of a Supersonic Gas Jet, by 
D. C. Pack, Q. J. Mech. and Appl. Math., 
vol. III, 1950, pp. 173-181. 

On Supersonic Flow of a Two-Dimen- 
sional Jet in Uniform Stream, by S. I. Pai, 
J. Aeronaut. Sci., vol. 19, Jan. 1952, pp. 
61-65. 

The Rapid Discharge of Gases from Ves- 
sels, by J. Kestin and J. S. Glass, Aircraft 
Eng., vol. 23, Oct. 1951, pp. 300-304. 

On the Interaction Between Multiple 
Jets and an Adjacent Surface, by J. Jonas, 
Aeronaut. Eng. Rev., vol. II, Jan. 1952, pp. 
21-25. 

The Boundary Layer in the Converging 
Nozzle of a Swirl Atomizer, by G. I. Tay- 


“lor, Q. J. Mech. and Appl. Math., vol. III, 


1950, pp. 129-139. 

The Application of Boundary-Layer 
Theory to a Swirling Liquid Flow Through 
a Nozzle, by A. M. Binnie and D. P. Har- 
ris, Q. J. Mech. and Appl. Math., vol. III, 
1950, pp. 89-106. 

Recent Investigation of Temperature 
Recovery and Heat Transmission on Cones 
and Cylinders in Axial Flow in the N.O.L. 
Aeroballistics Wind Tunnel, by G. R. 
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Eber, J. Aeronaut. Sci., vol. 19, Jan. 1952, 
pp. 1-6. 


Vapor Phase Oxidation of Hydrocai- 
bons, by F, H. Garner, Trans. Faraday 
Soc., vol. 47, Aug. 1951, pp. 877-899. 

Determination of Aldehydes in Combus- 
tion Products, by H. C. Bailey and J. H. 
Knox, J. Chem. Soc., vol. 75, Oct. 1951, pp. 
2741-2742. 

Material Transfer in a Turbulent Flame, 
hy V. J. Berry, D. M. Mason, and B. H. 
Sage, Jet Propulsion Lab. Progress Report 
20-148, Aug. 24, 1951, 36 pp. 

Application of the Phase Behavior of 
inary Systems to the Improvement. of 
Combustion, by H. B. Palmer, J. O. 
Hirschfelder, and C. A. Boyd, U. Wise. 
Naval Research Lab., Dept. of Chemistry 
(‘M-702, Jan. 2, 1952. 

The Structure of the Boundary Layer in 
he Burning of Solid Fuel, by N. P. Voz- 
nesenkiy and A. B. Chernyshev, Chem. 
ibstracts, vol. 46, Jan. 25, 1952, p. 703f. 

Summary of Combustion Studies for 
Generation of Inert Gas, by J. F. Hill, SAE 
Nat. Aeronaut. Meeting, Los Angeles, Oct. 
3-6, 1951, Preprint 660, 6 pp. 

The Design of a Combustion Tunnel, by 
R. MeCarthy and H. Emmons, Harvard 
Div. Appl. Sci., Combustion Tunnel Lab., 
Jan, 1952. 

Blowoff of Flames from Short Burner 
Ports, by C. W. Wilson and N. J. Hawkins, 
Ind. Eng. Chem., vol. 48, 1951, pp. 2129- 
2135. 

Calculated Data for the Combustion 
with Liquid Oxygen Diluted Alcohols and 
Paraffin in Rocket Motors, by I. C. Hutch- 
eon, R & M 2572, 1951, 13 pp. 

Heats of Combustion of Some Nitro 
Alcohols, by R. M. Currie, C. O. Bennett, 
and D. E. Holcomb, Ind. Eng. Chem., vol. 
44, Feb. 1952, pp. 329-331. 

Combustion Characteristics of Diesel 
Fuels as Measured in a Constant-Volume 
Bomb, by R. W. Hurn and K. J. Hughes, 
Q. Trans. Soc. Automotive Engrs., Jan. 
1952, pp. 24-35. 

Investigation of Hydrocarbon Ignition, 
by C. E. Frank and A. U. Blackham, 
NACA Tech. Note 2549, Jan. 1952, 33 pp. 

New Technique Simplifies Evaluation of 
Precombustion Reactions, by E. B. Rif- 
kin, J. Soc. Automotive Engrs., vol. 60, 
Feb. 1952, pp. 50-55. 

Hydrocarbons (a Symposium), Discus- 
sions Faraday Soc., The Aberdeen Univ- 

ersity Press, Ltd., no. 10, 1951, 339 + vii 
pp. 36 shillings. 

This collection of articles, published as 
part of an irregular series by the same or- 
ganization which produces the Transac- 
tions of the Faraday Society, is the result of 
a symposium on Hydrocarbons held at 
Oxford in April, 1951. It is of wider in- 
terest than the title would indicate, since 
many of the articles include general discus- 
sions of possible combustion mechanisms. 
Two general topics are covered: I. Hy- 
drocarbon Structure and Bond Properties 
(pp. 1-128), twelve articles; and II. Hy- 
drocarbon Reactions, subdivided into A. 
Thermal Reactions (pp. 129-235), eight 
articles; and B. Oxidation Reactions (pp. 
236-338), ten articles. 

Density Measurements in Gaseous Det- 


May-JuNE 1952 


Combustion 


onation Waves, by G. B. Kistiakowsky, J. 
Chem. Phys., vol. 19, Dec., 1951, pp. 1611— 
1612. 

Considérations théorique sur la propa- 
gation des flammes: I. Vitesses. II. 
Domaines de propagabilité, by P. Blanquet 
and M. F. Hoare, J. Chimie phys., vol. 48, 
Sept.—Oct. 1951, pp. 399-406. 

Effect of Preheat on Burning Velocity, 
by G. Dixon-Lewis, Fuel, vol. 30, 1951, pp. 
287-288. 

Propagation d’une déflagration dans un 
tube recouvert d’une pellicule d’huile, by 
R. Loison, Compt. rend., vol. 234, Jan. 28, 
1952, pp. 512-513. 

Théorie simplifiée de la hauteur des 
flammes laminaires de diffusion, by F. 
Guyomard, Compt. rend., vol. 234, Jan. 2, 
1952, pp. 67-69. 

Flammenfortpflanzung in Rohrstrecken. 
II, by H. Behrens, 7. Elektrochem., vol. 55, 
Dec. 1951, pp. 695-696. 

Basic Research on Combustion, Detona- 
tion and Shock Waves at the Laboratoire 
de Recherches Techniques de Saint Louis 
of the Direction des Etudes et Fabrication 
d’Armement, Office of Naval Research, 


London, Tech. Report ONRL-130-51, 
American Embassy, London, England, 
1952. 


Inflammation of Methyl and Ethyl Ni- 
trate Vapors: Characteristics of the Ex- 
plosion Process and Intermediate Reac- 
tions Involving Nitrogen Dioxide, by P. 
Gray and A. D. Yoffe, J. Chem. Soc., 1950, 
pp. 3180-3185. 

Experimental Study of the Radiation of 
Diffusion Flames (in French), by F. Guyo- 
mard, Compt. rend., vol. 233, July 16, 
1951, pp. 237-239. 

A Spectroscopic Investigation into the 
Structure of Diffusion Flames, by H. G. 
Wolfhard and W. G. Parker, Proc. Phys. 
Soc. (London), vol. 65, pt. 1, Jan. 1, 1952, 
pp. 2-19. 

Spectroscopic Studies of the Hydrogen- 
Oxygen Explosion Initiated by Flash Pho- 
tolysis of Nitrogen Dioxide, by R. G. W. 
Norrish and G. Porter, Proc. Roy. Soc., 
Series A, vol. 210, Jan. 22, 1952, pp. 439- 
460. 

Spectroscopic Study of Combustion, by 
G. H. Rothgery and J. T. Grey, USN. 
USAF Proj. Squid Tech. Memo. CAL-34, 
Sept. 28, 1951. 

Flame Radiation Research Joint Com- 
mittee: Reports of 1949 Trials at Ijmui- 
den (Holland), by J. J. Broeze, G. Ribaud, 
and O. A. Saunders, J. Inst. Fuel, vol. 24, 
no. 140, Suppl., 1951. 


Fuels, Propellants, 
and Materials 


Quantitative Evaluation of Rocket Pro- 
pellants, by S. S. Penner, Am. J. Phys., 
vol. 20, Jan. 1952, pp. 26-31. 

High Reaction Energy Substances for 
Jet and Rocket Propulsion, a Review of 
Existing Literature, by L. E. Swearingen 
and A. Veis, USAF-AMC/TR-6633, May 
1951. 

Rocket Propellants, by J. D. Clark, Ord- 
nance, vol. XXXVI, Jan.-Feb. 1952, pp. 
661-663. 

Jet Fuel Quality Requirements, by D. P. 
Barnard, Petroleum Processing, vol. 6, 


1951, pp. 1229-1232. 


The Heats of Combustion of Magnesium 
and Aluminum, by C. E. Holley, Jr., and 
E. J. Huber, Jr., J. Am. Chem. Soc., vol. 
73, Dec. 1951, pp. 5577-5579. 

Ni-Al-Titanium Carbide Cermets, by 
Cc. C. McBride, USAF-AMC TR-6458, 
Jan. 1951. 

Compacts of Pure Tungsten Carbide, by 
A. E. Williams, Metal Treatment, vol. 18, 
Oct., 1951, p. 45. 

New Ordnance Materials, by B. S. 
Mesick, Ordnance, vol. 36, Jan.-Feb. 1952, 
pp. 573-576. 

Elastic and Flow Properties of Dense, 
Pure Oxide Refractories, by J. F. Wygant, 
J. Am. Ceram. Soc., vol. 34, Dec. 1951, pp. 
374-380. 

Apparatus for Low-Temperature Ten- 
sion Tests of Metals, by R. J. Nosborg, 
ASTM Bull., Oct. 1951, pp. 41-45. 

The Behavior of Graphite Under Alter- 
nating Stress, by L. Green, J. Appl. Mech., 
vol. 18, Dec. 1951, pp. 345-349. 

Ceramics in Jet Propulsion, by F. 
Zwicky, Aviation Age, vol. 17, Jan. 1952, 
pp. 31-32. 

Cermets May Answer Jet Designers’ 
Prayers, Pt. I, by W. J. Koshuba and J. A. 
Stavrolakis, Zron Age, vol. 168, Nov. 1951, 
pp. 77-80; Pt. II, ron Age, vol. 168, Dec. 
6, 1951, pp. 154-158. 

A Preliminary Study of the Oxidation 
and Physical Properties of Ti C-Base Cer- 
mets, by C. C. McBride, H. M. Green- 
house, and T. S. Shevlin, J. Am. Ceram. 
Soc., vol. 35, Jan.1952, pp. 28-32. 

New Chromium Carbides Have High 
Temperature and High Corrosion Resist- 
ance, unsigned, Materials and Methods, vol. 
34, Dec. 1951, pp. 69-70. 

Silicon Carbide Refractories Used as 
Alternates for Special Service Alloys, by 
J. L. Everhart, Materials and Methods, vol. 
33, Nov., 1951, pp. 71-75. 

Refractory Bodies Composed of Boron 
and Titanium Carbides Bonded with Met- 
als, by J. A. Nelson, J. Electrochem. Soc., 
vol. 98, Dec. 1951, pp. 465-473. 

Boron Carbide as a Base Material for a 
Cermet, by H. J. Hamjian and W. G. Lid- 
man, J. Am. Ceram. Soc., vol. 35, Feb. 
1952, pp. 44-48. 

Elevated Temperature Properties of Ti- 
tanium Carbide Base Ceramals Containing 
Nickel or Iron, by A. L. Cooper and L, E. 
Colteryahn, Research Memo. 
E51110, Dec. 1951, 47 pp. 

Solar’s Ceramic Coatings, by G. C. 
Close, Aviation Age, vol. 16, Sept. 1951, pp. 
21-22. 

Trends in Powder Metallurgy, by H. W. 
Greenwood, Metallurgia, vol. 44, Dec. 
1951, pp. 291-294. 

Speedier Powder Quality Control Made 
Possible, unsigned, Vetallurgia, vol. 45, 
Jan. 1952, pp. 49-51. 

Are-Cast and Powdered Molybdenum 
Tested, by J. H. Bechtold, Zron Age, vol. 
167, Nov. 22, 1951, pp. 92-95. 

Mechanical Properties of Arc-Cast and 
Powder Metallurgy Molybdenum, by 
J. H. Bechtold, J. Electrochem. Soc., vol. 
98, Dee. 1951, pp. 495-504. 

Cast Molybdenum of High Purity, by 
G. W. P. Rengstorff, J. Metals, vol. 192, 
Feb. 1952, pp. 157 et seq. 

Effect of Strain Hardening and Temper- 
ing on 18-8 Stainless Steels, by P. Bastien, 
Metal Treatment, vol. 51, Oct. 1951. 

Impact Properties of Stainless Steel and 


173 


| 
j 
pl 
ils 
n- 
ed 
is, 
zh j 
| 
ne 
ns 
by 
le 
th F 
| 
Pr, 
| 
| 
h., 
al, 
yp. 
aft I 
as 
AVG 
ng 
LV: 
Il, 
yer 
| 
ure 
nes 
| 


- 9 Per Cent Nickel Steel After Exposure 
Under Stress to Liquid Nitrogen, by T. N. 
Armstrong, ASTM Bull., Oct. 1951, pp. 
85-86. 
Corrosion Aspects of Fusion Welded 
Aireraft High-Strength Aluminum Al- 
loys, by L. W. Smith, Corrosion, vol. 7, 
i Dee. 1951, pp. 423-437. 
; Sprayed Aluminum-Base Coatings for 
Aluminum Alloys, by F. A. Champion, 
Metal Indus., vo.. 79, Oct. 26, 1951, pp. 
355-358. 


Physical-Chemical Topics 


Autoignition by Rapid Compression, by 


C. Livengood and W. A. Leary, 
_ Eng. Chem., vol. 438, Dec. 1951, pp. 2797- 
2805. 


Heat Capacity Lag in Gases, by R. 
Walker, NACA Tech. Note 2537, Nov. 
4 1951, 40 pp. 

The Critical Constants of Polar Gases, 

by J. S. Rowlinson, J. Chem. Phys., vol. 

19, 1951, pp. 831-833. 

Theoretical Foundations of the Kinetics 
of Gas Reactions in Flow Systems (in 
Russian), by M. F. Nagiev, Doklady Akad. 

Nauk USSR, vol. 80, 1951, pp. 397-400. 
_ Spherical Interface. I. Thermody- 
: + namics, by F, P. Buff, J. Chem. Phys., vol. 
—-'19, Dee. 1951, pp. 1591-1594. 

Temperature of Evaporation of Drop- 
7 — lets (in Russian), by N. P. Tverskaya, 
Imest. Akad. Nauk USSR, Ser. Geog. Geo- 

phys., vol. 15, 1951, pp. 74-81. 
wee Experimental Values of the Surface 
| ¥ Tension of Supercooled Water, by P. T. 
- Hacher, NACA Tech. Note 2510, Oct. 
1951, 20 pp. 
7 Evaporation Rate of Liquid Helium, by 
A. Wexler, J. Appl. Phys., vol. 22, Dec. 
1951, pp. 1463-1470. 

Transport Phenomena in Polyatomic 
= by C. S. Wan Chang and G, E. 
Uhlenbeck, Mich. U. Eng. Research Inst., 
Report no. CM-681, July 10, 1951, 46 pp. 
Thermodynamic and Transport Prop- 
erties of Liquids. Problem I. Thermal 
Conductivities of Liquids, by E. Johnson, 
0 Jr., and W. Sheffy, USN fand AF Project 

Squid, TM PR-17, Oct. 1951. 

Absorption of Sound in Fluids, by J. J. 
Markham, R. T. Byer, and R. B. Lindsay, 
Revs. Modern Phys., vol. 23, Oct. 1951, pp. 
353-411. 

Behavior of Gas Discharge Plasma in 
High-Frequency Electromagnetic Fields, 
by L. Goldstein, Elec. Commun., vol. 28, 
: Dec. 1951, pp. 305-321. 

_ Thermodynamics of Carbon Disulfide 
nthesis, Equilibrium in Some Reactions 
Involving Hydrogen Sulfide, Sulfur Diox- 
ide, Carbon Dioxide, and Carbon, by A. J. 
Owen, K. W. Sykes, and D. J. D. Thomas, 


Trans. Faraday Soc., vol. 47, 1951, pp. 
419-428, 


Autoxidation of Hydrazine; Effect of 
Dissolved Metals and Deactivators, by 
L. F. Audrieth and P. H. Mohr, Ind. Eng. 
Chem., vol. 43, 1951, pp. 1774-1779. 

Some Semiempirical Methods for Esti- 
mating the Entropies of Oxides, Silicates, 
and Titanates, and the Prediction of Reac- 
tion Equilibria and Yields, by W. J. Knapp 
and W. D. Van Vorst, J. Am. Ceram. Soc., 
vol. 34, Dec. 1951, pp. 384-387. 

Absolute Density of Nitromethane Va- 
por Above Its Boiling Point; by R. M. Cor- 
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elli, Ann. chim. applicata, vol. 39, 1949, pp. 
591-596. 

The Critical Temperature and Critical 
Pressure of Nitrogen by D. White, A. S. 
Friedman, and H. L. Johnston, J. Am. 
Chem. Soc., vol. 73, Dec. 1951, pp. 5713- 
5715. 

Solubility of Fluorine Gas in Various 
Fluorocarbons and Concentrated Acids, by 
W. B. Burford, III, C. E. Weber, H. C. An- 
derson, M. B. Rogers, et al., Nat. Nuclear 
Energy Ser., Div. VII, 1, Prepn., Proper- 
ties, and Technol. of Fluorine and Org. 
Fluorine Compds., 1951, 188 + 97 pp. 

Thermal Data, Vapor Pressure, and 
Entropy of Chlorine Trifluoride, by J. W. 
Grisard, H. A. Bernhardt, and G. D. Oli- 
ver, J. Am. Chem. Soc., vol. 73, Dec. 1951, 
pp. 5725-5727. 

Preparation and Heats of Combustion 
of Some Amine Nitrates, by T. L. Cottrell 
and J. E. Gill, J. Chem. Soc., 1951, pp. 
1798-1800. 

Empirical Equations for the Calcula- 
tion of the Dissociation Energies of C-H 
and C-C Bonds in the Molecules of Satu- 
rated Hydrocarbons and Free Aliphatic 
Radicals, by V. V. Voevodskii, Doklady 
Akad. Nauk USSR, vol. 79, 1951, pp. 455- 
458. 

Vibrational Intensities. III. Carbon 
Dioxide and Nitrous Oxide, by D. F. Eg- 
gers, Jr., and B. L. Crawford, Jr., J. Chem. 
Phys., vol. 19, Dec. 1951, pp. 1554-1561. 

Heats of Solution and Reaction in 
Liquid Ammonia. VIII, by H. D. Mul- 
der with F. C. Schmidt, J. Am. Chem. Soc,, 
vol. 73, Dec. 1951, pp. 5575-5577. 

Crystal Structure of Nitric Acid, by V. 
Luzzati, Mém. services chim. état (Paris), 
vol. 35, 1950, pp. 7-59. 


Instrumentation and 
Experimental Techniques 


Ambient Temperature Independent 
Thermopiles for Radiation Pyrometry, by 
W. G. Fastie, J. Opt. Soc. Am., vol. 41 
Nov. 1951, pp. 823-829. 

Line Reversal Techniques in the Deter- 
mination of Temperature of Gun Flash and 
Other Rapid Transient Phenomena, by 
J. T. Agnew, Am. Soc. Mech. Eng., Paper 
No. 51—F-6, June 5, 1951. 

Problems of Fluid Measurement, by M. 
Schatzmann, Engineers’ Digest, vol. 13, 
Jan. 1952, pp. 15-16. 

Fundamentals of Orifice Measurement, 
by F. M. Partridge, /nstruments, vol. 25, 
Feb. 1952, pp. 207-209. 

An Induction Flowmeter Design Suita- 
ble for Radioactive Liquids, by W. G 
James, Rev. Sci. Instruments, vol. 22, Dec. 
1951, pp. 989-1002. 

On the Generation of Magnetism by 
Fluid Motion, by H. Bondi and T. Gold, 
Monthly Notices Roy. Astron. Soc., vol. 
110, 1950, pp. 607-611. 

How to Transfer Volatile Liquids Under 
Pressure, by R. J. Sweeney, Chem. Eng., 
vol. 58, Nov. 1951, pp. 158-161. 

Field Inspection of Boiler Tubes with 
Ultrasonic Reflectoscope, by J. A. Tash, 
Trans. Am. Soc. Mech. Engrs., vol. 74, 
Feb. 1952, pp. 201-206. 

Methods and Instruments for Measur- 
ing Smallest Gas Pressures (in German), 
by H. Schwarz, Arch. Tech. Messen, vol. 
188, Sept. 1951, pp. T96-T97. 


Measurement of Instantaneous Vector 
Air Velocity by Hot-Wire Methods, by 


C. E. Pearson, J. Aeronaut. Sci., vol. 19, 
Feb. 1952, pp. 73-82. 

Bibliography on the Measurement of 
Gas Temperature, by P. D. Freeze, Nat. 
Bur. Standards Circular 513, Aug. 20. 
1951 (Supt. of Documents, U.S. Gov. 
Printing Office, Wash., D.C. 15 cents). 

High-Temperature Arcs, by D. W 
Rudorff, Engineers’ Digest, vol. 13, Jan 
1952, pp. 6-8. 

The Panradiometer: an Absolute 
Measuring Instrument for Environmental! 
Radiation, by C. H. Richards, A. M. Stoll, 
and J. D. Hardy, Rev. Sci. Instruments, 
vol. 22, Dec. 1951, pp. 925-933. 

Radiation Measuring Instruments for 
the Infra-Red to Ultra-Violet Waveband, 
by A. C. Menzies, Proc. Inst. Elect. Engrs.., 
Pt. IT, vol. 98, Dec. 1951, pp. 771-780. 

Multiple Shielded High Temperature 
Probes; Comparison of Experimental and 
Calculated Errors, by E. M. Moffett, 
SAE Report T13, Jan. 1952, 20 pp. + 3 
plates. 

A New Method for Determining the 
Static Temperature of High-Velocity 
Gas Streams, by J. A. Clark, Trans. Am. 
Soc. Mech. Engrs., vol. 74, Feb. 1952, 
pp. 219-228. 

A Solar Furnace of 120 Inches Diameter 
for Studying Highly Refractory Sub- 
stances, by W. M. Conn, Rev. Sez. Instru- 
ments, vol. 23, Dec. 1951, pp. 945-951. 

Calorimeter for the Determination of 
the True Heat Capacity and of Latent 
Heats at High Temperatures, by V. A. 
Sokolov, Zhur. Tekh. Fiz., vol. 18, 1948, 
pp. 813-822. 

New Precision Calorimeter for Deter- 
mining Heats of Combustion, by A. Mag- 
nus and F. Becker. Z. physik. Chem., 
vol. 196, 1951, pp. 378-398. 

Semiconductors Capable of Withstand- 
ing High Temperatures (Refractory Ther- 
mistors), by N’Guyen Thien Chi and J 
Suchet, Ann. Radio-électricité, vol. 6, 
April 1951, pp. 99-105. 

Accurate Time for Scientific Observa- 
tions, by F. E. Fowler, Electronics, vol. 25, 
Jan. 1952, pp. 98-101. 

The Performance Equations of the Eddy 
Current Gyroscope, by R. Clay and R. 
Harrison, USN-BuOrd /NAVORD-R-1397, 
Oct. 26, 1951. 

The Impact of Servomechanism Tech- 
niques on ition, by A. Porter, 

ASME Paper No. 51—SA-22, March 20, 
1951, 22 pp. 


Flight, Ballistics, and 
Vehicle Development 
Life on Mars in View of Physiological 
Principles, by H. Strughold, CADO Tech. 


Data Digest, vol. 16, Nov. 1951. 
Dirigibility of Rockets, by H. Barten- 


“bach, J. Space Flight, vol. 3, Dec. 1951, 


pp. 1-2. 

The Artificial Satellite, by E. Burgess, 
The Engineer, vol. 192, Oct. 12, 1951, 
pp. 456-458. 

Advance Man to the Moon, unsigned, 
Aero Digest, vol. 63, Oct. 1951, pp. 25-32. 

Power Sources for Orbital Rockets, by 
L. J. Grant, J. Space Flight, vol. 3, Nov. 
1951, pp. 1-3. 

Preview of Space Flight, by L. 8. Black, 
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Aero Digest, vol. 63, Oct. 1951, pp. 17-24. 
The Brainwork Is Done, by R. Me- 

Larren, Aero Digest, vol. 63, Oct. 1951, 

pp. 34, 36, 38, 42, 44, 46, 50, 52, 54, 

5&8, 60, 63, 64. 

Shapes of Craters Formed in Plaster of 
Paris by Ultra-Speed Pellets, by J. S. 
Rhinehart and W. C. White, Am. J. 
Phys., vol. 20, Jan. 1952, pp. 14-18. 

Thermal Stresses in a Rotating Elastic 
Solid of Revolution, by C. E. Iliffe, 
ner, VOl. 178, Dec. 28, 1951, pp. 835-836. 


Slingshots or Precision Weapons, un- 
signed, Interavia, vol. VI, Oct. 1951, pp. 
541-547. 

Rocket Away!, unsigned, Jntleravia, 


vol. VI, Oct. 1951, pp. 556-557. 
Guided Anti-Aircraft Missiles, unsigned 
Interavia, vol. VI, Oct. 1951, pp. 554-555. 
The Pilotless Aircraft Type Missile, by 


SE. Weaver, SAE Report 731, Jan. 
1952, 9 pp. + 1 plate. 

USAF’s Guided Missile Test Center, 
Pt. Il, by R. E. Stockwell and R. Haw- 
thorne, Aviation Age, Dec. 1951, pp. 
40-43. 

What We Have Learned from V-2 


lirings, unsigned, Aviation Week, Nov. 26, 
1951, vol. 55, pp. 23-30. 

Some General Aspects of 
Testing, by J. H. Zabriskie, SAE 
733, Jan. 1952, 6 pp. 

Bumblebee Program, by J. E. Arnold, 
Aviation Age, vol. 17, Jan. 1952, pp. 31-35. 

Peenemunde, 1951, unsigned, Aviation 
Age, vol. 117, Feb. 1952, p. 33. 

Supersonic Flight from the Pilot’s View, 
by W. B. Bridgeman, Aeronaut. Eng. Rev., 
vol. 11, Feb. 1952, pp. 24-28. 

Atomic-Powered Aircraft, unsigned, 
Flight, vol. 60, Dec. 21, 1951, p. 780. 

P-W to Build Atomic Aircraft Engine, 
unsigned, Aviation Week, vol. 55, Dec. 
24, 1951, p. 18. 


Missile 
Report 


Other Classifications 
Meteorology, Astrophysics, Ete. ) 


The Scintillation of Starlight, by A. H. 
Mikeseli, J. Opt. Soc. Am., vol. 41, Oct. 
1951, pp. 689-695. 

Trends in U. 8. Air Force Research and 
Development, by D. L. Putt, Tech. Data 
Digest, vol. 17, Jan. 1952, pp. 7-12. 

Post-War Developments in Operational 
Research, by A. W. Swan, Engineer, vol. 
173, Dec. 21, 1951, pp. 812-814. 

High Altitude Research at the 
Physics Laboratory, by L. W. 
Bumblebee Report No. 153, APL, 


Applied 
Fraser, 
Johns 


Hopkins U., May 1951. 
Guided Missiles, a New Field for the 
Radio Engineer, by 8S. Ramo, Proc. Inst. 


Radio Engrs., vol. 40, Jan. 1952, p. 3. 

Methods Developed at the Naval 
Ordnance Test Station for Rapid Calcula- 
tion of Air Density, by N. R. Williams 
and K. L. Coulson, USN-BuOrd/- 
NAVORD-R-1910, Aug. 28, 1951. 

The Temperature of the Upper Atmos- 
phere, by D. R. Bates, Proc. Phys. Soc., 
(London), vol. B, Sept. 1951, pp. 895-821. 

Interpretation of the Polarization of 
Light Reflected by the Different Regions 
of the Surface of the Planet Mars (in 
French), by A. Dollfus, Compt. rend., 
vol. 233, Aug. 6, 1951, pp. 467-469. 

Meeting of the Space Medicine Associa- 
tion, unsigned, J. Aviation Med., vol. 23 
Dec. 1951, pp. 544-547, 
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plots graphs 


10 times faster 


than by hand 


operates from IBM or keyboard 
digital input 


A new problem has developed in modern science 
and engineering—the problem of manual plotting. 


The slowness, fatigue, human errors of 
manual plotting frequently cause serious time 
lags, increase expenses and overhead. 


The Teleplotter solves that problem. Plotting 
40 points per minute, it is faster, cheaper, 
more accurate than manual plotting. 


How it operates: 

Data is introduced into the 

Teleplotter by IBM equipment or 

by the Teleplotter keyboard. 

The Teleplotter Reading-Head, housing 
photo-electric “eyes,” travels over the 
graph paper, counts grid lines and 
spaces, prints a symbol when its counts 
coincide with the IBM or keyboard data. 


” 


Write for brochure, “Automatic Data Analysis, 
giving full details on the Teleplotter and 

other Telecomputing Instruments. Coupon on 
the right is for your convenience. 


4 
Vital data: 
Plotting speed—40 points per minute, 
3 approximately 10 times faster Ke 
than manual plotting 
st, Plotting area—26x55inchesorless 
2 Plotting accuracy — +.25 mm(paper 
= stretch does not affect accuracy) 24 
= Scale factors —five factors, with t's 
3 independent selection for aA 
each axis tsy 
“a Plots X-Y data on linear or “a 
logarithmic paper 
pes Time required to set up Teleplotter bes 
for operation — negligible 


lew 


2 
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TELECOMPUTING 


CORPORATION, BURBANK, a 


“Turn hours into minutes with Automatic Data Analysis” 


structural test data; fractionating 


Applications: . 

Experimental uses 
Plots wind tunnei data; static and a 
tower design data, etc. Z 
Plots results of manual calculations “di 
and results of automatic electronic 24 
calculating machines to detect 23 


random errors. 


Theoretical uses 
Plots theoretical curves and results 
of theoretical calculations. 

General uses 

Plots insurance tables; stock market, 
statistical and financial records; 
production control data; results 

of engineering studies and analyses; 
management control data. 

Plots results of continuous 

process operations for supervision 
and control purposes. 


Mr. John H. Weaver, Engineering , Dept. —RT-5 
Dear Sir: Please send me ‘‘Automatic Data 
Analysis.” 


Name 


Company 


Street address 


City and State 


A new achievement in aufomatic data analysis a 
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ere is a better lining material 
for uncooled rocket motors 


Exceptional resistance to the high 
temperature and severe erosive 
conditions developed in uncooled 
rocket motors is offered by spe- 
cial silicon carbide liners and 
throats. In comparative tests they 
have outlasted metals by several 


times. 

Information is available upon 
request to Dept. RS-52, Refractor- 
ies Div., The Carborundum Co., 


Perth Amboy, N. 
CARBORUNDUM | 


Trade Mark 


“‘Carborundum” is a registered trademark indicating manufacture by The Carborundum Company. 


GROVE 
SELF-ACTUATED Get the Facts about 


CONTROLS AND REGULATORS 


your Future 
for 
pressure liquids and gases wit h 
7 


= 


Pressure Reducing 
Back Pressure Discover the 
Relief and greater opportuni- 
mea Manual Valves ties offered engi- 
neers by the 
greatest diversity of 
projects of any air- 
craft company in the 
East! Write today 
for fact-packed 
brochure. 


fe Special light weight units 
for air-borne service. 


Special materials for corrosive. 


or low-temperature applications. SEND FOR ENGINEERING BROCHURE 


THE GLENN L. MARTIN COMPANY 
GROVE 6 NTR LS if Personnel Dept. « Section A « Baltimore 3, Md. LY 
0 . s Please send me your brochure describing engineering opportuni- X 
ties at Martin. 
Gus 
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Ulustration is artist's conception of Air 
Force B-61 Matador pilotless bomber. 
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GUIDED MISSILES 


The Air Force Missile Test Center, used by all our 
military services as a long-range proving ground, 
stretches thousands of miles from Florida, out 
over the Bahamas, into the South Atlantic. 


— bomber roars away from its launching 
stand, picks up speed, zooms into the blue. Set- 
ting its course for a far-off target in the ocean, it 
rockets over a chain of tiny islands where men and 
machines check its flight, its behavior, the operation of 
its guidance and control systems. It’s a vital part of our 
air power of the future—aeronautical research and de- 
velopment laying the foundation for continued U. S. 
air supremacy! 


Operated by the USAF’s Air Research and Develop- 
ment Command, the Missile Test Center is geared up 
to test the wide variety of missiles, rockets and pilot- 
less aircraft vital to modern air power. It reached its 
full stature with the recent completion of down-range 
observation stations. And the dramatic B-61 pilotless 
bomber, the Matador, designed and produced by 
Martin as part of its diversified missiles program, was 
the first to use the completed range. THE Gienn L. 
Martin Company, Baltimore 3, Maryland. | 


AIRCRAFT 


Builders of Dependable Aircraft Since 1909 


DEVELOPERS AND MANUFACTURERS OF: Navy P5M-1 Marlin Navy Viking high-altitude research rockets * Air Force XB-51 
seaplanes « Air Force B-57A Canberra night intruder bombers « developmental tactical bomber * Martin airliners * Guided 


Air Force B-61 Matador pilotless bombers * Navy P4M-1 missiles © Electronic fire control & radar systems « LEADERS IN 
Mercator patrol planes « Navy KDM-1 Plover target drones « Building Air Power to Guard the Peace, Air Transport to Serve It. 
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2° CONTROL @ INDUSTRIAL SOUND O 
q 
PHYSICISTS AND SENIOR 3c 
4 
RESEARCH ENGINEERS "rocket 
gtlemced by 
4 
mutffi 
POSITIONS NOW OPEN mufflers 
c 
,@ @ In the laboratory, in the test cell, on the va 
eeee 7 airfield, hundreds of INDUSTRIAL SOUND ms 
q~ CONTROL installations are subduing the 
noise, heat and gas velocities generated dur- 
Senior Engineers and Physicists having }« ing toting of the tig jet. - 
d d Whatever your noise problem, ISC’s skilled > 
outstanding academic bac groun an | engineering, design, and installation “know r 
of: } fe) how” — gained through years of practical 
experience in the fie of: » 7 ea | qu experience — can help you, and is at your “ 
} instant service, 
4° Po) We welcome the challenge of the unusual 
@ Microwave Techniques | Zz problem. For full information c 
> z 
@ Moving Target Indication Yo 
@ Servomechanisms Cementation (Muffelite) Limited, London 
e Applied Physics : Les Travaux Souterrains, Paris “ 
@ Gyroscopic Equipment 
4 
@ Computers }> 45 Granby Street, Hartford, Conn. 
@ Pulse Techniques 2119 So. Sepulveda Bivd., Los Angeles, Calif. 
@ Radar € 
@ Fire Control SOUND CONTROL 
@ Circuit Analysis 
@ Autopilot Design 
@ Applied Mathematics ; | CURRAN ENGINEERING CO. 
@ Electronic Subminiaturization © 
@ Instrument Design { 
4 
@ Automatic Production Equipment Manufacturer of 
@ Test Equipment } 
@ Electronic Design MECHANICAL COMPONENTS “ 
@ Flight Test Instrumentation from 
METALS, CERAMICS, AND PHENOLICS 
4 
are offered excellent working conditions | ueennts 
and opportunities foradvancementinour | Consultants and Specialists of 
4 
Aerophysics Laboratory. Salaries are | ROCKET IGNITER ASSEMBLIES 
commensurate with ability, experience | and 
and background. Send information as | LONGITUDINAL SHAPED CHARGE 
4 
to = education, —_— and work | CUTTERS 
preference to: : 
NORTH AMERICAN AVIATION, INC. § CENCO' Process for 
HIGH TEMPERATURE-HIGH DIELECTRIC 
___ Aerophysics Laboratories INSULATING OF METALLIC ASSEMBLIES 
Box No. S-4 
12214 South Lakewood Blvd. 
= Downey, California 7614 Santa Monica Blvd. 
! Los Angeles 46 California. 
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* New Guided Missiles Count On 
Greer Accumulators to Supply Power 
Actuate Controls and Mechanisms 


High Pressure Sphericals 


There‘s a Greer guided missile accumulator for almost every 
requirement, or Greer will build to your exact specifications. 


a It's got a target by the tail! 


HIS GUIDED MISSILE may go through some mighty tight turns 
= to stay on target... and that means G loads go up, fast. But 
the Greer Accumulator in the circuit will take it. 


Field Office: 298 Commercial Bldg., Dayton,O. © Thomson Engineering Service, 708 Hemphill St.. Fort Worth 4,Tex, * Harold E. Webb, 918 N. Kenilworth Ave., Glendale 2, Calif 


SPECIAL PRODUCTS DIVISION - GREER HYDRAULICS INC. 


For Greer Accumulators are designed to operate under the high 
G loads caused by tight maneuvers and high accelerations. 
Their volumetric efficiency is extremely high. They are prac- 
tically impervious to temperature; function under operating 
pressures as high as 6000 psi; can vary in capacity from 2 
cubic inches to 25 gallons. 

Uses for Greer Accumulators in a guided missile are many, 
varied, and important. They actuate Servo controls and mech- 
anisms, function as a source of instantaneous energy, power 
supply for primary and secondary circuits, a pressurized reser- 
voir, pressure-transfer barrier,and eliminate pressure pulsations. 
These compact. lightweight power sources can be designed to 
fit almost any given space. They are unusually flexible in shape, 
capacity and pressure, and can be planned as an integral part 
of the missile itself. 


Developing accumulators for guided missiles is just one interest 
of the Greer Special Products Division. If you are confronted 
with a problem relating to the field of aviation, call Greer en- 
gineers for a seca There is no cost, no obligation. 


GUIDED MISSILE 


454 Eighteenth Street, Brooklyn 15, New York 
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of 
‘noe Jet Propulsion Laboratory 


. . has several openings in the following engineering fields: 


angular 


accelerometers 


Preliminary Design 


Engineering Analysis “~~ 
Aerodynamics 
Aerodynamic Heating 
Heat Transfer 


Thermodynamics 


Dynamics 
Structural Studies 
Statistical Studies 
+ Structural Design 
Propulsion System Development 
Test Engineers—Wind Tunnel 
Experience in supersonic 
wind-tunnel tests desirable 
Excellent opportunity exists to learn all phases related to 


missiles, rockets and wind tunnel test operations. Apply 
giving details pertaining to academic background and work 
Please request 


experience to 


LABORATORY - 


4800 Oak Grove Drive : 
Pasadena 2, California 


Z) 


LABORATORIES 
Angeles 64, Calif... 


A... accelerometers 


are made by Statham Lab- 
oratories for measurement in 
ranges as low as + 1.5 radians 
per second per second. 


The design permits close 
balance against linear 
acceleration effects with a 
high degree of mechanical 
shock resistance and leads to 
a damping characteristic 
relatively insensitive 

to temperature. 

The transducer element, an 
unbonded strain gage bridge, 
provides an electrical output 
proportional to applied 
angular acceleration for 
recording or telemetering in 
conventional a.c. or 


d.c. circuits. 


SCINTILLA MAGNETO DIVISION 


BENDIX AVIATION CORPORATION 


iretanerers of Ignition Systems for Jet, 


@ Piston Power Plants, Fuel Injection Equipment for Railway, 


Marine and Industrial Diesel Engines, Electrical Connectors, 
Ignition Analyzers, Moldings and Special Components. 


Turbine and | 
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fighter 


with long-range 


radar eyes 


Designed and built for Navy carriers, 
the Douglas F3D Skyknight provides our 
fleets with round-the-clock protection. 
Attack, patrol, reconnaissance, or escort, 
Skyknight can handle them all. 

Aided by its radar eyes, the Skyknight 
can search out distant targets 24 hours 
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the Douglas Skyknight 


a day. The pilot of this unique two-man. 
twin-jet, long-range fighter—guided by 
his radar operator—comes in on targets 
with split-hair accuracy .. . to hit with 
both rockets and bullets. And although 
Skyknight approaches sonic speeds. its 
hydraulic flaps can slow it down for 


combat maneuvers or carrier landings. 

The carrier-based F3D Skyknight. 
now in volume production, is typical of 
Douglas leadership in aviation. Planes 
that can be mass-produced to fly further 
and faster with a bigger payload is the 
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The United States Army Field Forces today 
is training its first Guided Missiles Group 
with Fairchild Missiles. In firing these ad- 
vanced type anti-aircraft missiles, the Army 
Field Forces is preparing now for the day 
when missile batteries will defend cities and 
vital military installations. 


Firing on the desert missile range at 
Ft. Bliss, Texas, officers and enlisted cadres 
are learning the skills and techniques neces- 
sary for the tactical application of these 


Fal AND AIRPLANE CORPORATION 


Fairchild Aircraf , Hagerstown, Maryland @ 


ADVANCED TRAINER... 


FOR THE FIRST GUIDED MISSILES GROUP, U.S. ARMY 


The Fairchild Plant at Wyandanch, L. 1.—The World's First Privately Built Missile Plant 


new weapons under conditions similar to 
those in actual combat. 

Fairchild’s Guided Missiles Division also is 
providing similar anti-aircraft missiles for the 
United States Navy and the United States 
Air Force. Its advanced engineering and 
technical facilities are being devoted to the 
design and development of new missiles and 
improved versions of current missiles to pro- 
vide our Armed Forces with the latest and 
best possible weapons. rs 


AIRCHILD Divition 


WYANDANCH, L.I., N.Y. 


Fairchild Engine and Stratos Divisions, Farmingdale, Long Island, New York 
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IN SIZE... 
GIANTS 
IN PERFORMANCE! 


ACTUAL 
SIZE 


TOP BOTTOM 
VIEW VIEW 


LOW LOSS FACTOR OF .015 at 

1 MC. : 

e INSULATION RESISTANCE (min.) 
50,000 Megohms 

e ARC RESISTANCE, ASTM SECS. 
250+ 

© OPERATING TEMPERATURE 
135°C (limitation due to contact 
metal) 

© Q-value at 40°C—50% RH—not 
less than 1000 

® Q-value at 40°C—90% RH min.— 

not less than 333 


WRITE FOR DATA SHEETS 


Mycalex 8-Pin Sub-Miniature Tube Sockets are 
fully described in the new Data Sheets. Other 
tal are available on Mycalex Insulation 


for every electronic or electrical application. 


SINCE 1919 
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MYCALEX CORPORATION OF AMERICA 


Executive Offices: 30 ROCKEFELLER PLAZA, NEW YORK 20 —Piant & General Offices: CLIFTON, NJ. 


,8-PIN 
‘Y SUB- 
MINIATURE 
TUBE 
3 


OCKETS 


MYCALEX 410 
SOCKET BODY 
CONTACTS 


MATL SEE TABLE 


130 + 003 OD. 
CENTER SHIELO 
BRASS CAD PLATE 


New MYCALEX 410 Sub-Miniature Tube Sockets are designed 
for use in electronic and electrical equipment where space is at 
a premium. Because they are extremely compact, these sockets 
offer a ready solution to numerous design problems involving 
spatial limitations. Installation is simple, mounting being ac- 
complished without screws or rivets in shaped chassis holes. 


Improved electrical performance and greater mechanical pro- 
tection for the tube than are available with ordinary insulat- 
ing materials are afforded by this socket through the use of 
MYCALEX 410 glass-bonded mica. MYCALEX 410 is rated 
Grade L-4B insulation under N.M.E.S. JAN-I-10. It offers 
superior electri- 
cal and mechani- 
cal properties in 
combination 
with practical 
cost per unit. 


MYCALEX TUBE SOCKET CORPORATION 
Under Exclusive License of 
MYCALEX CORPORATION OF AMERICA 


30 ROCKEFELLER PLAZA, NEW YORK 20, N.Y, 


Owners of ‘MYCALEX’ Patents and Trade-Marks 
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ANHYDROUS AMMONIA ° 
ITROSYL CHLORIDE 


n Products 
NITROGEN SOLUTIONS oN 


Representative Solvay N 
NITRATE SODIUM NITRATE 


ROCKET DEVELOPERS * 
4 
SUBJECT: SHIPMENT OF NITROGEN TROXIDE 
Nitrogen Tetroxide is no longer @ 
quid propellants, you 
xiadizer- Solvay 15 
9 ship it in 1.C.C.- approved cyl- 
ners. Free technical 
If you would Like to test Nitrogen 
> 
Tetroxide, we will gladly furnish cost and 
shipping information. 
HE SOLV AY PR —— 
MMONIUM 
METHANOL 


